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ABSTRACT 


Volume  1  of  this  report  discusses  the  radiative  and  collisional 
processes  which  determine  the  recombination  rate  in  a  plasma.  A  gen¬ 
eral  expression  is  derived  for  the  recombination  rate,  and  numerical 
values  are  given  for  a  hydrogen  plasma  in  the  temperature  range  250 °K 
to  64,000®K  and  for  all  electron  densities.  These  results  are  applied  to 
an  expanding  hydrogen  plasma  and  to  the  expanding  debris  from  a  nuclear 
explosion,  and  indicate  that  in  an  explosion  occurring  in  a  vacuum,  sig¬ 
nificant  recombination  may  take  place  at  times  later  than  a  millisecond 
after  the  detonation. 

Volume  II  contains  seven  papers  dealing  with  the  physics  of  recombi¬ 
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a  number  of  atoms  and  ions  obtained  from  a  theory  of  atomic  screening 
constants. 
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THE  STATIONARY  PLASMA* 


2.1  The  Basic  Equations 

In  order  to  Investigate  the  non ‘equilibrium  properties  of  a 
partially  ionized  plasma,  it  is  necessary  to  consider  the  electronic  and 
atomic  processes  which  affect  the  populations  of  the  individual  atomic 
energy  levels.  A  particular  level  may  be  populated  by  spontaneous  and 
stimulated  transitions  from  other  discrete  levels,  by  Impact  induced 
transitions  from  bound  states,  by  radiative  recombination  of  free  elec¬ 
trons  and  by  three-body  recombination  processes  involving  free  electrons. 
It  may  be  depopulated  by  the  inverse  processes  of  radiative  transitions 
out  of  the  level,  impact  induced  transitions  out  of  the  level,  photo¬ 
ionization  and  impact  induced  ionization. 


The  increase  in  population  per  second  of  a  particular  level  (p)  is 
given  by  the  differences  between  the  rates  of  population  and  depopula¬ 
tion.  The  equations  for  the  populations  in  each  of  the  levels  will 
have  the  form 


J  ^  P 


^jp  \  ■  %P 


(2-1) 


In  this  equation  n(p)  is  the  population  of  the  p^^’  level.  Tlie  expres¬ 
sion  n(J)  appearing  in  the  summation  represents  the  rate  at  which 
atoms  are  entering  the  p^^  level  from  the  level.  Tlic  term  b^ 
represents  the  rate  of  capture  of  free  electrons  into  the  p**^'  level 


and  the  expression  a  n(p)  represents  the  rate  at  which  atoms  are 

PP 


*Thi8  chapter  is  repeated  In  Volume  II  to  make  It  available  in  an  un¬ 
classified  document. 
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leaving  the  level  for  other  levels  and  the  ionized  state.  The  a's 
and  b's  include  the  effects  of  both  radiative  and  collision  processes, 
and  they  are  functions  of  the  density  and  temperature  of  the  free  elec¬ 
trons.  Explicit  expressions  for  the  coefficients  are  developed  in 
Appendix  A. 

2.2  Comparison  of  the  Radiative  and  Collisional  Processes 

Tables  giving  the  coefficients  for  the  radiative  and  collisional 
processes  are  collected  in  Appendix  A.  The  tables  for  hydrogen  are 
essentially  complete.  For  other  elements  the  tables  represent  a  start 
on  the  cookie te  assembly  of  data  needed  for  estimating  recombination  in 
complex  plasmas.  For  the  present  they  serve  as  an  indication  of  the 
compilation  and  are  useful  in  showing  the  pattern  of  variation  of  the 
coefficients. 

An  examination  of  the  tables  shows  that  the  spontaneous  transition 
probabilities  decrease  rapidly  with  increasing  excitation,  in  contrast 
to  the  collisional  excitation  and  ionization  probabilities,  which  in¬ 
crease  rapidly.  Thus,  from  the  tables  for  hydrogen  in  Appendix  A,  it 

12  -3 

may  be  observed  that  even  if  the  electron  density  is  as  low  as  10  cm 
and  the  electron  temperature  is  as  low  as  1000“K,  electron  impact  ioniza¬ 
tion  of  level  n  =  10  is  more  rapid  than  tlic  loss  by  spontaneous  transi- 

20  -3 

tion.  For  an  electron  density  of  10  cm  and  an  electron  temperature 
of  1000*’K,  electron  impact  ionization  of  the  first  excited  level  is  more 


2 


rapid  than  Che  loss  by  a  spontaneous  transition.  Electronic  excita¬ 
tion  and  de-excitation  processes  are  still  more  efficient  so  chat  a 
Boltzmann  distribution  of  Che  highly  excited  states  is  to  be  expected 
except  at  low  electron  densities  and  Cenqperatures. 

2.3  Previous  Calculations  of  Recombination  Rates 

(2-2) 

Baker  and  Menzel  have  applied  equations  of  the  form  of 

Equation  (2-1)  Co  an  investigation  of  the  relative  intensities  of 
spectral  lines  in  gaseous  nebulae  in  which  the  ground  state  of  atomic 
hydrogen  is  ionized  by  stellar  radiation.  In  these  circumstances,  a 
steady  state  prevails  so  that 

—  n(p)  =  0 

and  Equation  (2-1)  reduces  to  a  set  of  algebraic  equations,  the  solu¬ 
tion  of  which  is  equivalent  to  tlie  inversion  of  a  matrix.  Bakarand 
Menzel  do  not  distinguish  between  the  individual  sub-levels  corres¬ 
ponding  Co  a  given  principal  quantum  number  and  consequently  assume 
impliciCy  chat 

n(n,i)  =  )  n(n,i)  (2-2) 

n 


3 


(2-3) 

Burgess  '  has  repeated  the  calculations  taking  account  of  the 
different  sub-levels  and  he  gives  the  equilibrium  populations  for  electron 
temperatures  of  10,000°K  and  20,000‘^K.  He  finds  that  the  ^-levels  are 
far  from  being  populated  according  to  their  statistical  weights  and  that 
for  an  optically  thin  plasma,  the  p-levels  (i*l)  are  much  under-populated 
due  to  (he  large  probability  of  np  -  Is  transitions. 

In  both  Investigations,  all  collision  processes  are  Ignored.  Al¬ 
though  this  Is  a  permissible  approximation  for  the  low-lying  levels  of 
hydrogen  when  the  electron  density  and  the  electron  temperature  are  not 
high  (as  Is  the  case  in  must  gaseous  nebulae),  It  will  ultimately  be 
Invalid  as  n  Increases. 

Since  collision  processes  which  redistribute  the  sub-levels  corres¬ 
ponding  to  a  given  n,  such  as 

e  +  H(2p)  -♦  e  +  H(2s)  (2-3) 

+  H(2p)  -  h'*'  +  H(2s)  (2-4) 

are  very  efficient,  the  effect  of  Including  collision  processes  will 
tend  to  ensure  Equation  (2-2).  The  more  refined  treatment  of  Burgess  is 
therefore  actually  less  accurate  than  that  of  Baker  and  Menzel  when  the 
electron  temperature  Is  high  and  when  the  principal  quantum  number  n 
Is  large. 

Glovanelll^ '  has  included  electron  collision  processes  In  his 
treatment,  which  again  assumes  statistical  equilibrium.  In  his  first 
paper  In  which  he  assumes  Equation  (2-2,))ic  shows  that  Saha's  ionization 


4 


formula  applies  to  the  populations  of  the  high-lying  levels,  provided  the 
electron  teiiq>erature  is  employed.  This  occurs  because  the  most  efficient 
depopulating  mechanism  is  electron  impact  ionization  and  the  most  effi¬ 
cient  populating  mechanism  is  the  inverse  process  of  three-body  recom¬ 
bination  and  these  processes  are  very  rapid.  The  same  situation  must 
apply  even  in  non-equilibrium  cases. 

(2-5) 

In  a  second  paper  Giovanelli  discusses  the  populations  of  the 

low-lying  levels.  The  study  contains  an  inconsistency,  which  may  be 
serious,  in  that  while  he  distinguishes  between  the  individual  sub- 
levels,  he  ignores  redistribution  processes  such  as  Equation  (2-3). 

His  results  are  therefore  invalid  at  high  electron  densities  and  tempera¬ 
tures.  We  note  in  passing  that  his  results  are  also  invalid  at  low  elec¬ 
tron  densities  and  temperatures  in  that  he  ignores  the  decay  of  the 
metastable  2s  state  by  the  two-quantum  process 

H(2s)  -  Hds)  +  hv^  +  hv^  (2-5) 

Discussions  of  the  hydrogen  plasma  in  non-equilibrium  situations 
have  been  given  by  D'Angelo'  and  by  Bates  and  Kingston^  In 

both  investigations,  heavy  particle  collisions  are  ignored  and  it  is 
assumed  that  the  plasma  is  optically  thin.  Both  investigations  assume 
Equation  (2-2)  and  ignore  processes  such  as  Equation  (2-3).  This  is 
equivalent  to  supposing  the  Equation  (2-2)  is  very  fast. 

D'Angelo  obtains  an  approximate  solution  of  Equation  (2-1)  by  attemp¬ 
ting  to  follow  through  in  detail  the  reaction  paths  of  any  individual 
electron,  making  the  assuiqption  that  recombination  is  effective  only  for 
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those  paths  which  terminate  in  the  ground  state.  He  ignores  the  effect 

of  electron  impact  excitation  and  de-excitation.  He  gives  results  for 

three  electron  temperatures,  lOOO^K,  3000°K  and  10,000°K  and  for  electron 

12  13  -3 

densities  in  the  range  from  10  to  10  cm  only.  His  investigation 
reveals  the  importance  of  three-body  recombination  processes  (as  did  the 
earlier  work  of  Giovanelli)  but  his  method  of  solution  is  laborious  to 
apply  and  the  accuracy  of  the  numerical  values  of  the  effective  recombina¬ 
tion  coefficients  is  uncertain. 

The  method  employed  by  Bates  and  Kingston  is  simple  and  more  accurate. 
They  argue  that  the  rates  at  which  atoms  in  a  particular  excited  state 
are  produced  or  destroyed  are  very  much  greater  than  Che  rate  at  which 
Che  very  small  number  density  of  atoms  in  Chat  state  changes.  It  follows 
Chat  for  an  optically  Chin  plasma  except  for  Che  ground  state 

^  n(p)  =0  (p  =  2,  3,  ...)  (2-6) 

so  Chat  recombination  is  effective  only  into  Che  ground  state.  Then 
writing 

S“<‘>  (2-7) 

(where  n(l)  signifies  the  ground  state) 

defines  Che  effective  recombination  coefficient  or  ,.(n  ,T  ),  which  will 

eff  e’  e  ’ 

depend  upon  Che  electron  density  as  well  as  Che  electron  temperature. 

Their  method  reduces  the  sec  of  coupled  differential  equations,  Equation 
(2-1),  to  a  set  of  coupled  algebraic  equations.  The  number  of  levels  re¬ 
quiring  individual  examination  is  restricted  since  Che  populations 
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of  the  highly  excited  levels  are  given  by  the  Saha  formula,  the  level 
at  which  this  becomes  true  being  provided  by  the  calculations. 

2. 4  The  Number  of  Discrete  Levels 

In  none  of  the  investigations  we  have  described  has  the  question  of 
the  number  of  levels  that  should  be  included  in  the  equations  been  exam¬ 
ined.  It  is  clearly  not  infinite  because  ultimately  the  discrete  levels 
will  overlap  the  continuum. 

Ivanov-Kholodnyi,  Nikol'skii  and  Gulyaev '  have  discussed  the 
lowering  of  the  continuum  due  to  pre-ionization.  The  treatment  is  not 
entirely  satisfactory  but  should  suffice  for  our  purposes.  They  find 
that  the  number  n^  of  levels  of  different  principal  quantum  number  which 
remain  discrete  is  given  by 

21.65  -  £n(n  ) 

fn(n^  +  1)  - g - S_  (2-8) 

Values  of  n^  have  been  computed  from  Equation  (2-8)  for  a  wide  range 
of  electron  densities  and  they  are  presented  in  Figure  2-1.  The  values 
are  quite  large  and  it  seems  probable  that  Stark  broadening  of  the  lines 
will  cause  a  greater  effective  depression  of  the  levels.  The  Stark 
broadening  has  been  discussed  by  GrJem,Kolb  and  Shen^  and  it  may 

be  necessary  to  carry  out  detailed  computations  based  on  their  theory. 
However,  the  fact  that  the  highly-excited  levels  are  in  essentially 
thermal  equilibrium  and  play  only  a  small  role  in  determining  the  effect¬ 
ive  recombination  suggests  chat  the  solution  of  the  equations  is  insensitive 
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Figure  2*1.  Highest  Bound  State  of  an  Electron  i.n  Hydrogen 
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to  the  nuniber  of  levels  Included.  The  Stark  broadening  will  be  Inqpor- 
tant,  however.  In  determining  the  opacity,  which  In  turn  effects  the 
recombination. 

2.5  Computation  of  RecoBd>lnatlon  Rate  In  Optically  Thin  Hydrogen 

(2-7) 

Adopting  the  method  of  Bates  and  Kingston  ,  we  first  change 
the  variables  In  Equation  (2-1)  from  n(p)  to  d(p)  by  Introducing  the 
relationship 

n(p)  -  n^  n(X‘^)  0)^  (h^/2nmkTg)^^^  exp(Ip/kT^)  e(p)  (2-9) 

(where  m  is  bmss  of  electron,  k  Is  Boltzmann's  constant,  h  Planck's 
constant  and  1^  the  Ionization  potential  of  this  p^^  level.) 

The  new  variables  d(p)  are  the  ratios  of  populations  In  the  p^^  levels 
to  the  populations  derived  from  the  Saha  relationship.  At  equilibrium 
d(p)  Is  unity. 


The  Introduction  of  (2-9)  Into  Equation  (2-1)  gives  an  equation  of 
the  form* 


b 

P 


a 

PP 


0(p) 


p=  1,  2,  3,  ... 


(2-10) 


Assuming  with  Bates  and  Kingston  that 


_  0 

dt  " 


p  =  2,  3,  4, . . . 


(2-11) 


*the  a's  do  not  have  same  meanings  In  Equation  (2-10)  as  In  Equation 
(2-1) 
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we  obtain 


®jp  •‘p 


a  0(p) 

PP 

p  =  2|  3|  4| . . . 


(2-12) 


with  the  equation  for  the  ground  state  where  p  »  1  treated  as  a  special 
case: 


^  0(1)  =  y  ajj  0(J)  +  -  a^j  0(1)  (2-13) 

jVp 

If  the  series  of  equations  in  (2-12)  are  truncated  at  some  level  k,  the 
Equations  (2-12)  are  in  the  form  of  k-1  linear  equations  with  k  unknowns, 
the  unknowns  being  taken  as  the  0's.  By  inversion  of  the  matrix,  a 
solution  may  be  obtained  for  all  0's  in  terms  of  one  unknown,  0(1).  Be¬ 
cause  of  the  linearity  of  the  equations,  the  solution  will  have  the  form 

0(p)  *  *0  + 

p  =  2,  3...  k  (2-14) 

The  right  side  of  Equation  (2-13)  can  then  be  evaluated  also  in  terms  of 
the  single  unknown  0(1)  and  taking  account  of  Equations  (2-9)  and  (2-7), 
one  suy  obtain  the  effective  recombination  coefficient  in  the  form: 

•  ‘l  <2-‘5) 

Values  of  a^  and  a^  in  Equation  (2-15)  have  been  obtained  for  an 
optically  thin  hydrogen  plassa  by  machine  methods  for  electron  tempera¬ 
tures  ranging  from  2S0'’K  to  64,000''K  and  for  all  electron  densities. 
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indices  give  the  power  of  10  by  which  the  entries  must  be  multiplied. 
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Indices  give  the  power  of  10  by  which  the  entries  must  be  multiplied. 


levels  up  to  n  =  10  being  Included  in  the  confutation.  Preliminary 
values  of  a^  and  aj^  haw  been  provided  by  Bates  and  are  given  In  Tables 
2-la  and  2-lb,  respectively. 

There  are  certain  Irregularities  In  the  values  as  given  In  the 
tables  which  suggest  there  might  be  an  error  In  parts  of  the  tabulation. 
Until  a  further  check  of  the  computations  has  been  coofleted  the  values 
In  the  table  should  be  regarded  as  uncertain  by  a  factor  of  2. 

2.6  Discussion  of  Computation  and  Results 

In  order  to  Illustrate  the  significance  of  the  differential 
Equations  (2-10),  the  algebraic  Equations  (2-12),  and  the  solution  as 
given  In  Tables  2-la  and  2-lb,  It  will  be  worthwhile  to  consider  the 
changes  which  take  place  In  the  populations  of  the  various  excited  levels 
of  a  hydrogen  plasma.  Suppose  that  for  an  optically  thin  hydrogen  plasma, 
the  temperature  and  number  density  of  the  free  electrons  are  maintained 
constant  by  some  arbitrary  means.  Under  these  conditions,  the  rate  of 
change  of  the  populations  of  all  the  bound  states  are  given  by  the 
differential  Equations  (2-1),  and  for  given  Initial  conditions  which 
assign  a  value  to  the  population  of  each  of  the  bound  states  at  the  out¬ 
set,  a  solution  of  the  Equations  (2-1)  will  give  the  populations  of  each 
of  the  states  for  all  subsequent  time.  The  general  character  of  this 
solution  Is  Illustrated  schematically  in  Figure  2-2. 

In  the  figure,  the  populations  of  the  several  states  are  described  by 
the  variable  0,  which  Is  the  ratio  of  the  number  In  a  state  to  the 
number  In  that  state  under  conditions  of  thermodynamic  equilibrium. 

Since  we  are  supposing  that  the  number  density  and  temperature 
of  the  free  electrons  remains  constant. 
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6(p)  will  be  proportional  to  n(p) .  In  Figure  2-2,  It  Is  supposed  that 
the  ground  state  Is  unpopulated  at  the  outset,  but  the  Initial  popula¬ 
tions  of  the  various  excited  states  are  taken  at  values  which  have  no 
particular  significance.  As  indicated  In  the  figure,  at  the  outset  there 
will  be  a  very  rapid  change  in  the  populations  of  the  upper  states.  In 
particular,  the  value  of  6  for  the  upper  states  will  change  nuch  store 
rapidly  than  that  of  the  ground  state.  After  a  short  time,  the  rapid 
changes  In  6  of  the  upper  states  will  decline,  and  a  pseudo-steady  state 
will  occur  such  that.  If  viewed  on  a  short  time  scale,  all  the  populations 
remain  constant.  There  will,  however,  be  a  non-zero  rate  of  Increase  In 
the  ground  state,  which  will  make  Itself  felt  over  a  longer  period  of 
tlsie.  This  pseudo-steady  state  condition  Is  representated  In  Figure  2-2 
as  occurring  Just  before  the  first  tlsie  break,  and  It  is  this  condition 
which  is  described  by  the  algebraic  Equations,  2-12  and  by  their  solu¬ 
tion  obtained  taking  the  recombination  coefficient  or  ■  a^.  At  an  appre¬ 
ciably  later  period  of  tlsie,  as  represented  after  the  first  time  break 
In  Figure  2-2,  the  Increase  In  6  for  the  ground  state  will  be  appreciable, 
and  there  will  be  a  corresponding  change  In  the  pseudo-steady  state  values 
of  the  populations  of  the  various  higher  states.  This  situation  Is  re¬ 
presented  In  the  middle  time-period  of  Figure  2-2,  and  is  also  described 
by  the  solution  to  the  algebraic  Equations,  2-12,  but  now  since  0(1)  Is 
not  equal  to  zero.  It  Is  necessary  to  take  the  effective  recombination 
rate,  as  o  -  a^  -  a^6(l). 

The  populations  of  the  various  states  continue  to  change  as  recomb¬ 
ination  proceeds,  and  the  right  hand  section  of  Figure  2-2  indicates 
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the  general  nature  of  these  changes  on  a  much  different  time  scale  than 
the  portions  of  the  figures  at  the  left  and  the  center.  The  populations 
will  eventually  approach  asymptotic  values  at  which  they  are  constant , 
and  at  which  time  the  net  recombination  ceases.  This  condition,  which 
Is  indicated  at  the  extreme  right  of  the  figure  represents  the  steady 
state  under  conditions  of  radiation  loss.  It  can  be  seen  that  for  this 
condition,  the  population  of  the  ground  state  exceeds  that  which  It  would 
have  for  therinpdynamlc  equilibrium.  The  population  of  the  ground  state 
can  be  obtained  readily  from  the  Tables  2-la  and  2-lb  by  setting  a  * 
a  -  a,  0(1)  -  0. 

In  simplifying  the  differential  Equations  2-10,  to  the  corresponding 
algebraic  Equations,.  2-12,  one  gives  up  any  possibility  of  obtaining 
Information  about  the  Initial  variaMon  shown  at  the  extreme  left  of 
Figure  2-2.  Except  for  this  initial  response,  which  Is  unimportant  for 
the  applications  considered  here,  the  effective  recombination  rate 
and  the  populations  of  all  the  states  can  be  obtained  throughout  the 
period  Illustrated  In  the  second  and  third  time  blocks  In  the  Figure 
from  the  algebraic  equations. 

It  should  be  observed  that  In  the  proceedure  followed  here,  one 
supposes  that  the  temperature  and  number  density  of  the  free  electrons 
and  the  conditions  on  light  absorption  as  described  by  the  X'P^i'unBters 
are  given  conditions  of  the  problem.  To  do  otherwise  would  make  the 
equations  (either  the  differential  Equations  2-10  or  the  algebraic 
Equations  2-12)  non-linear  and  would  vastly  Increase  the  difficulties  of 
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obtaining  a  solution.  The  solution  obtained  here  can  nonetheless  be 
applied  to  important  problems  in  which  the  temperature,  density  and 
light  absorption  in  the  plasma  change  with  time  because  the  pseudo¬ 
steady  state  to  which  Equation  (2-12)  applies  is  achieved  in  times 
short  compared  to  those  for  significant  variation  in  the  plasma  properties. 
This  matter  will  be  discussed  further  in  Section  3  of  Volume  I. 

2.7  Plasma  Other  Than  Hvdroaen 

The  general  plasma  has  several  ionized  components  each  of  which 
can  exist  in  several  states  of  ionization.  Nevertheless,  the  general 
plasma  may  be  treated  in  a  manner  similar  to  that  employed  for  the  hy¬ 
drogen  plassia,  but  with  certain  modifications. 

Just  as  for  the  hydrogen  plasma,  it  is  necessary  to  collect  to¬ 
gether  all  the  radiative  and  collislonal  transition  probabilities.  Un¬ 
fortunately,  no  systematic  tabulation  of  these  probabilities  has  been 
given  even  for  ground  states  of  systems  and  the  collection  of  the  data 
is  a  major  task.  We  have  made  a  beginning  in  assembling  the  data  and 
the  tabulated  value  are  included  in  Appendix  A. 

The  basic  equations  are  similar  to  (2-1)  and  their  solution  pre¬ 
sents  no  problem  provided  heavy  particle  collisions  can  be  ignored. 

With  this  proviso,  there  is  a  set  of  equations  of  the  form  (2-1)  for 
each  component  of  the  plasma  and  these  sets  of  equations  are  uncoupled. 

The  quasi -steady  state  method  may  then  be  applied  to  each  component 
separately. 
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It  does,  in  fact,  appear  reasonable  to  neglect  heavy  particle 
collisions  except  for  charge  transfer  processes  involving  the  ground 
states.  This  exception  arises  because,  although  charge  transfer  into 
the  ground  state  of  a  system  may  be  inefficient,  it  does  couple  compo¬ 
nents  together  and  it  has  a  permanent  recombining  effect.  It  does  not 
introduce  any  difficulties  into  solving  the  equations  since  the  ground 
state  equation  receives  special  consideration  in  any  event. 


18 


REFERENCES 


2-1  Ivanov-Kholodnyl,  G.  S.,  et  al.  "Ionization  and  Excitation  of 
Hydrogen  1.  Elementary  Processes  for  the  Upper  Levels",  Soviet 
Astronomy-AJ  754  (1961). 

2-2  Baker,  J.  G.  and  Menzel,  D.  H.,  Astrophys,  J.  No.  1  (July, 
1937). 

2-3  Burgess,  A.,  "The  Hydrogen  Recombination  Spectrum",  Mon.  Nat. 
Roy.  Astron.  Soc.  118.  478  (1958). 

2-4  Giovanelli,  R.  G.,  "Hydrogen  Atmospheres  in  the  Absence  of 
Thermodynamic  Equilibrium" 

I.  "The  Population  of  the  Upper  Atomic  States,  Australian 
J.  Sci.  Research  Al,  275  (1948) 

2-5  Giovanelli,  R.  G.,  "The  Population  of  the  Lower  Atomic  States" 
Australia  J.  Sci.  Research  M,  289  (1948). 

2-6  D'Angelo,  N.,  "Recombination  of  Ions  and  Electrons",  Phys .  Rev. 
121.  505  (1961). 

2-7  Bates,  D.  R.  and  Kingston.  A.  E.,  "Recombination  through 
Electron-Electron  Collisions",  Nature  189.  652  (1961). 

2-8  Griem,  H.  R.,  Kolb,  A.  C.,  and  Shen,  K.  Y.,  "Stark  Broadening 
of  Hydrogen  Lines  in  a  Plasma",  Phys.  Rev.  116.  4  (1959). 


APPENDIX  A 

A.  Dalgarno,  A.  Naqvt,  E.  Ashley  and  H.  E.  Stubbs 
ATOMIC  PROCESSES  IN  A  RECOMBINING  PLASMA 

The  time  dependence  of  the  population  of  each  state  of  a  recombining 
plasma  can  be  obtained  from  a  consideration  of  the  several  processes  by 
which  the  atoms  enter  and  leave  a  given  state.  Both  radiative  and 
colllsional  processes  must  be  considered. 

A. 1  Radiative  Transitions  Between  Bound  States 

The  rate  of  depopulation  of  the  state  p  by  spontaneous  transitions 
to  lower  states  q,  may  be  written 

n(p)^  A(p,q)  (A-1) 

q<p 

where  A(p,q)  is  the  probability  of  a  spontaneous  transition  from  the 
state  p  to  the  state  q  measured  in  sec 

If  the  plasma  is  not  optically  thin,  stimulated  emissions  and  absorp¬ 
tions  occur.  For  the  rate  of  upward  transitions  which  depopulate  the 
state  p,  we  write 

Y  “r  A(r,p)  (A-2) 

Gj  / _  r 

. 

where  is  the  statistical  weight  of  the  state  p  and  X(p,r)  is  a  factor 
which  takes  account  of  the  opacity.  In  particular  for  a  plasma  in  ther¬ 
modynamic  equilibrium  with  the  radiation  field 


iO 


X  (r 


,p)  -  |exp  (hv^p/kT)j 


-1 


(A-3) 


where  energy  of  the  r  -  p  transition.  For  an  optically  thin 

plasma, X vanishes.  The  rate  of  stimulated  downward  transitions  can  be 
written 


n(p)  A  X (p,q)  A(p,q) 


(A-4) 


The  p  state  can  also  be  populated  by  bound-bound  radiative  transi¬ 
tions  and  we  have  for  the  rates  of  population  by  the  three  processes 


^  n(r)  A(r,p) 
r>p 

^  n(r)  X(r,p)  A(r,p) 
r>p 

“p  Z!  X  (q,p)  A(p,q). 


(A-5) 


(A-6) 


(A-7) 


<j<p  q 

The  probability  of  a  spontaneous  transition  from  a  level  p  to  a 
lower  level  q  may  be  expressed  in  the  form 


A(p,q)  -  — y  S(p,q) 
3hc  u 


(A-8) 


where  v  is  the  frequency  of  the  emitted  radiation,  is  the  statistical 
weight  of  the  p^^  level  and  S(p,q)  is  the  transition  strength,  defined 


according  to 


S(p,q)  -  |<  J  ^j'Fq  >  | 


(A-9) 
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where ^  and  ^  are  the  configurations  of  the  Initial  and  final  states 
P  q 

and  r  is  the  position  vector  of  the  J  electron.  In  the  case  of  Russell- 
Saunders  coupling,  S(p,q)  reduces  to  the  form 


S(p.q)  -/(L)  /(M) 


(A- 10) 


where  j  (L)  is  the  relative  strength  of  a  multiplet  within  a  transition 

array,  /(M)  is  the  relative  strength  of  a  line  within  a  multiplet  and 

0  is  a  transition  Integral  Involving  radial  wave  functions  only.  Tables 

of  :/iL)  and  (M)  are  available  for  most  multiplets  (Goldberg^^ 

White  and  Eliason^^  Russell^^  ,  Allen^^  and  if  necessary,  they 

fA-5) 

can  easily  be  extended  following  Rohrlich^  ').  The  transition  integral 
0  presents  greater  difficulties.  It  R(n.i|r)  is  the  initial  radial  wave 
function  of  the  active  electron  and  R(n'i'|r)  is  the  final  wave  function. 


‘“I 


CO  I 

R(ni:  |r)r\(n'i:'  |r)  dr 


(A-11) 


where  is  the  lower  of  the  azimuthal  quantum  numbers and  £'  which 
characterize  the  active  electrons. 


2 

For  hydrogen,  values  of  a  have  been  computed  by  Green, 

Rush  and  Chandler^^  and  the  corresponding  values  of  A(nil,  for 

n  and  n*  less  than  or  equal  to  8  are  presented  in  Table  A-1.  Of  interest 
also  are  the  total  transition  probabilities  A(n,n')  obtained  by  summing 
A(ni;,  n'i')  over  all  possible  values  of  2  and  £'  consistent  with  n  and 
n',  assuming  that  the  sub-levels  are  distributed  according  to  their 
statistical  weights.  The  total  transition  probabilities  so  obtained  are 


I 

I 

I 

I 


I 


I 


I 

3 
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24 
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presented  in  Table  A-2.  There  are  ocaslonal  differences  from  previous 

(A- 7) 

tabulations  of  Menzel  and  Pekerls 

For  more  complex  systems,  approximations  must  be  made.  For  many 

cases  and  especially  for  the  excited  states  of  light  elements,  the 

f  A*8) 

Coulond)  approximation  of  Bates  and  Damgaard'^  '  should  be  adequate.  In 

this  approximation  the  active  electron  Is  regarded  as  moving  In  the 

field  of  an  effective  nuclear  charge  being  determined  from  the  observed 

Ionization  potentials  for  the  states.  Bates  and  Damgaard  have  presented 
2 

tables  from  which  a  may  be  determined  for  transitions  Involving  small 
values  of  ^ and  /  and  we  have  used  them  to  calculate  the  transition 
probabilities  of  the  more  important  lines  of  Li  I,  Li  II,  C  I,  C  II, 

C  III,  N  I,  N  II,  N  III,  0  I,  0  II  and  0  III.  The  results,  illustrative 
of  the  Information  which  may  be  obtained  by  the  Bates  and  Damgaard 
methods,  are  presented  in  Table  A*3.  Also,  Table  A-4  lists  the  results 
of  calculations  by  this  method  of  the  transition  probabilities  of 
resonance  lines  of  a  number  of  atoms  of  interest. 

For  heavier  elements  It  may  not  be  possible  to  apply  the  Coulomb 
approximation  since  the  Ionization  potentials  are  not  always  known  and 
a  less  empirical  approximation  scheme  Is  required  A  theory  using 
screening  constants,  which  is  described  in  the  following  appendices, 
is  capable  of  dealing  with  these  cases. 
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TABLE  A-2 

TOTAL  T1UH8ITIOH  TOOBABIUTIES  TOR  HTDROGEH  IN 


1 


2 

4.696  (8) 

- 

- 

- 

- 

3 

5.572  (7) 

4.410  (7) 

- 

- 

- 

4 

1.278  (7) 

8.419  (6) 

8.986  (6) 

- 

- 

S 

4.122  (6) 

2.530  (6) 

2.201  (6) 

2.699  (6) 

- 

6 

1.643  (6) 

9.732  (5) 

7.783  (5) 

7.710  (5) 

1.025 

(6) 

7 

7.564  (5) 

4.389  (5) 

3.358  (5) 

3.041  (5) 

3.251 

(5) 

8 

3.867  (5) 

2.215  (5) 

1.650  (5) 

1.423  (5) 

1.381 

(5) 

9 

2.141  (5) 

1.216  (5) 

8.904  (4) 

7.455  (4) 

6.904 

(4) 

10 

1.262  (5) 

7.122  (4) 

5.156  (4) 

4.232  (4) 

3.798 

(4) 

11 

7.829  (4) 

4.397  (4) 

3.156  (4) 

2.555  <4) 

2.245 

(4) 

12 

5.063  (4) 

2.834  (4) 

2.021  (4) 

1.620  (4) 

1.402 

(4) 

13 

3.391  <4) 

1.893  (4) 

1.343  (4) 

1.068  (4) 

9.143 

(3) 

14 

2.340  (4) 

1.303  (4) 

9.211  (3) 

7.284  (3) 

6.181 

(3) 

IS 

1.656  (4) 

9.210  (3) 

6.490  (3) 

5.108  (3) 

4.306 

(3) 

16 

1.199  (4) 

6.658  (3) 

4.680  (3) 

3.669  (3) 

3.007 

(3) 

17 

8.852  (3) 

4.910  (3) 

3.444  (3) 

2.692  (3) 

2.248  (3) 

18 

6.650  (3) 

3.685  (3) 

2.580  (3) 

2.012  (3) 

1.674 

(3) 

19 

5.074  (3) 

2.809  (3) 

1.964  <3) 

1.528  (3) 

1.268 

(3) 

20 

3.925  (3) 

2.172  (3) 

1.517  (3) 

1.178  (3) 

9.746 

(2) 

A(n,n') 

6  7  8  9  10  '  11 


4.558 

(5) 

. 

. 

. 

1.560 

(5) 

2.271 

(5) 

- 

- 

- 

7.061 

(4) 

8.232 

(4) 

1.232 

(5) 

- 

- 

3.686 

(4) 

3.903 

v4) 

4.674 

(4) 

7.147 

(4) 

- 

2.108 

(4) 

2.116 

(4) 

2.299 

(4) 

2.811 

(4) 

4.374 

(4) 

1.288 

(4) 

1.250 

(4) 

1.286 

(4) 

1.426 

(4) 

1.773 

(4) 

2.797 

( 

8.267 

(3) 

7.841 

(3) 

7.799 

(3) 

8.187 

(3) 

9.226 

(3) 

1.163 

( 

5.523 

(3) 

5.153 

(3) 

5.007 

(3) 

5.007 

(3) 

5.414 

(3) 

6.182 

( 

3.813 

(3) 

3.514 

(3) 

3.357 

(3) 

3.323 

(3) 

3.422 

(3) 

3.876 

( 

2.705 

(3) 

2.469 

(3) 

2.329 

(3) 

2.267 

(3) 

2.278 

(3) 

2.376 

( 

1.965(3) 

1.780 

(3) 

1.663 

(3) 

1.597 

(3) 

1.577 

(3) 

1.605 

( 

1.456 

(3) 

1.311 

(3) 

1.215 

(3) 

1.156 

(3) 

1.126 

(3) 

1.126 

< 

1.099 

(3) 

9.844 

(2) 

9.065 

(2) 

8.550 

(2) 

8.248 

(2) 

8.136 

( 

8.419 

(2) 

7.512 

(2) 

6.882 

(2) 

6.449 

(2) 

6.170 

(2) 

6.023 

< 

The  bT«ck«ti  giva  th«  povtrs  of  10  by  which  th«  ancrio  me  ba  Bulcipllad. 


DROGEN  IN  lac 


-1 


8 


9  10  ’  11 


12 


13  14  IS  16  17  18  19 


0 

1.232 

(5) 

. 

- 

- 

- 

►) 

4.674 

(4) 

7.147 

(4) 

- 

- 

- 

0 

2.299 

(4) 

2.811 

(4) 

4.374 

(4) 

- 

• 

0 

1.286 

(4) 

1.426 

(4) 

1.773 

(4) 

2.797 

(4) 

- 

1) 

7.799 

(3) 

8.187 

(3) 

9.226 

(3) 

1.163 

(4) 

1.856 

(4) 

1) 

5.007 

(3) 

5.007 

(3) 

5.414 

(3) 

6.182 

(3) 

7.880 

(3) 

1.270 

(4) 

1) 

3.357 

(3) 

3.323 

(3) 

3.422 

(3) 

3.876 

(3) 

4.268 

(3) 

5.493 

(3) 

8.928 

(3) 

) 

2.329 

(3) 

2.267 

(3) 

2.278 

(3) 

2.376 

(3) 

2.595 

(3) 

3.024 

(3) 

3.923 

(3) 

6.425 

(3) 

•  * 

) 

1.663 

(3) 

1.597 

(3) 

1.577 

(3) 

1.605 

(3) 

1.692 

(3) 

1.865 

(3) 

2.191 

(3) 

2.863 

(3) 

4. 

717 

(3) 

) 

1.215 

(3) 

1.156 

(3) 

1.126 

(3) 

1.126 

(3) 

1.158 

(3) 

1.232 

(3) 

1.368 

(3) 

1.619 

(3) 

2. 

128 

(3) 

3.528 

(3) 

) 

9.065 

(2) 

8.550 

(2) 

8.248 

(2) 

8.136 

(2) 

8.219 

(2) 

8.327 

(2) 

9.139 

(2) 

1.022 

(3) 

1. 

217 

(3) 

1.609 

(3) 

2.679 

(3) 

) 

6.882 

(2) 

6.449 

(2) 

6.170 

(2) 

6.023 

(2) 

6.002 

(2) 

6.116 

(2) 

6.394 

(2) 

6.899 

(2) 

7. 

763 

(2) 

9.291 

(2) 

1.234 

(3) 
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TABLE  A-3 


BOUND-BOUND  TRANSITION  PROBABILITIES  A(p,q)  FOR  THE  IMPORTANT 
LINES  OF  Lil,  LlII,  Cl,  ClI,  CIII,  NI,  Nil,  NIII,  01,  011  AND  OIII 

Lil  Lines 


Transition 

^2 

•^1 

2 

O’ 

A(10^sec 

2p^P°-2s^S 

1/2 

1/2 

6709.8 

5.424 

3.63 

3/2 

1/2 

6709.6 

5.424 

3.63 

3p^P°-2s^S 

1/2 

1/2 

3233.6 

0.01641 

0.098* 

3/2 

1/2 

3233.6 

0.01641 

0.098* 

4p^r°-2s^S 

1/2 

1/2 

2742.0 

0.01392 

0.136* 

3/2 

1/2 

2742.0 

0.01392 

0.136* 

5p^P°-2s^S 

1/2 

1/2 

2563.0 

3.90x10’^ 

0.0468* 

3/2 

3/2 

2563.0 

3.90x10"^ 

0.0468* 

Lill  Lines 


Transition 

X(X) 

B 

- _ 

2 

a 

AdO^sec 

ls2p^P°-ls^  ^S 

'.v9.3 

6 

3.72.:10'^ 

l.OO* 

ls3p^P®-ls^  ^S 

178.0 

6 

9.647x10*^ 

691 

ls4p^P°-ls^  ^S 

171.6 

6 

3.219x10’^ 

257 

Cl  Lines 

2p(^P°)38^p"-2p^  ^P 

1657 

18 

0.01397 

1.24* 

2s2p^  ^D°-2s^2p^  ^P 

1561 

30 

00.9563 

101 

2s2p^  ^P°-2s^2p^  ^P 

1329 

18 

0.01284 

2.21* 

2p(^P°)48^P°2p^  ■*P 

1280 

18 

0.001488 

0.287* 

^below  907.  maximum  value 
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TABLE  A-3  (continued) 


CII  Ltnea 


'i’raiisl  1  ion 

X(X) 

>’<M) 

or^ 

A(10^sec  S 

282p^  ^D-28^2p^P° 

1335 

1.096 

2s2p^  ^S-2s^2p^P® 

1037 

1.047 

282p^  ^P-28^2p^P° 

904 

0.8418 

28^3s^S-28^2p^P° 

858 

6 

0.062 

59.4* 

28^(^S)3d^D-28^2p^P° 

687 

0.04008 

1500 

28^(^S)48^S-28^2p^P® 

636 

6 

0.00817 

20.8* 

28^(^S)4d^D28^2p^P° 

595 

60 

0.01357 

78.0* 

28^(^S)58^S28^2p^P° 

577 

6 

0.002743 

8.64* 

28^(^S)Sd^D-28^2p^P° 

560 

60 

0.005464 

37.7* 

cm  Lines 


28(^S)2p^P°-28^ 

977 

6 

0.5723 

248 

28(^S)3p^P®-28^ 

386 

6 

0.05579 

392 

28(^S)4p^P°-28^ 

310 

6 

0.006529 

88.5 

NI  Lines 


2s^2p^(^P)38S- 
28^2p^  ^S° 

1200 

12 

0.007382 

0.860* 

2p^(^P)48S-2p^  ^S° 

964 

12 

0.0007751 

0.175* 

2p^(^P)3dS-2p^4S° 

953 

120 

0.04008 

93.7 

2s2p^  S-2s^2p^  ^S® 

1135 

1135 

0.5631 

2p^(^P)58S-2p^  V 

911 

12 

0.0003190 

0.0853* 

*Below  90%  raaxlnun  value 
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TABLE  A-3  (continued) 
Nil  Lines 


Transition 

x(X) 

- 2 - 

a 

OOgjOj 

2s2p'’  ^D°-28‘^sp‘^  •’P 

1085 

30 

0.3732 

148 

2s2p^  ^P°-2s^2p^  ^P 

916 

18 

0.3884 

195 

2p(^P°)3s^P°-2p^  ^P 

671 

18 

0,01671 

22.3 

2s2p^  ^S°-2s^2p^  ^P 

645 

12 

0,01833 

55.2 

2p(^P°)3d^D°-2p^  ^P 

534 

135 

0.00729 

109 

NIII  Lines  (C.A.) 


2s2p^  ^D-2s^2p^P° 

991 

0.1947 

2s2p^  ^S-2s^2p^P° 

764 

0,2135 

2s2p^  ^P-2s^2p^P° 

686 

0,1986 

2s^(^S)3s^S-2s^2p^P° 

452 

6 

0.01173 

77.1 

2s^(^S)3d^D-2s^2p^P° 

374 

60 

0.01028 

238 

2s^(^S74s^S-2s^2p^P° 

332 

6 

0  001995 

33.0 

2s2p(^P°)3p^P- 

2s^2p^P° 

324 

0.006543 

2s^(^S)4d^D-2s^2p^P° 

309 

60 

0.00149 

61.2 

01  Lines 


2p^(V)3s^S°-2p^  ^P 

1303 

12 

0.01261 

4.62 

2p^(^S°)48^S°-2p^  ^P 

1040 

12 

0.001584 

1.14 

2p^(^S°)3dV-2p^  ^P 

1027 

0.03764 

2p^(V)3sV-2p^  ^P 

989 

15 

0.0007629 

0.160* 

2p^(^S°)28^S°-2p^  ^P 

977 

12 

0.004858 

4.20 

*Below  90%  maximum  value 
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TABLE  A-3  (continued) 


on  Lines  (C.A.) 


Transition 

1 

uX) 

afw 

2 

ff 

AClO^sec  S 

2s2p^  S-2s^2p^ 

834 

0.2670 

2s^2p^(^P)3aS- 

2s22p^  V 

539 

12 

0.006351 

8.2* 

2p^(^P)3dS-2p^  ^S° 

430 

120 

0.005050 

128  j 

2p^(^P)3dS-2p^  ^S° 

430 

0,005356 

1 

01 II  Lines 


2s2p^  ^D°-2s^2p^  ^P 

834 

30 

0.1417 

98,4* 

282p^  ^P°-2s^2p^  ^P 

703 

18 

0  1538 

179 

282p^  ^S°-2s^2p^  ^P 

508 

12 

0,1643 

1014 

2s^2p(^P°)3s^P°- 

2s^2p^  ^P 

374 

18 

0  006293 

48.6* 

2s^2p(^P°)3d^D°- 

2s^2p^  ^P 

306 

135 

0,005453 

346 

*Bclow  90%  maximum  value 
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TABLE  A-4 

TRANSITION  PROBABILITIES  FOR  RESONANCE  LINES 


Element 

Transition 
(Upper  -  Lower) 

^2-^1 

X(X) 

2 

O’ 

AdO^sec 

Rb 

I 

4p^(^S)5p  ^P®- 
4p^(^S)58  ^S 

1/2  -  1/2 

7947.60 

8.82 

35.5 

Sr 

I 

58(^S)5p  ^P°- 
5s^  ^S 

1  -  0 

4607.331 

4.72 

195 

Sr 

II 

4p^(^S)5p  V-  1 
4p^(^S)58  ^S 

1/2  -  1/2 

4077.714 

4.92 

147 

Zr 

I 

4d^58(aS)5p  ^F°- 

4  -  4 

6121.49 

3.19 

52.5 

2  2  3 

4d  58  ^F . 

2  -  2 

6134.58 

49.6 

Zr 

II 

4d^(a^F)5p  V-  5 
4d^(a^F)58  S 

1/2  -  4  1/2 

3391.96 

3.59 

186 

Nb 

I 

4d^58(a^F)5p  ^F°- 
4d^(a^D)58 

5  1/2-4  1/2 

4058.933 

3.39 

410 

Mo 

I 

4d^(a^S)5p  ^P°- 

4  -  3 

3798  359 

3.35 

123 

4d^(a^S)58  ^S 

3-3 

3864.115 

117 

2-3 

3902.968 

114 

Tc 

I 

5d^58(a^S)5p  ^P°- 

3  1/2-2  1/2 

4298.2 

3.15 

160 

5d^58^  ®S 

2  1/2-2  1/2 

4263.4 

164 

1  1/2-2  1/2 

4239.4 

167 

Ru 

I 

4d^(aS)5p  V- 

7  .  A  .  S 

Entire 

Multlplet 

3799.124 

3.27 

120 

4d'(a^F)58  P 
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TABLE  A-4  (continued) 


Element 

Transition 
(Upper  -  Lower) 

^2  '  ‘•l 

uX) 

2 

CT 

Rh  1 

4d®(^F)5p  V- 

3  1/2-4  1/2 

3692,357 

3,21 

4d®(^F)58 

2  1/2-3  1/2 

3657  987 

Ag  1 

4d^°(^S)5p  ^P°- 

1  1/2-  1/2 

3280  682 

2.95 

4d^°(^S)58 

1/2  -  1/2 

3382.890 

Cd  I 

4d^°580^S)5p 
4d'°582  4 

1  -  0 

2288  7 

1  99 

In  I 

58^(^S)68  ^S- 

1/2-  1  1/2 

4511  310 

0  890 

58^(^S)5p 

1/2  -  1/2 

4101  764 

Sn  1 

58^5p(^P°l)6s  ^P* 
58^5p^  ^P^ 

*-1-2 

3175  046 

0  202 

0  -  1 

3034  120 

Cs  I 

5p^(^S)6p  ^P°- 

1  1/2-  1/2 

8521  10 

9  94 

5p^(.^S)68 

1/2  -  1/2 

8943  50 

Ba  1 

68(^S)6p  ^P®- 

68^  's 

1  -  0 

5535  484 

5  78 

Ba  11 

5p^(^S)6p  ^P°- 

1  1/2-  1/2 

4554  033 

5  96 

5p^(4)68 

1/2  -  1/2 

4934.086 

U  I 

5d68(a^O)6p  ^D®- 

2  1/2-2  1/2 

7270.11 

5.10 

2  2 

5d68  D 


AOO^sec 

115 

1082 

169 

154 

671 

39  2 
26,1 
10  6 
29.2 
32.5 
28  1 
138 

127 

100 

50.0 
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TkSLE  A-A  (continued) 


Element 

Trnniltlon 
(Upper  -  Lower) 

J  -  J 

2  1 

x(X) 

2 

0 

A(10^eec 

Ih  1 

4d®(^P)5p  V- 

3  1/2-4  1/2 

3692. 3S7 

3.21 

115 

4d®(^P)5i  S 

2  1/2-3  1/2 

3657.987 

1082 

Ag  I 

4d^°(^S)5p  ^P®- 

1  1/2-  1/2 

3280  682 

2.95 

169 

4d^°(^S)5e 

1/2  •  1/2 

3382.890 

154 

Cd  I 

4d^°5eO^S)5p  ^P® 

-  1-0 

2288.7 

1.99 

671 

4d'°5.2  4 

In  I 

5e^(4)6e  V 

1/2-  1  1/2 

4511.310 

0.890 

39  2 

5e^(^S)5p  V 

1/2  -  1/2 

4101.764 

26.1 

Sn  I 

5e^5p(^pji)6e 

*-1-2 

3175.046 

0  202 

10.6 

5e^5p^  ^P^ 

0-1 

3034.120 

29.2 

Cs  I 

5p®(^S)6p  ^P®- 

1  1/2-  1/2 

8521 . 10 

9.94 

32.5 

5p^(^S)6o 

1/2  -  1/2 

8943.50 

28.1 

Bn  I 

6o(^S)6p  ^P®- 

1-0 

5535.484 

5.78 

138 

6e^ 

Bn  11 

5p^(^S)6p  ^P^- 

1  1/2-  1/2 

4554.033 

5.96 

127 

5p®(^S)6e 

1/2  -  1/2 

4934.086 

100 

U  1 

5d6e(e^D)6p  ^D®- 

2  1/2-2  1/2 

7270.11 

5.10 

50.0 

5<I6e* 

A.  2  lUdl«tlv«  Boun<t-rr««  Tr>n«lttona 


In  order  to  calculate  the  coefficient  a(p)  for  radiative  capture  of 
electrons  into  state  p 

+  e  -»  X(p)  +  hv 

we  must  make  some  assumption  about  the  electron  velocity  distribution. 
Since  the  electrons  exchange  energy  readily  in  elastic  collisions,  they 
may  be  described  by  a  Maxwellian  distribution  characterised  by  an  elec¬ 
tron  temperature  The  recombination  coefficient  is  accordingly 

a<pii.)  -  \  .xp(ip/ki,) 

*  J  Q^(P)  (hv)^  exp(-hv/kT^)  [l  +X(v)]  d(hv) 

I 

^  (A- 12) 

where  is  the  ionization  potential  of  the  p^^  level,  hv  is  energy  of 
the  radiation  photon  and  Q^(p)  is  the  photoionization  cross  section 
corresponding  to  the  process 

X(p)  +  hv  -  +  e  (A- 13) 

The  rate  of  population  by  radiative  capture  is 

a(p|T^)  n(X*)  (A-14) 

^ere  n  and  n(X^)  are  respectively  the  electron  and  X^  ion  number  den- 
sities. 
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The  rate  of  depopulation  by  photolonlzatlon  may  be  written 


n(p)  P(P) 


where 


8Jt 

P(P)  -  -J-J —  / 

he  ^ 

P 


(A-15) 


The  quantal  description  of  bound-free  transitions  is  entirely  similar 
to  that  of  bound-bound  transitions  except  that  the  radial  wave  function 
R(n'^'|r)  in  (A- 11)  becomes  R(e'C'|r),  describing  an  electron  of  positive 
energy  e*  moving  in  the  field  of  the  ion.  Thus,  the  cross  section  for 
the  transition  of  an  electron  from  the  bound  state  (ni!)  to  the  continuum 
state  of  energy  c '  is 

Ancea 

Q^(ni!;€’)  - 


'I' 


|g(n^;e'JJ')  | 


(A-16) 


where  I  is  the  ionization  potential,  C.,  are  simple  algebraic  quantities 
which  depend  upon  the  initial  and  final  atomic  configurations  and 
g(ni;e'X')  is  a  radial  transition  integral 

g(ni;€‘i')  ■  r^R(n.2jr)  R(e'i|r)dr  (A-17) 

o 

Although  exact  analytical  expressions  are  available  for  the  matrix 

elements  describing  the  photoionization  of  the  excited  states  of  hydrogen 
(A-9) 

(Gordon  ),  they  are  complicated  and  have  not  yet  been  given  a  numerical 
form  except  for  the  lowest  lying  levels.  Giovanelli^^  employs  an 
approxiisate  formula  due  to  Gaunt for  the  photoionization  cross 


section  of  level  p  at  a  frequency  v 

4  10 


which  in  numerical  form  is 


-  .  .  2.82x10^®  .  2.  . 

QyCp)  =  - -  (cn  )  (A- 19) 

P  V 

This  choice  has  also  been  made  by  Elwert^^  and  Ivanov-Kholodnyi 
et  but  it  is  of  questionable  accuracy.  Recently,  Burgess 

has  presented  tables  useful  for  low  velocities  of  the  ejected  electron, 
from  which  more  accurate  photo ionization  cross  sections  may  be  obtained 
for  a  hydrogen  plasma.  His  calculations  distinguish  between  the  various 
sub-levels  and  apply  to  all  levels  with  p  S  12. 


Using  his  calculations,  the  radiative  recombination  coefficients  may 
be  calculated  according  to  (A- 12)  and  Burgess  has  given  values  of  oe(nX|T^) 
for  electron  temperatures  of  10,000°K  and  20,000^K.  He  do  not  wish  to 
distinguish  between  the  individual  sub-levels  (except  possibly  in  the 
latest  stages  of  the  plasma  decay).  For  a  hydrogenlc  plasma  of  nuclear 
charge  Z,  after  suitably  averaging  over  the  sub-levels,  we  obtain  the 
formula (Seaton 

a(p,Z|T^)  =  DZ  Xp^^^  Sp(X)  (A- 20) 


(A- 21) 
(A- 22) 

(A- 23) 


37 


with 


(x)  exp  (x'i  El(x) 


(A-2A  ) 


S*"  (x)  0.1728  x‘^^  f'"  (u 


.  l)-^»(u 


^  t) 


8^“^  (X)  -  0.0A88  X 


r"  (o 

^  o 


,  l)-''’(o^ 


A  15.788  „2  . 

^  “  k¥“  — T"  = 

pep 


u> 


-  (A-25) 

t  j  u  +  l)e  ’^^clu  (A-26) 

-V 

£ -  dv  (A-27) 

V 


Seatoii^^  has  published  tables  oi  xS^^^(x),  xS^^^.^x)  and  xS^“\x),  iroro 
which  the  radiative  reconibinatlun  cueiiicients  a(p,2|T^)  may  be  calculated 
readily.  Values  obtained  lor  hydrogen  (?.  1)  in  this  way  have  been  provided 
by  Bates  and  are  shown  in  Table  A-5.  Seaton  also  presents  tables  similar 
to  those  for  the  xS'  ''(x)  Irom  whicli  the  total  recombination  coefficient 
to  all  levels 


a 


t 


may  be  calculated  directly. 


(A -28 ) 


For  mure  complex  systems  the  only  method  generally  applicable  is  an 
extension  of  the  Coulomb  approximation.  Tliis  yields  formulas  identical  to 
(A-16)  and  (A-17)  except  that  ii  is  replaced  by  an  effective  (|Uanlum  number 
n*  z/i/  X  being  the  residual  charge  on  the  atom  when  the  electron 

has  been  ejected. 
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TABLE  A-5 


I 
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The  indicles  give  the  poeers  of  10  by  which  the  entries  aust  be  nultiplied. 


According  to  Burgess  and  Seaton 


lA-16) 


g(n*r;t-V)  .  ^  )  cos  n  [n*  +  ,i(e')  +  X(n*i;e'X')]  (A-29) 

r(n\e) 


wlie  re 


G(n*i;€' ^') 


(A-30) 


and 


Iffc  n* 


>* 
fc 


2 


(i(e')  is  the  continuum  phase  shift  obtained  by  extrapolation  of  the  quantum 
defects  energy  levels  of  the  bound  states. 


n* 


and 


V  / 

^££'  ^  n*  ■  .2 

n* 


(n*  >  f  2) 


(A-32) 


(A-33) 


The  remaining  functions  G^^i(n*),  v^^,(n*)  and  coefficients 

' *  ^XX'  been  tabulated  by  Burgess  and  Seaton  for  1  ^  n*  £  12 

and  f,  !'  £  3.  For  values  of  n*  <  i  +  2,  additional  tables  are  presented 
which  permit  the  calculation  of  g(n*X;£'X')  for  transitions  from  specific 
bound  states  to  the  continuum. 


As  an  example  of  the  application  of  Burgess  and  Seaton' s^^ 
formulas  and  tables  to  the  treatment  of  complex  atomic  systems,  they  have 
been  used  to  calculate  the  photo ionization  cross  sections  and  radiative 
recombination  coefficients  for  some  of  the  observed  states  of  aluminum. 
The  ground  state  configuration  of  Al  1  is  3s^3p  and  the  observed  excited 
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2 

states,  with  the  exception  of  one  3s3p  configuration,  are  configurations 
of  the  form  3s^n£  where  3  <  n  <  11  and  0<f<3(cf.  C.  Moore  ^ .  The 
latter  configurations  consist  of  a  single  (nf)  electron  outside  a  closed 
shell  for  which  the  C^,  assume  the  particularly  simple  form  (selection 
rules  on  the  azimuthal  quantum  number  allow  only  transitions  where 
i'  --  f  ±  l)i 

_  f  + 1  —  SL  g  ^  #  V 

^f+1  2i+l  ’  ^i-1  2f+l  (A-34) 

Table  A-6  gives  the  Burgess  and  Seaton  coefficients  for  the  allowed 
transitions  In  A1 1 .  These  coefficients  were  used  to  calculate  the  func¬ 
tional  parameters  ^^^ •(”*)»  individual  transitions. 

The  results,  together  with  the  ionization  potentials  and  quantum  defects 

I 

of  the  bound  states  «*■«  listed  in  Table  A-7.  The  ^(n‘'  ,f),  also 

given  In  Table  A-7,  were  determined  from  the  curves  of  Figure  A-1  which 
shows  plots  of  the  continuous  phase  shifts  for  each(n.C)  series. 

For  small  positive  electron  energies  e',  the  phase  shift  may  be  taken 
as  a  linear  function  of  e*.  The  extrapolated  quantum  defects  determined 
from  the  curves  of  Figure  A-1,  lead  to  the  following  approximate  expres¬ 
sions  for  u(e'): 


H(€')  =  1.755 

(ns  series) 

n(£')  =  1.264  -  0.065e' 

(np  series) 

(A-35) 

H(€')  =  0.985  +  3.3c’ 

(nd  series) 

H(c')  -  0.044  +  0.26c' 

(nf  series) 
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•4 
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Of^ 
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2.452 

1.309 
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1.163 

1.901 

1.163 

1.779 

(0.660) 

(3.100) 

(0.732) 

(2.752) 
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2.342 
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2 

fhotoionizatlon  eroit  ••ctions  were  calculated  for  all  but  the  3§3p 
and  n£  levels,  for  ejected  electron  energies  c*  lying  in  the  interval 
(ftydbergs)  0  5  e'-  0.5  for  the  nd  levels,  0  5  c'l.  0.3  for  the  ns  levels 
and  0  5.  e'l.  0.2  for  the  np  eeries.  For  electron  energies  in  excess  of 
the  maxiDua  values  indicated,  it  was  found  that  the  aj^roxisMte  foraulaa 
led  to  incorrect  results.  This  SMy  in  part  be  due  to  the  aseusiption  of 
a  linear  dependence  for  the  however,  it  should  be  noted  that  the 

approxinations  Inherent  in  the  derivation  of  (A-29)  Unit  the  validity 
of  the  calculations  to  small  ejedted  electron  energies. 

The  photo ionization  cross  sections  obtained  are  listed  in  Table  A-8. 

2  2 

The  3s  3d  and  3s  4d  cross  sections  were  calculated  under  the  assusqption 
that  n*  j(+2  since  explicit  expressions  are  not  tabulated  by  Burgess  and 
Seaton for  nd  “•  c'p  snd  nd  *•  c'f  transitions  when  n*  <  ^+2.  The 
resulting  cross  sections  are  therefore  of  lesser  accuracy  than  the  remain¬ 
ing  cross  sections  shown. 


The  radiative  recoad>ination  coefficient  amy  be  calculated  from  Bq.  (A-12) 
once  the  photoionization  croas  sections  arc  known  as  functions  of  the  ejected 
electron  energy.  Thus,  in  terms  of  electron  energy  for  a  thin  plasma 


a(n/|T.)  .  ^  d.  (A-36) 

Instead  of  evaluating  the  integral  of  (A-36),  the  recombination  coefficient 
can  be  determined  approximately  by  use  of  the  ei^ression 


4r 


a(ae|T^)  *  V  (lU) 


(A-37) 
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where  (ni)  the  crosa  section  for  radiative  recombination 
detailed  balancing 


Qg  (ni)  »  (2i+l) 


2  ^ 

me  e 


and  V  Is  the  average  electron  speed 


%  (n/) 


Substituting  (A-39)  and  (AOS)  in  (A-37),  we  find 


ate/ll,)  .  (5J5-)  (n/) 

^  IflC 


6  2 

-  2.046x10®  (2i+l)  (1  +€)  ~i/2 


is  given  by 


(AOS) 


(A09) 


(A-40) 


The  radiative  recosdiination  coefficients  calculated  from  (A-40)  are 
given  in  Table  A-9  for  electron  tesiperatures  of  6190”K  and  24»760”K. 


TABLE  A-9 


RAOUTZVE  REOnaiHATIOM  COEFFICIENTS 
F(»  A1  I  0£(n£|T  )tn  cm  / sec) 


T  -  6190®K  T  -  24,760®K 

0  6 


The  brackets  give  the  power  of  10  by  which  the 
entries  must  ',0  multiplied. 
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A.  3  Electron  Collision  Induced  Transitions  B«tW8en  Bound  Sft«i 


Let  (Ptq)  be  the  cross  section  for  the  excitation  to  the  q  state 
by  the  impact  of  an  electron  of  energy  e  with  an  atom  in  the  p  state. 

The  rate  of  depopulation  of  the  p^^  state  by  electron  impact  is 


n 

e 


n(p) 


^C(p.r) 


r  >  p 


(A-41) 


where  the  collision  probability  C(p,r)  depends  upon  the  electron  temp¬ 
erature  T  and  is  obtained  by  integrating  the  product  of  the  electron 
velocity  and  the  cross  section  over  the  Maxwellian  distribution  accord¬ 
ing  to 


C(p,r) 


(KT  ) 
e 


-3/2 


1 


r 


00 

q^Cp.r) 

-  I 

P 


c  exp(-€/kT^)  d€  (A-A2) 


The  probability  of  de-excitation  to  a  lower  level  is  similarly 


n 

e 


C(p.q) 


q 


(A-43) 


Q 


The  direct  and  Inverse  collision  probabilities  are  related  through  the 
principle  of  detailed  balancing: 


excitation  and  de-excitation  transitions  induced  by  electron-atoa 

(a_4v 

collisions  (c.f.  Allen'  %  none  of  which  is  entirely  satsifactory. 
Those  quantum  mechanical  treatMnts  in  which  some  reliability  may  be 
placed  (e.g.  Mott  and  Massey Bates  et  al^^”^^\  Milford  and  co- 
workers^^' McCrea  and  Me  Klrgan^^’^^\  and  Me  Cerrol^^'^^^  h«ve 


been  based  on  the  Born  approximation,  and  therefore,  are  valid  only 
for  high  electron  impact  energies. 


A  Blnpllfled  formila  daacrlblng  the  croaa-aectlon  for  axeitatlon  froH  a 

/A»A\ 

q  level  to  a  p  level  Induced  by  the  form  (c.f.  Allen^  ') 

QJq.p)  -  1.44  X  10“xa^  (ch)^f(q.p)  -i—  -  i)  (c«^  (A-47) 

q>p 

where  e  denotes  the  lispacting  electron  energy,  and  e  _  is  the  ex> 

H»P 

citation  energy,  both  expresastd  in  gy^berga  (e  ■  c  >  cl.  f<q,p)  reprcscnta 

•qtP  q  P, 

the  oscillator  strength  for  absorption,  and  is  related  to  the  spontaneous 
transition  probability  A(Ptq)  by 

3  u> 

f(q.p)  “  -^2  2  nf"  (A-48) 

8k  e  V  q 

Equations  (A-47)  and  (A-48)  euty  be  used  to  calculate  collisional  tranaition 
probabilities  for  the  atomic  systems  for  which  the  Spontaneous  transition 
probabilities  are  given  in  Tables  A-1,  A-2,  and  A-3.  If  Q^(p,q)  denotes 
the  corresponding  collisional  de-excitation  cross  section,  then  (A-47)  ‘ 
and  (A-48),  with  the  aid  of  the  following  expression  (c.f.  Allen^^*^^): 

Vp  ^c^**’**^  “  Vq  ^e^**'**^  (A-49) 

Wpi(p.q)  -  «^f(q.p)  (A-50)' 


yields 


Q^(p.q) 


5.78  X  10“  xa?  (ch) 


u 

-£ 

w 

P 


(A-51) 


Ixposssimis  sisdlar  to  (A-47)  and  (A-Sl)  with  values  of  A(p,q  )  obtained 
by  the-  forsHilas  of  Hansel  and  Pekeris^^*^^  for  hydrogen  were  employed 
by  Ciovanelli^*"^®^  and  Ivanov-Kholodnyi  et  al^^*^^^. 


S2 


Perhaps  the  nost  accurate  method  of  computing  colllslonal  transition 

(A-23i 

probabilities  is  the  classical  method  of  Gryzlnski  \  which  represents 
an  improvement  over  the  original  formulation  of  Thomson^^'^^\  in  that 
it  takes  account  of  the  fact  that  the  atomic  electrons  are  moving.  The 
cross* sections  so  desired  fall  off  too  rapidly  at  high  electron  energies, 
as  fc  ^  instead  of  the  correct  (c  ^  In  e),  but  their  accuracy  j^obebly 
suffices,  the  effect  of  excitation  and  de*excitation  not  being  a  pri* 
raary  one. 


Considering  in  sosm  detail  the  classical  transfer  of  axmentum  during 
a  coulomb  collision  between  an  incident  electron  of  energy  e  and  an  atomic  elec¬ 
tron  having  binding  energy  Gryzlnski  derived  the  following  approxismte  ex- 
preusior  for  the  differential  cross-section  corresponding  to  a 

loss  in  energy  of  the  incident  electron  of  magnitude  Ac ; 


o(<^)  « 


s-4  f"-*— 

I?  ♦  'iv/ 


c(Ae) 


^  €  -  € 


3 


^  >  e 


(A-52) 

The  total  cross  section  for  a  collision  involving  an  energy  loss  subject 

to  the  restriction  €  <  c  <  e  is  then  given  by 

p  -  ~  q 

P 

■  |a(^5e)d(^)  •  Q(€  )  -  Q(c  )  (A-53) 

€  p  q  J  p  q 

*q 

where  QCe^)  the  cross  section  for  a  collision  with  electron  loss  greater 

than  e.  Is  defined  by  Equation  (A-^*)).  Thus  idien  Ae  >  e  ,  as  in  tiM  case 
P  P 

for  hydrogenic  orbits » 


1/2 


53 


1 

3  «  ^ 


Q(«p) 


3.53  K  IQ 


•16 


^  \  V  ’  V  / 


V 


[i“  +  ;*  6  -  ^*)  ■  ] 


*-n1/2 


(cm  )  12€  2  « 

P 


(A- 54) 


The  collision  trsnsition  pri^ability  is  given  by 

.1/2 


«P.-0-(lfei)  V  ‘  (A- 


55) 


Inserting  *p)  leads  to  the  following  collision  prob¬ 

abilities  for  hydrogen,  where  c  >  l/p3  and  c  ■  l/q^ 

C(P.0  ■  x®4^2  ■*§”  /*€•«?  {“^*^'^'^2  ^  cm^»ec  ^ 

e  q  I  qp  T. 


(A- 56) 


where 


and 


<«l^  -  P^ 


y  ■  «p  g 


/i 


j  [2  -(g+  y(l  +  4g)]  (y  -  l)^'*(g  +1)  v  <  1  ♦  a 

.1/2  ,  ^  .3/2 
g  (g  +  y) 

(A-5« 

^  j  c-3g  ♦  y(3  ♦  4g)]  y'^^ 

.3/2 


,1/2 


(A- 57) 


f 

»r»' 


(g  +  y)' 


y  >  1  +  g 
(A-59) 

C(p.q)  has  been  calculated  for  hydrogen  by  Bates  using  equation 
(A- 56).  The  results  are  given  In  Table  (A>10)  for  p  &  q  less  than 
or  equal  to  10. 
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For  non -hydrogenie  systeins  having  known  energy  levels  the  more 
general  expressions  of  (A-52)  and  (A-53)  must  be  used  for  the  calcula¬ 
tion  of  collision  probabilities. 


i 

1 

1 
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.4  Electron  Collision  Induced  Bound-Free  Trantttlotii 


Let  q(p,e)  dc  be  the  cross  section  for  the  formation  of  electrons 
with  energy  In  the  range  c,  c  +  de  by  collision  of  an  electron  having 
energy  E  with  an  atom  In  the  p  state.  The  total  Ionisation  cross  section 
Is 

E-Ip 

Qe(p)  -  J  qE(P,«)  d«  (A-60) 

o 

and  the  probability  of  Ionization  by  collision  Is 


D(P) 


(kT^) 


-3/2 


I 


qg(p)E 


exp(-E/kTg) 


(A-61) 


The  rate  of  depopulation  of  the  (p)  state  due  to  electron  Impact  Ioniza¬ 
tion  Is 

n(p)  D(p)  (A-62) 

The  Inverse  process  Is  three-body  recombination: 


X'*’  +  6  +  0-*  X(p)  +  e. 


(A-63) 


The  efficiency  of  the  three-body  process  may  be  written  In  the  form 


F(P) 


«2  exp(-«j/kT^)  exp(-«2/Wg) 


X  Y(p|«j,e2)dej  d«2 


(A-64) 


where  y(p|<]^*2)  probability  that  an  atom  Is  found  In  the  (p)  state 

as  a  result  of  the  collision  with  a  positive  Ion  and  two  electrons  with 
enargles  and  C2  respectively.  The  probability  is  related  to  the  In¬ 
verse  process  of  electron  iaq>act  ionization  by  detailed  balancing.  Thus 
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ISmau. 

<Ip  +  «i  +  *2)  ■  “5 -  *1*2  Y(pI*i.  *2) 

h  u)p 

Where  ou^  Is  the  statistical  weight  of  the  continuum  electron  and 

B  ■  I  +  c,  +  e,  .  (A-66) 

p  1  2 

The  rate  of  population  of  level  (p)  by  the  three-body  recombination 
process  (Equation  A-63)  Is  then 

n^  nCX"^)  P(p)  (A-67) 

Alternatively,  the  probability  of  three-body  recombination  F(p)  may 
be  obtained  by  a  slispler  method  noted  by  Ivanov-Kholodnyl,  et  al^^ 

They  remark  that  In  thermodynamic  equilibrium,  the  number  of  particles 
leaving  a  given  quantum  level  via  any  specific  process  Is  equal  to  the 
number  of  particles  entering  this  state  via  the  Inverse  process.  It 
follows  that 

F(p)  -  9(p)  D(p)  (A-68) 
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vrtiere 


,2n«k  T  v3/2 

cp(p)  .  p2  expdp/kT^) 


(A-69) 


As  in  the  case  of  excitation,  the  cross  sections  decrease  too  rapidly 
as  the  energy  of  Impact  Increases.  Since  Ionization  and  Its  Inverse 
are  the  controlling  processes,  the  error  may  be  serious  and  we  have 
therefore  given  some  attention  to  the  problem  of  how  we  may  obtain  Im¬ 
proved  accuracy,  whilst  still  retaining  the  universal  applicability  and 
the  simplicity  of  the  Gryzlnskl  method.  A  new  universal  method  that 
has  resulted  Is  described  below. 


Again  for  reasons  similar  to  those  mentioned  above  for  bound- 
bound  collision  probabilities,  the  most  accurate  method  for  computing 
electron  impact  ionization  probabilities,  which  is  of  universal  ap¬ 
plicability  is  that  of  Gryzinski^^”^^\  Giovanelli^^"^®^  and 
Elwert^^"^^^  used  the  classical  Thomson  fonoula  which  can  be  expected 
to  give  somewhat  less  accurate  results  than  that  of  Grysinski. 
Ivanov-Kholodnyi  et  al^^~^^^  employed  approximate  expressions  simi¬ 
lar  to  those  for  the  collisional  excitation  cross-sections. 


O 


or 


The  classical  Thomson  cross-section  (cf.  Thomson 
Fowler for  the  ionisation  of  an  atom  having  n^  valence  elec¬ 
trons  by  an  electron  with  energy  € ,  with  the  corresponding  ejection 
of  an  electron  from  a  p-state  into  an  energy  interval  t,  C  +  d^  is 

2 

Anna 

Q.  (P.C)  di  =  - - - H 

‘  ‘<5  -  %>' 

The  Thomson  total  ionisation  cross-section  is  then 


(P) 


Anna 

ja 

c 


-  ■■  ■ 
(C  ^  Cp)2 


=  An  a  n 
o  o 


“  s) 

c«2 

P 


(A-71) 


a^  is  the  radius  of  the  first  Bohn  orbit  in  hydrogen. 


Ai  pointed  out  in  the  previous  section,  the  Grysinskl  formulss 
represent  sn  laprovesient  over  the  ThooMon  derivation  in  that  accouht 
is  taken  of  the  motion  of  the  bound  atomic  electrons.  According  to 
Grysinski  the  total  ionisation  cross-section  for  an  electron-electron 
collision  in  which  the  impacting  electron  suffers  an  energy  loss  AE 


greater  than  the  threshold  ionisation  energy  s  is  given  by 

P 


f^roax 

q(c  )  m  }  9(At)  d  (AE) 

o  • 


(A- 72) 


where  AE^^  is  the  maximum  permissible  transfer  of  energy  consistent 
with  conservation  of  momentum  and  energy,  and  a  (AE)  is  the  differ¬ 
ential  cross-section  described  by  Equation  (A-52).  (The  expression 
for  Q(c)  given  in  (A-54)  was  obtained  from  the  integration  indicated 
by  (A-72)).  Thus  the  croes-section  for  ionisation  from  the  p  level 
is  Just  (A-S4).  The  ionisation  transition  probability  for  a  Haxwel- 
lian  electron  energy  distribution  is  given  by 


D(p|T^) 


(A-  73) 


Making  the  substitutions  x  •  e/c^  we  find  from  (A-S4)  and  (A-73)  that 


(cm^/sec)  X  >  2 


:5 


X  <  2 


In  the  ceee  of  hydrogen,  «  ■  1/p^  *0  thpt 


r¥  f  (50’' -  X }  dx 


“‘"Iv 


(A- 75) 


i  I"  »-■’  {'*'^^1^  «><>'»)*■  x>J 


Values  of  D(p|  in  hydrogen  (provided  by  Bates)  are  given  in  Table 
A-11  for  p  less  chan  or  equal  to  10.  For  the  more  general  ease,  Equa¬ 
tion  (A- 74)  oust  be  used  to  compute  D(nl  |  T^). 


As  in  the  case  of  excitation,  the  cross-sections  decrease  too 
rapidly  as  the  energy  of  impact  increases.  Since  ionisation  and  its 
inverse  are  the  controlling  processes  at  the  higher  energy  levels, 
the  error  may  be  serious  and  ve  have  therefore  given  scam  attention 
to  Che  prdilem  of  how  to  obtain  improved  accuracy,  %rhilst  still  re¬ 
taining  the  universal  applicability  and  the  simplicity  of  the  Grysin- 
ski  sKthod. 


Bewricing  Equation  (A-71)  for  the  nicmson  ionisation  crMS-section 
in  the  form 


«) 


Asa 


2  1 


0-?) 


(A-?6) 


and  introducing  the  diswnsicmless  parameter  X  ■  s/l 


TABLE  A-11 


I 


(A-47)  becomes 


/  2  1^,  1  ^ 


(A-77) 


^e  rlght-hend  side  of  (A-77)  is  e  universel  function  of  X,  independ¬ 
ent  of  the  perticuler  ionising  process. 


It  is  known  from  the  asymptotic  form  of  the  Born  epproximetion 
(cf.  Mott  end  Massey that  for  large  values  of  X,  behaves 
asymptotically  according  to 

Qj  -  X*^  In  X  (A-78) 

so  that  (A-77)  is  incorrect  in  the  limit  of  high  impact  energies. 

The  modification  of  (A-71)  by  Grysinski also  has  this  defect. 


Ue  shall  attempt  to  improve  (A-77)  by  introducing  a  correcting 
factor  f  (X)  such  that 


I  Q<(U)  . 


(A- 79) 


«  4«  a/  g(X) 


and  determine  f  (X)  by  appeal  to  experimental  data. 

In  figure  (A-2),  all  the  available  experiswntal  data^^”^^ 

on  the  electron  impact  single  ionisation  of  atoms  has  been 
2 

used  to  determine  I  *•  *  function  of  X,  adopting  for  n^  the 

nusiber  of  electrons  in  the  valence  shell.  If  f  (X)  exists  such  that 
(A-79)is  true,  this  procedure  would  produce  a  single  curve.  Instead, 
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there  results  distinct  curves  which  are  similar  In  shape  but  have  an 
appreciable  spread  In  magnitude  of  the  croas-scctlons.  This  spread 
represents  the  possible  error  that  would  be  a  consequence  of  the  use 
of  the  mean  of  the  curves  of  Figure  (A-2)  as  a  universal  function. 

Ve  may  reduce  this  factor  by  taking  partial  account  of  the  elec¬ 
tron  structure  of  the  atomic  systems.  It  being  an  oversimplification 
to  adopt  for  n  almply  the  nusd>er  of  valence  electrons  n^.  The  signifi¬ 
cant  quantity  Is  actually  the  equivalent  number  of  dispersion  elec¬ 
trons  for  transitions  into  the  contlnum. 


To  determine  n,  we  may  proceed  in  two  alternative  ways.  Ue  may 
choose  n  to  ensure  all  the  curves  have  the  same  magnitude  at  large 
energies  or  we  may  observe  from  (A- 71)  that  the  maximum  value  of 
is  given  by 


Q 


max 

1 


(A- 80) 


and  determine  n  from 


n 


Q 


max 

«  a  ^ 
o 


(A-81) 


The  latter  procedure  is  the  simpler  and  we  adopt  It.  It  leads  to  the 
curves  shown  In  Figure  (A-3),  all  of  the  normalised  experimental  data 
for  Ionisation  of  neutral  atoms  lies  within  the  shaded  area.  No 
curve  now  departs  from  the  swan  of  all  the  curves  by  a  factor  of  awre 
than  about  lOX  (cf  shaded  area  of  Figure  A-3)  and  we  shall  select  the 
mean  as  a  universal  curve.  The  normalised  curve  for  Ionisation  of 


69 


0IDA40-I370 
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Figure  A-3.  Uaivereal  Curve  for  Electron  Impact  Ionization 

Cross  Sections 


of  He^  Ion  falls  Just  outside  of  this  area  and  is  also  shovm  In  Figure 
A-3. 

The  corresponding  values  of  n  are  given  In  Table  (A-12a).  We  see 
that  there  Is  a  pattern  to  n  which  allows  us  to  predict  It  with  some 
reliability  at  least  for  the  lighter  elaaents.  Thus,  since  n  appears  to 
be  an  Increasing  function  of  atomic  number,  we  can  assign  Integer  values 
of  n  to  appropriate  ranges  of  the  atomic  number  as  sho«m  In  Table  (A>12b). 


From  the  universal  curve  of  Figure  (A-3),  we  can  proceed  to  the  de¬ 
termination  of  a  universal  Ionization  probability  curve  as  a  function  of 
electron  temperature  T^.  Thus  from  Equation  (A-S8)  we  find 


D(T  ) 
e 


(kT^) 


-3/2 


f"  Q^(X,I)  dX 


(A-82) 


Using  Equation  (A- 77) 

D(T  )  =  (kT  )"^^"  4n  na  ^  T*  Xg(X)e"^^'‘’^e  dX  (A-83) 

e  \**™/  *  ®  Ji 

where  T^  =  T^/1  Is  the  electron  temperature  measured  In  units  of  the 
Ionization  potential.  Finally 

0(T^)^^3/2  ^  4«  dX  (A-84) 

The  right  hand  side  of  which  Is  a  function  of  the  electron  temperature 
T^  (mly  and  may  be  computed  once  for  all  atomic  systeois  for  any  given 
value  of  T^,  using  the  universal  curve  g(X)  given  In  Fisure  (A-3). 
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X4BLI  A-12a 


EQUIVALBMT  NUMBER  OP  DISRISXON  BUCTBONS, 
n  OETBIHIMR)  Wi  EQUATION  (A-ll) 


Element 

H 

He 

1^ 

B 

mm 

Kr 

Xe 

Ng 

% 

1 

2 

— 

6 

6 

Di 

2 

n 

0.8 

1.34 

1.75 

2.24 

4.20 

6.05 

7.14 

3.74 

1 


TABLE  A-12b 
ESTIMAnD  VALU  OF  n 


n 

H  -  Be 

1 

B  -  Ne 

2 

Me  •  A 

4 

K  •  Kr 

6 

Ce  -  Xe 

7 

73 


The  function 


DT 

. ^  la  shown  as  a  function  of  T'  In  Figure  (A-4). 

n  e 

The  Ionisation  probability  for  any  element  for  any  temperature  may  be 

read  off  the  graph,  the  appropriate  abscissa  being  obtained  by  dividing 

T^  by  the  Imlsatlon  potential  1. 


1 

I 

! 
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94^  m,mo  4o,ooo 

r/  *  T^/I 


Figure  A-4.  IkiiverMl  XonlMtlon  ProlMiblllty  Curve 
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A. 5  Heavy  Particle  Colltatons 


Consider  now  Cransltlona  produced  by  collisions  with  heavy  particles 
which  we  label  Y^.  Let  (p«r)  be  the  cross  section  for  the  excitation 
to  the  r^^  state  by  the  iaipact  of  a  heavy  particle  with  an  atom 
X  in  the  p^^  state,  the  energy  of  reletive  sntion  being  e.  Hie  rate  of 
depopulation  of  the  p^^  state  by  heavy  particle,  impact  is 

n(p)  5"  n(Y  )  ^  P  (p.r)  (A-85) 

T  ^  r>p  ^ 

where  the  collision  probability  Pj.(p,r)  is  obtained  by  multiplying  the 
cross  section  by  the  relative  velocity  of  the  collision  and  integrating 
over  ion  velocity  distribution  f(v)  to  give 


P^(P.r)  ■  J  **''  (A-86) 

If  f(v)  is  the  Mexwellien  velocity  distribution  corresponding  to  temp> 
ereture  T, 


00 


P^(P.r) 


-  ~  (W) 

xp- 


'3/2 


I 


QjCp.r) 


e  exp(>e/kT)  de 


I  - 
r 


(A-87) 


where  is  the  reduced  mess. 


The  probebility  of  de-excitetion  to  a  lower  level  is  siisilerly 
given  by 

n(p)  n(Y^)  Y 

1  qCp 


76 


The  rates  of  population  are 


^  n(Yj^)  ^  n(q)  (q,p)  (A-89) 

1  q<p 


i  r>p 


The  ionisation  end  recombination  involving  collision  with  ions  may 
be  formulated  in  a  manner  similar  to  that  for  electron  collisions.  The 
rate  of  depopulation  by  ion  impact  induced  ionisation  may  be  written 


n(p)  {n(x'*’)  Q(p)  +  Z  "(Yj)  Qj(P))  <^-91) 


J 

OD 


Q(P)  •  -"-"T  /  Q»iP)  exp  (-B/2kT)  dE 


(A-92) 


Q  (p)  is  Che  total  Ionisation  cross  section  of  the  state  (p)  by  Ions  of 
G 

energy  B. 

The  rate  of  population  by  three>body  recon^lnatlon 

+  Yj  +  e  -  X(p)  +  Yj  (A-93) 

may  be  written 

nCx”*”)  n(Yj)  B(p)  (A-94) 

where  R(p)  Is  an  expression  exactly  similar  to  Equation  (A-64) 
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since  for  allowed  traneltlona  the  croas  sect  Iona  for  the  impact  of 
an  electron  and  a  singly  charged  positive  Ion  of  the  SMie  velocity  are 
of  much  the  same  alse,  excitation  and  Ionization  by  Ion  liqpact  Is  usually 
very  Inefficient  compared  to  electron  Impact.  For  exaa^le,  three*body 
recombination 

+  X"*"  +  e  -►  X‘*‘  +  X  (A-95) 


will  be  negligible  compared  to 

x''”  +  e  +  e-*X-«'e 


(A>96) 


E;:ceptions  may  occur  when  the  energy  defect  of  the  process  Is 
small  as  In 


H(n,i)  -►  H*  +  H(n,/') 


(A-97) 


and  for  transitions  between  pairs  of  highly  excited  states.  The  effects 
of  Ion  Impacts  Is  to  strengthen  the  argument  that  a  BoltssMnn  equili¬ 
brium  prevails  for  the  excited  states.  The  position  regarding  charge 
transfer  reactions  Is  analogous,  charge  transfer  being  rapid  at  low 
velocities  only  when  the  energy  defect  Is  low. 


We  merely  wish  to  provde  a  few  numbers  In  order  to  show  quantitatively 
that,  apart  from  re-distributing  processes,  (and  transitions  between 
pairs  of  highly  excited  states),  electron  collisions  arc  more  efficient 
than  heavy  particle  collisions. 
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(A- 34) 

Seaton'  has  calculated  the  rate  coefficient  for  the  re-distribution 

process 

h'*'  +  H(2s)  h'*'  H(2p)  (A-98) 

o  •4>  3 

at  10,000  K,  and  he  obtains  the  very  large  value  of  about  5  x  10  cm 

sec  Re-distribution  may  also  occur  by  charge  transfer 

h'*’  +  H(2s)-r  H(2p)  +  h'*'  (A-99) 

but  Boyd  and  Oalgarno^^  have  shown  that  (A-99)  is  much  less  rapid 
than  (A-98) . 

Excitation  and  ionisation  proceed  much  more  slowly  and  may  be 
neglected.  Thus,  the  maximum  values  of  the  cross  sections  for  the 
direct  processes 

h'*’  +  H(ls)  -►  +  H(2s)  (A- 100) 

h'*’  +  H(ls)  -►  +  H(2p)  (A-101) 

-17  2 

are  respectively  about  5  x  10  cm  at  about  20  kev  impact  energy, 

(Bates^^  and  1  x  10  cm^  at  about  30  kev  (Bates^^  and 

the  cross  sections  are  decreasing  rapidly  with  decreasing  energy. 

Similarly,  it  swy  be  seen  from  the  calculations  of  Bates  and 
(A*38  A* 39) 

Dalgarno'  *  '  that  excitation  accompanying  charge  transfer  such 

as 

h'*’  ♦  H(ls)  H(2p)  (A-102) 


.•9 


Is  extrsmely  Inefficient  for  temperatures  below  several  kev. 

Thus,  for  a  hydrogen  plasma  or  other  singly  Ionized  plasma  of  a 
single  component,  processes  due  to  Ion  Impact  may  be  Ignored.  When  more 
than  one  component  exists  In  a  plasma,  however.  It  may  be  necessary  to 
consider  heavy  particle  processes,  for  charge  transfer  processes  are 
unique  In  coupling  together  the  two  components.  This  fact  may  be  seen 
from  Equation  (A- 106),  by  noting  that  If  the  heavy  body  coefficients 
?!,  Q^,  and  are  set  equal  to  zero,  all  the  terms  containing  n(Y^), 
the  1. umber  density  of  other  Ions  vanish.  Account  must  be  taken  not 
only  of  ordinary  charge  transfer 

+  y'*'  X"*"  +  y"^  (A- 103) 

but  also  of  radiative  charge  transfer 

x"””  +  y'*'  -•  x'*’  +  y"*^  +  hv  (A- 104) 

which  may  be  more  rapid  than  the  previous  reaction. 

An  Investigation  of  the  possible  rates  of  these  cherge  transfer 
reactions,  one  of  the  participants  of  which  Is  an  Ionic  system  in  Its 
ground  state.  Is  necessary,  and  as  an  example  of  the  latter  process 

He”*^  H  He'*'  He'*'  hv  (A- 105) 

may  be  treated  with  high  accuracy. 

Special  features  also  arise  In  the  case  of  a  single  cosqwnent  plasma 
In  which  the  cosqwnent  may  exist  In  more  than  one  state  of  Ionization. 
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A.6  Th«  Basic  Eauatlons 


The  increase  in  population  per  second  of  a  particular  level  (p)  Is 
given  by  the  difference  between  the  rates  of  population  and  depopulation 
that  we  have  written  in  the  earlier  sections.  Thus 

r  ^  A(p,q)^l  +  X(p,q)^  ^  X(p,r)A(r ,p)| 

+  Ue  C(p,r)  +  Z  C(p,q)  +  D(p)J 
+  Z  n(Y^)  j^Zp^(p,r)  +  EPj(p,q)  +  Qj^(p)J 
+  3(p)|  n(p) 

•  Z|A(r,p)  [l  +  X  (r,p)]  +  n^  C(r,p)  +  ^ 

-(qo%)  A(pq)  +  n^  C(q^p)  n(Y^)Pj^(q,p)|  n(q) 

cKp)  XCn”**)  x"*"  y^Cp) 

*  \  \  "iCY*)  nCx”^)  Rj(p)  (A-106) 

We  arrive  at  a  set  of  coupled,  differential  equations  of  the  form 

MeI  +  «  „(p)  .  y  .  n  +  b  (A-107) 
dt  pp  Z-  pj  J  p  ' 

J  ^  ? 

where  the  a's  and  b's  can  be  identified  by  coaiparison  with  the 
preceding  equation. 
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APPENDIX  B 
A.  Naqvl 

ATOiXC  tlUMSniON  PROBABILITIES 
B.l  INTRODUCTION 

Considar  two  anargy  lavala,  1  and  2,  of  any  aton,  and  lat  thalr 
«M«tloa,  E^  and  £2,  ba  auah  that  E2  >  E^.  Tranaitiona  fra*  2  to  1  raault 
In  anlaalon,  and  tha  ravaraa  tranaltlona,  fro*  1  to  2,  ara,  aaong  othar 
proeaaaaa,  cauaad  by  tha  abaorptlon  of  radiation,  of  fraquancy  V2j^  given 
by 

Ez  “  ^  \  ^  (B-l) 


*1  m\ 

Tha  aaMunt  of  anargy  (In  arg  cn*  aac~  )  aalttad  In  a  apactral  line 
la  given  by 

B(X>\)  (B-2) 

3 

wham  n2  la  tha  nuabar  of  atoaw  par  ea  in  laval  2.  A(2,l)  and  B(2,l) 
ara,  raapactlvaly,  tha  probabllltlaa  (la  aac'^)  of  apontanaoua  tranaltlona 
and  tha  tranaltlona  atUulatad  by  tha  praaanca  of  an  axtamal  radiation 
flald.  I2j^  la  tha  Intanalty  of  thla  axtamal  radiation  flald  at  tha 
fraquancy  V2^. 

Similarly  tha  amount  of  anargy  (In  arga  cm’'’  aac*  )  abaorbad  In  a 
apactral  llna  la  glvan  by 
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(B-3) 


«h«r«  B(l,2)  is  th«  probability  of  tho  upward  transition,  from  laval  1  to 
laval  2,  with  tha  absorption  of  radiation  of  fraquancy  V2|^* 

Tha  thraa  transition  probabllltlas  ara  ralatad  to  aaeh  othar  by 
tha  folUnrlng  rolatlons  dua  to  Blnataln. 

(B-*) 


Tt"?" 


(B-5) 


whara  and  §2  ara  tha  statistical  weights  of  tha  two  lavals.  It  Is 
tharafora  sufficient  to  determine  one  of  the  three  transition  probabilities, 
usually  A(2,l). 

Tha  thaory  for  calculating  atoadc  traaaltlon  probiAllltleo  f<nr 

alactrlc  dlpola  radiation  Is  well  known  (Condon  and  niortley^*^\ 

(B>2) 

SlatMr'  Several  reviews  have  been  written  In  recent  yews  also 

(Garstang^^^^,  The  prohabillty  A(2,l)  of  a  apontMoous 

radlatl^  traaaltlon  la  glvm  by 


S6 


(B-6) 


iA«r«  g2  it  tht  ttatlttlcAl  weight  of  level  2,  1  it  the  wavelength  of  the 
epectrel  line  emitted,  S(2,l)  ia  known  et  the  line  strength  end  for  LS 
coupling  Mn  be  written  ee 


%ihere  J*<«  end  f<M)  ere  feetore  which  eriee  ee  e  result  of  integratim 
over  the  angular  pert  of  the  wave  function.  ]f(L)  is  known  as  the  rela¬ 
tive  strength  of  a  line  within  the  nultiplet  to  which  it  belongs  md 
^(M)  is  known  as  the  relative  nultiplet  strength.  Both  of  these  quanti¬ 
ties  can  be  calculated  in  a  closed  fom  for  any  given  line  end  any  given 
nultiplet.  ^(L)  has  been  tebulated  by  Mkite  and  Bliason^^*^^  and  by 
Bussell^ fW  has  been  tebuUted  by  GeUberg^"^  These 

tables  have  also  been  reproduced  by  Allen^*^^  Becently  Rehrlich^^^^^ 
has  derived  general  expresaiene,  baaed  on  Raceh's  nsthod,  fron  which 
nultiplet  strengths  for  traneitlone  net  covered  by  Goldberg  can  be  obtained. 


The  electric-dipole  tr«sition  integral  a  ie  given  by 

00 


idiere  the  P'e  are  r  tines  the  redial  wove  fuMtiena,  and  is  the  greater 
of  the  two  asiauthal  quantun  nunbere  and  ij*  calculation  of  the 
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•btolut*  transition  probabilitiss  thorafero  d^anda  upon  tha  ealculatim 

2 

of  tha  transition  intatral  a  ,  ohidi  in  turn  dopmda  an  appropriata 
ehoiea  of  tha  radial  wava  functions. 

B.2  THE  SADUL  HAVE  FUMCTI(»IS 

Tha  radial  oava  functions,  P(ni:r),  ara  tha  solutions  of  tha  radial 
part  of  tha  non«ralativistic  Schroadingar  aquation, 

+  <•-»> 

whara  E  is  tha  aigan  valua  of  tha  Haailltonian  for  a  givan  stationary  stata 
and  V(r)  is  tha  cantral  fiald  in  uhich  tha  alactron  novas.  Hots  that 
atonic  units  arc  usad  throu^iout  this  raport. 

Tha  iaportaaca  of  a  kaowlodgo  of  tha  radial  nsva  functions  cannot 
ba  ovar  anphasisad.  Tha  radial  nova  function,  in  a  outsail,  glvas  us 
an  insight  into  all  othar  physical  propartias  of  tha  aton.  A  knowladga  of 
tha  radial  uava  functions  is  roquirod  not  only  in  tha  calculation  of  tha 
transition  probabilitiss  but  also  photoionisation  cross  sections,  elastic 
scattaring  cross*saccions,  collision  cross-sactions  for  eraitation  of 
atons,  x*ray  scattarii^  and  nunorous  nolocular  and  solid  state  propartias, 
anong  others. 

Tha  nost  coansaly  usad  nathod  is  tha  self  •consistent  fiald  nathod, 
davalopad  by  Hartraa  and  axtandad  by  Fock  to  include  exchange  effects  (for 
a  recant  ravioir,  sea  Hartraa These  ara  obtained  by  nuMnrical 


M 


solution  of  s  systsm  of  sinultaneous ,  inhomogeneous,  non>lineer,  integro 
differentiel  equetions,  the  so-celled  Fock  equetions.  Efforts  heve  elso 
been  isede  to  include  the  effects  due  to  configuretion  mixing  end  correle- 
tion  of  electrons.  (Sec  Hertree^^^^^),  The  procedure  is  very  leborious, 
end  has  to  be  performed  seperetely  for  each  configuretion  of  each  atom. 
Since  numerical  methods  are  used  the  resultc  are  given  in  a  tabular  form, 
and  cannot  be  used  for  further  analytical  calculations.  It  is,  therefore, 
not  possible  to  develop  analytical  methods  of  great  generality,  but  every 
problem  must  be  solved  separately.  Thus  although  the  theory  was  developed 
in  all  essential  details  over  30  years  ago,  a  list  of  atoms  and  their  con¬ 
figurations  for  which  the  self-consistent  field  radial  wave  functions  arc 
available  is  pitifully  small.  In  particular  it  is  clear  that  if  one  wants 
to  calculate  atomic  parameters,  such  as  transition  probabilities,  for  a 
large  number  of  atoms  one  must  develop  methods  of  considerable  generality 
which  can,  at  one  time,  be  applied  to  a  large  nua^er  of  eases. 

On  account  of  their  obvious  usefulness,  many  attempts  have  been 
made  to  obtain  analytic  wave  functions  (for  example  Zener' 

Eckert Slater^®”^^^ ;  Mores,  Young  end  Haurwlts^®**'*^  Tubls^®"^®^; 
Norse  and  Yilsms^*"^^^;  Breens^*"*'®  Roothan,  Sachs  and 

Weiss^®*^^).  Soam  of  these  wore  merely  ettaapts  to  obtain  analytical 
expressions  of  which  the  parameters  were  deteneined  by  comparison  with 
the  tabular  wave  functions  of  the  self-consistent  field  method. 

Since  the  advent  of  qusnttaa  theory,  it  has  been  very  tempting  to 
use  the  exact  results  obtained  for  hydrogen,  for  the  more  complex  atoms 
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also,  by  th«  •ubatltutlon  of  olthor  on  "offactivo"  principal  quantum 

mimbar  n*  for  n,  or  a  "aeraanad"  nuelaar  eharga  (■  Z  •  s^)  for  tha 

O  Q» 

actual  miclaar  cbarga  Z.  Kara  a^  danotaa  tha  acraaning  of  M  alactron  CH 
by  all  tha  alacttama  ki  tha  atom.  Tha  aeraanteg  paramat«ra  hava 
baan  datarminad  usually  by  a  aaad'ampirical  nathod  (Pauling^**^^^ ; 

Pauling  and  Shaman^^^^^),  although  a  variational  traatnant  has  baan 
appliad  to  aoma  sinpla  casas.  Tha  intarast  in  tha  usa  of  hydrogonic  i»va 
functions  with  a  acraanad  nuelaar  charge  has  baan  ravivad  through  a  nathod 
for  calculating  tha  acraaning  paramatara  davalopad  by  Lay  tar 

1.3  THE  HYIAOGEMIC  HAVE  FUNCTIOMS 


Tha  acraanad  hydroganie  radial  wava  functions,  briefly  rafarrad  to 
aa  tha  "hydroganie  wava  fimetiont",  ara  tha  aolutiona  of  Equation  1-9  with 
V(r)  raplacad  tha  Goulamh  potantial  dua  to  a  acraanad  nuelaar  eharga 

E^  (■  Z  •  share  a^j  is  a  acraaning  paraaatar  and  Z  la  tha  atoadc 

nuakar.  Thay  ara  given  by 


PC«t  7nl  jrjs 


H 

<»-10) 


nhara  n*  danotaa  (n-i-l),  tha  nunbar  of  nodes  of  tha  fimct^  Md  |P|  is 
a  eonfluant  Iqrpargaanatric  kmation  dofinod  below,  lacauaa  of  tha  w^ 
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Equation  B-10,  wa  ahall  find  it  eonvmlant  to  aubatltuta  In 
tha  abova  aquation. 


=  iZ-4s^. 


n 


(B-ll) 


nhora  Y  *<10010  2Z/n  and  (■  Zs^/n)  la  a  raduead  aeramli^  poranatar. 
Wa  BOW  hava 


(B-12) 


whara 


I/a 

=  rja+oi—l 

Hi*  [(a *>)«;!  J 


and 


w  >  -t  ; 


(B-13) 


lF|  (  ('♦•a  ^ 


(B-14) 


idiora 


'Inf  -  (  ')^  ' 


(B-15) 


91 


It  li  to  b«  noted  thet  different  for  different  ni  end  that 

the  redial  wave  functions  (Equations  B*10  or  B-12)  for  the  sasM  n  but  dlf~ 
ferent  I  are  not  strictly  orthogonal  a  It  Is  possible  to  sake  corrections 
due  to  the  lack  of  orthogonality  of  the  wave  functions,  but  they  are  not 
considered  in  this  report. 


For  subsequent  use  In  calculating  the  Slater  F  integrals  we  write 
below  an  expression  for  the  square  of  the  radial  wave  function  (Equations 
B"12  and  B^IA). 


#• 

an 


(B-16> 


where 


i  f-i 


( B*  n » 


B.4  RESUME  OF  LAYZER'S  SCREENING  THEMIY 
(B''23) 

Layser  has  shown  that  the  eigen  values  of  energy  will  possess 

the  correct  quadratic  and  linear  Z«dependence  If  the  one-electron  radial 
wave  functions  setlsfy  the  condition  thet, 

P  (nlyjr)=  p''  ("ly.O+OMyiO 
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where  the  super script  H  denotes  e  hydrogenic  wove  function  end  0(n/  ^ :r ) 
denotes  s  term  of  order  one  In  end 

0  if) 

\r) 

so  thet,  et  leest,  for  lerge  "y ,  0(ni  y  ix)  term  cen  be  neglected.  Note 
thet  our  notetion  differe  slightly  from  thet  used  by  Lsyser. 


Neglecting  0(ni  y  it)  terms,  end  epplylng  the  verletlonel  principle 
Leyser  hes  Shown  thet  the  screening  peremeter  n^,  which  gives  the  screen^ 
Ing  of  the  electron  a  by  ell  the  other  electrons  In  the  etom,  Is  given  by 


•n-.- 


4  ~^V. 

^  I*  iv., 


(B-20) 


where  denotes  the  metrix  elements  of  the  CoulosA  Interectlon  energy 
end  ^  denotes  the  number  of  electron  in  the  subshell  a.  He  cen  use 
severel  different  epproKlsmtlons  for  the  Coulomb  Interectlon,  for  exasple 
we  msy  or  mny  not  Include  ccmflguretlon  odLxlng  end  we  may  either  consider 
the  everege  energy  of  e  given  conflguretlon  or  we  msy  edd  to  It  the  con¬ 
tributions  which  come  from  different  ^ectroscoplc  terms  (cheracterlsed 
by  quantum  numbers  S  and  L),  balonglng  to  the  same  configuration. 


In  the  work  reported  here,  we  have  not  considered  conflguretlon 
mixing,  but  we  consider  the  energy  at  Individual  tMms  as  woll  as  the 
averi^e  energy  of  the  conflguretlon.  When  required,  the  screening  con¬ 
stants  which  are  dependent  on  SL  cen  be  labelled  explicitly  as 
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whereas  those  based  on  the  average  energy  will  be  labelled  as 

The  term  energies  as  functions  of  the  screening  paresieters  ere 

given  by 


w(t) 

where  the  susnetion  ia  over  ell  the  subshells  in  the  atom.  This  is  a  very 
simple  fonaule  for  calculating  the  term  energies  from  the  screening  perem^ 
eters. 


B.5  THE  COUL(»ffi  IMTERACTlCHiS  OF  ELECTRONS 


The  most  conseonly  uaed  epproKimotion  for  complex  etosis  (see,  for 
example,  Condon  end  Shortley^  is  the  sopereble  wove  function  epprox> 
imetion  in  which  each  electron  is  eeauMid  to  be  in  a  stationery  state  in 
the  field  of  the  nucleus  end  ell  the  other  electrons.  The  wove  function 
of  an  N>electron  etosM  is  expressed  in  the  form  of  a  determinant,  smde  up 
of  the  wave  functions  of  individual  electrons,  ea  follows. 


VI, (o')  u^ca‘)  . 

. 


(B-22) 


UjCrf)  ^xLa!‘) . 
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I 

I 

I 


where  u  is  the  weve  function  of  the  eingle  electron  problem,  e^  denotes 
the  eet  of  quentum  numbere  (n^  m^)  end  eubecrlpts  of  u,  from  1  to 

N,  denote  the  electron  to  which  thle  weve  function  belongs.  The  deter- 
minantel  form  ineuree  the  entieymmetry  property  of  the  weve  functions. 

For  e  meny-electron  weve  function  of  the  ebove  form,  Slater 
has  shown  that  the  matrix  elements  of  the  mutual  Coulomb  interactions  of 
the  electrons,  *ums  of  the  matrix  alements  of  the  corresponding 

two  electron  problem. 


Ignoring  configuration  mixing,  the  average  electrostatic  energy 
of  a  configuration  may  be  written  as 


(B-23) 


where 


and  the  two  electr^  interaction  (a,^}  is  given  by 

where  the  prisM  over  the  summation  sign  indicates  that  k  does  not  assume 
all  integral  values,  but  for  all  even  int^iral  values  and  for  either 
all  avm  or  all  odd  integral  values,  denotes  the  lesser  of  and  i^. 
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Th«  eoafflelanti  £  <i^  /^)  and  g  (i^  i^)  aria*  fron  the  angular 

and  tha  apin  part  of  th*  wav*  function*.  Thaa*  can  ba  caleulatad  by  the 

wall  known  laatrlx  nathoda.  For  datalla  of  thaaa  calculation*  th*  raadar 

la  referred  to  Condon  and  Shortlay'^  '  '  and  to  the  more  powerful  nathoda, 
(B«2S) 

due  to  Racah'  ,  baaed  on  tensor  operators.  These  coefficients  are 
Hated  In  Tablas  B*1  and  B-2. 

k  k 

F  and  G  are  Integrals  over  the  radial  part  of  the  wave  functions, 
also  known  as  Slater  Integrals.  They  are  defined  by 


k!' (n*  A 

with 

and 

S'  ~  R  y'pV  (B.2g) 
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TABLE  6>1 


COEFFICIENTS  OF  SUTER  INTEGRALS  FOR  EQUIVALENT  ELECTR(»fS 
(For  •  tabulated  value  of  all  non-vaniahlng  coefficients  are  listed) 


f^(f^  tJ 

_ O-flL 


_ SLJL 


0 

1 

2 

3 


1 

1 

1 

1 


0 

0 

0 

-2/25 

0 

0 

-2/63 

-2/63 

0 

-4/195 

-2/143 

-100/5577 

TABLE  B.2 


COEFFICIENTS  OF  SLATER  INTEGRALS  FOR  NON-EQUIVALENT  ELECTRONS 
(Ftor  a  tabulated  coi^ination  (i  ,  i.)  all  non -vanishing  coef¬ 
ficients  are  listed)  ^ 


0  0  I  -1/2 

0  110 
0  2  1  0 

0  3  10 

111  -1/6 
12  10 
13  10 

221  -1/10 

2  3  10 

3  3  1 


0  0 
-1/6  0 
0  -1/10 
0  0 
0  -1/15 

-1/15  0 

0  -3/70 

0  -1/35 

-3/70  0 

0  -2/105 


0  0  0  0 

0  0  0  0 

0  0  0  0 

-1/14  0  0  0 

0  0  0  0 

-3/70  0  0  0 

0  -2/63  0  0 

0  -1/35  0  0 

-2/105  0  -5/231  0 

0  -1/77  0  -50/3003 
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i^tre  and  danota  raspactlvaly  tha  aaiallar  and  the  greater  of  the 
two  radlua  vector  a,  and  r2>  of  the  two  electrona. 

The  electrostatic  energy  for  the  different  tame  belonging  any 
configuration  can  now  be  written  as 


where  for  some  of  the  sin^lar  cases  is  as  follows: 

(i)  All  shells  completely  filled,  or  there  is  a  single  incompletely 
filled  shell  with  either  one  electron  or  one  hole.  In  these  cases  the 
configuration  conaists  of  a  single  term  and, 


0 


(ii)  There  ia  a  single  incompletely  filled  shell  (designated  by 
n^  with  m  electrons 

Note  that  the  screening  paraawtars  are  functions  of  the  partials 
k  k 

of  the  Slater  integrals  F  and  G  with  respect  to  the  appropriate  screened 
nuclear  change  and  F  and  G  ,  which  we  evaluate  using  the  hydrogenic 
wave  function  (Equation  r'-12),  arc  thaswelves  functions  of  the  screened 
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nuclear  charges,  and  therefore  of  the  screening  paraateter  Is 

tharefere  readily  seen  that  the  rlglit«hand  side  of  Equation  2-20  Is  a  non¬ 
linear  function  of  all  the  screening  parameters,  n2 . adien 

N  Is  now  the  number  of  distinct  subshells  of  electrons  characterised  by 
n  and  i  quantum  nuadMrs. 


I 

I 

I 

I 
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B  6  THE  SIATER  INTEGRALS 


From  Equations  (B>26).  (B>27)  and  (B>28)  it  Is  easy  to  see  that 

Ic  k 

F  and  G  are  synaietrlc  in  the  two  electrons.  However,  on  Mcount  of 
the  occurrence  of r^  and  In  this  equation,  It  Is  necessary  to  break 
up  the  Integral  Into  two  parts,  which  for  the  F  Integral  we  write  as 
follows: 


(B-30) 


with 


Jo  Jo 


(B-31) 


and 


oD 


00 


=  rC^A^.r.) 


(a-32) 


f *■  (<K  1 ;  r)  r  T**  [  P  ( y\  e  y ;  r)] 


(B-33) 


Note  that  Fj^*^  and 
It  Is  easy  to  prove  that 


are  not  synnetrlcal  In  the  two  electrons. 
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(B-34) 


For  the  G  Integral,  theae  two  parta  turn  out  to  be  equal  and  we  have 

> 

fVj.V.  _  0  1  1  _____ 


S  (viifiviikjY.VO=^ 


,“  /  "f*-  K 

rr 


O  Jc> 


(B-35) 


.  n )  f, )  PfUifiY 


The  calculatlona  of  the  Slater  F  and  G  Integrals,  using  hydro - 
genic  wave  functions,  is  extresttly  laborious,  for  electron  shells  whose 
wave  functions  possess  more  than  one  or  two  nodes.  The  author  is  not 
aware  of  any  published  work  where  these  have  been  calculated  as  functions 
of  the  screened  nuclear  charges  of  the  two  electrons  Involved,  nor  are 
any  general  formulae  for  calculating  them  available  In  literature.  We 
have,  therefore,  developed  general  formulae  which  are  suitable  for  hand 
computations  as  well  as  computations  using  digital  computers. 


where 


(B.36) 


(s-37) 
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y.  in|42Wi4K 

i 

•ZZp  lf(h.t«A)(^) 

C"' 

fT  i 

'ff  («,{. . « Jo = • 


(B-38) 


(B-39) 


(B-40) 


(B-41) 


(lU><i+i-K+t+\)!  t  in, + J»>J-\  -vH-i-'t)'. 

'  tM>'.  -t '.  (,iyi,+  jr)t-\+li-i-»'>i 

«nd  o^j  Is  defined  by  Equation  (^17).  A  synbol  such  as  denotes 

Che  lesser  of  the  two  quantities  inside  Che  bracket. 

_W  -  .  ^  X  a  a. 

G^Cn,U^^',Y^>)=9,ft^  I  (M2, 

X _ 1  f  CW>>it«o(^lw 


uhere 


(B-43) 
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(B-44) 


V- 


‘T-  (B-45) 

D  *  -I,  \  (B-46) 


For  tha  cace  of  t%io  oqulvslont  cloc  Croat  and 
wt  have  Cha  following  tliq)la  axpraaalont, 


are  Identical; 


(B-47) 


where  t  and  D  are  now  given  by 


(B-48) 


and 

0  ,  al  +  \ 


(B-49) 


and  4  (t,D)  la  given  by  Equation  (B-44  ). 

A  nuaber  of  foroulae  uaing  the  above  equations  are  given  in 
Table  B~3.  Machine  conputationa  for  Cp 

coeplated  for  all  two  electron  coabinationa  cheek  aarked  in  Table  B-4« 
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TABLE  B-3 


SLATER  INTBCatALS 

(Th«  formulae  listed  here  are  for  F*  (op;  *nd  G*  (op;  X^l),  defined 

by  Equation  (b-50),  where  subscript  op  are  omitted  for 

brevity. ) 

F‘’6vs)=-^['-  0+XJ^(hY%  3X-tO] 

F‘0s,;lP)=  ^[\-  Ci+xf(3)'-^\)] 

fY3s,>4;  =  -^  [  1  -  ('+x)''(>+i2  5^49?  138  xIh^x+'o^+I 

-j^ [l  -  (i4lOCto)t%‘+5V4(&X+i)] 

-  6+w'’ (350/ r318X'^4XS-a  f 
+  Io8X%17X+3;] 
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TABLE  B-3  (continutd) 


I 

i 

I 

I 

I 

I 


F  (3p7  3Jj=r  -  (li-X)  ^\o  X 

-4-  loX  ■^\)J 


-J  r 


k  <  .  'I/V' 


4'53oX4Uoyi-ail^  +5) 


Cr" )t  0+W^(s  x^~  >s  y-t (%;> 

G0s,3s)=  jX()+xJ  (5x-3.)dVwx-i2.x:4ai; 
G'6sAe=  m  y^(i+y; 

&6s>)»>  y  y^(\+yy^(5-fcX*_  ?HX-»  l^sJ 
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TABLE  B-3  (continued) 


{ 

I  GrCi‘Afe)=  +  O'**  ^ 

*  ^  f  0-^X)  t4-JoX%>5-9Sy-483)^4t^ 

I 

2  i 

I  G(is^3dj  =  IS  X 

I  G'Cis.vlj=  a  )L^(HX~)\m')[-i.}>x+^) 

G'(3s,}<|>  i. x^(n- xj%'^vj^-^5'9s  fMV\ri''-ynX^iT) 

A  '^^i.  X-7oX*42?-6X'’-3ioXl  1 87) 

G°(ir;ip>  3X^(itx3"(i««>x-  "ttt  X 

^xV>+X)"(isx"'--4-s-  y+3ff) 
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TABLE  B>3  (continued) 


-u 


\i(,X 


\  u  ^  V 

t  O-^y)  0'^^X-Ui^XiiS^) 
G\ilf>,^A>  X^(i+r)  (2^o5■y’l  21  «.X-^^.ltJ) 
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TABLE  B-4 

THE  TWO  EUCTKM  COMBIHATIONS  FOR  WHICH  SLATER  INTEGRALS  HAVE  BEEN  CALCULATED 
(MaTked  by  X.  Becauac  of  syowtry,  only  the  combinations  above  the  diagonal  are  so  marked.) 
l»  2s  3s  4s  5s  6s  7s  8s  9s|2p  3p  4p  5p  6p  7p  8p  9p|3d  4d  5d  6d  7d  8d  9dl4f  5f  6f jSt. 


The  set  of  simultaneous  equations  (b-20)  contain  partial  derl- 

Ic  Ic 

vatlves  of  F  (a,p)  and  G  (a,^)  with  respect  to^  as  well 
k  k 

though  F  (a,^)  and  G  (oc,^)  are  synmetrlc  In  the  two  electrons  Of  and 

their  partlals  are  not  symmetric .  From  the  preceedlng  expressions 
k  k 

for  F  and  G  the  following  general  properties  of  these  tiK>  Integrals 
can  be  derived  very  simply, 

.Tfa  '  V) 


(B-50) 

(B-51) 


where  R  stands  for  either  Integral  F  or  G  and 


(5-52) 


The  following  expression  for  the  partial  derivatives  are  also  derived 


easily, 

T>Vk 

V,0 


2  \  X^) 


(&*53) 


(r-54) 
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It  Is  therefore  clear  that  both  partial  derivatives  of  F  and 
G  Integrals  can  be  written  as  functions  of  a  single  parameter 

The  case  for  tvo  equivalent  electrons  Is  a  particularly  simple 
one;  the  partial  of  F  (a, a)  la  merely  a  nusdier,  i^lch  isi  In  feet,  the 

value  of  this  Integral  ior'Y  =  ,  the  case  of  hydrogen  atom  (with 

VO.  n^ 

atomic  number  unity)  . 

=4  '  THE  SCREENING  PARAMETERS 


We  now  rewrite  Equation  (li-20)  explicitly  In  terms  of  the 
partial  derivatives  of  the  Slater  Integrals 


where 


Is  merely  a  number,  Independent  of  any  of  the 

and  . 


(3*57) 


(B-58) 


(B-59) 


f^C'i"  SL)  and  g^C'/"  SL)  denote  the  coefficients  In  which  the  SL  depen¬ 
dent  contributions  to  the  electrostatic  energy  are  explicitly  included. 

The  sysdx)!  denotes  the  orbital  quantum  nuad»ers  of  all  the  incompletely 
filled  stalls  Including  the  parentages,  etc. 
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We  note  thet. 


and  for  the  caee  Z  -»  oo 

Y  .  •y“ 

/•<?-*  ^<<^  =  "R7  ('-«) 

The  nethod  for  confuting  the  ecreenlng  paranetere  ie  ae  followa. 

"b 

We  flrat  calculate  a  aet  of  n'a  ulth  X  .  =  .  We  call  theae  aaya|>totic 

n'a  or  .  With  theae  'a,  and  none  auitably  chooaen  large  value  of 
Z,  we  compute  a  new  aet  of  X  ,  and  uae  theae  to  recalculate  n  'a  by 

Qp  Cl 

Iteration.  We  then  decrease  Z  by  aucceaslve  atepa  and  keep  repeating  the 
same  Iterative  procedure.  Thus  we  work  along  an  lao-electronlc  sequence 
and  obtain  the  screening  paraamters  for  any  desired  number  of  atoms  and 
Ions  belonging  to  this  Iso •electronic  sequence.  Since  the  procedure  Is 
very  time  consuming,  a  program  for  computations  on  a  digital  computor 
has  been  developed* 


The  asymptotic  screening  constants  for  two  electron  cooiblna- 
tlons,  and  using  the  average  energy  of  a  configuration  only,  can  however 
be  computed  quite  easily.  They  are  defined  as  follows: 


(S-62) 


and 


for  a  4  fi 
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<B-63) 


The  asymptotic  screening  constants  for  the  average  energy  of  a 
configuration  can  then  be  constructed  from  these  two  electron  screening 
constants  as  follows: 


(B-64) 


b.8  the  calculation  of  the  transition  integrals 


The  only  case  of  interest  for  an  electric  dipole  transition  is 
when  the  azimuthal  quantum  number  of  the  active  electron  changes  by  one 
unit.  Therefore,  the  calculation  of  the  following  integral  is  of  suf¬ 
ficient  generality: 


)tii  -{-I 

Jo  (B-65) 

=  9.9;.Y.  K 


where 


.0* 


ir-'-s 


(B-66) 


A  sisdlar  integral  with  been  evaluated  by  Gordoii^^^^\ 

The  necessary  forsulas  for  the  evaluation  of  Integrals  involving 
confluent  hypergeonetric  functions  have  been  given,  assmg  other,  by 
Landau  and  Lifshits^^^^^  in  their  mathesHtieal  anpeadix. 
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The  p  index  een  be  reduced  by  using  the  foroule 


(B-67) 


wherees  the  e  Index  cen  be  reduced  by  using  another  foroule  . 


I 

I 

1 

i 


Successive  applications  of  Equations  (6-67)  and  (b-68)  to  (b-*66)  yield 

3ji(+l  ,-Vll  )-  -J.a  M|  J  -Yr 

.  rtV  y>f)-  5  [Yi  "r+'>  Ya  (♦+ 


J,.  can  be  evaluated  in  closed  fora  by  another  forsula 
OlO  -da  Y'  j  ®^|4di-Yr 

3;r  («‘..*v),c-»  aRO(V,+Y)  (Y»->0 

^Y.Yv/CYv-lrO*) 

where  2^^  «  hypergeoaetric  function  given,  In  general, 

iF.C«Ac,f9=  It  51. r- 

^  o^Cq’»0(a4a)b(HiX^a)  3^ . 

C(C+|)(C4J9  3\  I 


(B-70) 


(B-71) 
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The  series  termlnetes  after  a  finite  nus^r  of  terms  irtien  a  and/or  b  are 
negative  integers.  On  substitution  of  (B-70)  in  (B-69)  and  simplifica¬ 


tion  «e  get  p 

A  -M  OH'.  1  <ai 


where 


I  V  »  - 


(B-73) 


The  following  relation  nay  be  used  to  reduce  the  nuad)er  of  hypergeometric 
functions  from  three  to  two  in  (»*72)9if  desired. 

(c  -  la.  -  Cb-  (»■>  -»  ^*■^(8-74) 

where 

Combining  (B-6S)  with  (B-72)  we  get 

(B-75) 


_  a  cti  iF  (-  »>t-  0 

Vx+T,  ^ 
-vTf,)  aP* 
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2  T.  2 

which  for  th«  spocial  camY’i*;!  *Yt.' TT  *  rodueoi  to  tho  following 


fB  *26) 

equation  firit  darivnd  by  Gordon^  ' . 


Ml  MO 


1+1 


rha  following  special  cases  can  be  obtained  easily 

V  t 

vV'i K)  =32  0-77) 

C.  (.v»f.J'iV,VO  *32(w-0^''5f  Vi’’  ["Vr  ■>»(»»■**•) 

T.(af,x‘  iYiYO  T  If-"  t'f****'-’^ 


(B-78) 


CB-79) 


rs-go) 


-  i  (>+a« 

CYx-Y.^''^ 

O.-h)'"-'** 
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B.9  ALTERNATE  EXPRESSION  FOR  TRANSITION  INTE(HtALS 


For  th«  8«k«  of  gonorallty  we  shall  evaluate  the  following 


Integral , 


‘  •r’’ 


(B-81) 


For  Coulomb  potential  wave  functlonst  we  shall  later  put  q  equal  to  sero. 
p  equals  unity  for  the  electric  dipole  case,  and  It  equals  2,  3,  4,  etc. 
for  higher  multipole  cases.  Substituting  from  Equation  (b-12)  we  get. 


(B>82) 


This  Integration  can  be  performed  using  certain  standard  results  In  tlM 

(B~28) 

theory  of  confluent  hypergeomstrlc  functions,  (see  for  exaeple  Slater  , 
pegs  54).  isrforniai  the  integration  we  obtain,.  .  I  « 


'  2.^  L  oir-i> hoh- vj  V 


*A  slsdlar  result  has  also  been  obtained  by  Frof.  D.  H.  Msnsel  (wi^ublished) . 
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where  the  expreeelon  In  the  squere  bracket  following  the  euBBetlon  sign 
can  be  readily  aeen  to  be  equal  to  hyper geoae trie  func¬ 

tion  defined  by  Equation  (^-71).  The  conputatlon  of  the  expreailon  on 
the  right  hand  aide  of  Che  above  equation  la  quite  atralght  forward  In 
principle,  but  becoaaa  very  laborloua  aa  n^  and  n^  beeoae  large. 


He  ahall  now  proceed  on  to  the  facial  caae  we  are  Intereated 
In,  naaaly  the  caae  where  p  >  1  and  q  >  0.  Fbr  thla  caae, 

iF,  (-«r,<,+-(’x+'++t  >  V, 

^ C  J  On?)'. 

'^r  i'  ~ 

I'Vf+v.) 

_  2.V  dt  j  »lr«-  'X'lr 

*  («,++,.*».»  O*  Ov+V'--' 

I.WnA'.yvOr 

fv 

«'  0  (B.SS) 

« 

y\  <B-«6) 

’  (V^>3^ 
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to  say  that  tha  aaat  axpraaslons  can  also  ba  obtaiaad  by  intagra* 
tlon  of  Equation  (b~81)  with  tba  radial  iMva  functions  glvan  by  Equation 
(B-12)  coad)lnad  with  Equation (b >14).  Sosib  further  slapllflcatlons  of  tha 
above  expression  can  be  brought  about  as  follows.  The  conson  demndnator 
of  all  tha  tarns  under  tha  two  sunnatlOM  la  ^ 

tha  consion  daoosdnator  and  parfornlng  blnonlal  expansion  of  the  resulting 
tern,  nanely  (k  “2  wa  obtain 

J  ^  ^  .i 


' ip>. 


hi 


Eeplaclng  p  by  another  Index  u,  such  that  u  ■  p  J,  wa  gat, 

^  ^3  ^“*3' 


L+i)i  ( -t  -  j  ) 


Yvn'" 


Tm  ^  Vi  / 


=  T  'y  CM. yhL)p - -5— ^ 

u« 


(B-M) 
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.h.,.  i  i 

4i  Ji.  ^ 


C^“pl 


(B-89) 


u  takes  on  all  integral  values  fr<»  0  to  (n|  +  nj).  Is  a  function  of 
the  quantum  nuad)ers  of  the  two  electrons  under  consideration,  but  it  is 
independent  of  the  screened  nuclear  charges,  and^2*  Hhen  n^  and 
n^  are  so  large  that  desk  calculations  are  too  tedius,  the  calculation 
of  on  electronic  computors  is  a  trivial  matter,  giving  us  analyti* 
cal  expressions  of  in  terms  of  as  follows. 


W|-fVlvi-7_ 


(B-90) 


Y'- 


It  should  be  noted  that  whereas  is  a  syaaccric  function  of  the  two 
electrons,  if*  is  not. 


Table  2-S  gives  expressions  for  for  a  selmetod  number  of 
transitions . 
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TA&E  B-5  DXFOLI  IIRIGRALS 
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I,(Tr.U)  - 


1  “ 

I 

I  I,(3f-^-)  = 

I 

t  I,  (4f.2*)  ■ 

I,(5f,2>.)  = 

= 


TABJ  B-5  (eont.) 


‘‘S') 


7jf 
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I(2f,5«) 

I(3f,5.) 


TAKiE  B-5  (eont.) 
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TABLES *5  (eont.) 


--2e">(»-^>’ t  ^  AOiZj^^ 

-™(a- "'^a-V 

t2i4(/ 1  ''’j‘^)(a'i‘/ 

♦W20(3^*^S  »oj.Y)4.^J.) -78?20/^ 
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TAn,I  B-S  (cone.) 


y  JO  -  ioij.  ’ 

.3Z(7^S 

t384(6J^<jS  Woj.*j» 
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TAm.SB-5  (cont.) 


O’30'* 


-70(3jf^. 
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1 

I 

1 


TAH.B  B  -5  (cont.) 


^('8^ 

- 14(^  ’  y!?f 

+3«.o(s^'ji  t\8^‘jf  ••%/)(5'''>)' 

-Z44.4(30^?j;+/^')(5,.^) 

•t59>5b^?^f| 
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TAK.B  b>5  (eont.) 


I 

i 

i 

I 

! 


.32(3SjJ  1 350^^ 

-9b(p  5as^'^.t6io^jN28<^|»'^)(jf^v)* 
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i  UWf)  - 

I 

I  = 

I 

I  I(6<l,2f)  - 

I 

l(7J,Zf)  = 

um  - 


UM^)  “ 


TAH.E  B>5  (eont.) 
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TA&B  B-5  (eont . ) 


I(6<1.4f)  - 
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TABLE  B-5  (eont.) 
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TAM<BB-5  (eont.) 
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TA&B  B*5  (eont.) 


132 


I 


B*5  (eont.) 


■  {^*^f  ■(V3‘)(>^l'^ 

m 

tg(2’#  +‘3jLjf7v*X}‘*>f 


B .  1 0  Results  of  Conmutat ions'* 


(a)  The  Normalization  Factor,  Q 


Q  .  is  defined  by  Equation  (B*13) 
n£ 

1/2 


L2n(ni!-l):J 


1 _ 

(Yi+1): 


(B-13) 


table  of 


The  following  recursion  relations  are  useful  in  extending  a 

1/2 


2i<2i-M) 


*n,i-l 


'n/ 


■[ 

[n«l 


nl/2 
Eli  iLtili  o 
n  n-i  ^n-1,2 


(B-91) 


(B-92) 


(n+i)(n+i-l)  22(22+1)  ‘^n-1,2-1 


(B-93) 


A  few  special  eases  of  sre  given  below: 


^,0 


•1/2 


(B-94) 


i.l  6  L 


1/2 


^.1  “  6  I  2  J 


^n.n-l  •  (2n)! 


.-1/2 


(B-96). 


A  staple  procedure  in  calculating  is  to  start  with  which 
-1/2 

is  always  2  ,  and  successively  calculate  Q^2>  using  Equation  (B-91.) 

until  one  arrives  at  Equation  (B-96)  can  be  used  as  a  final  check 


*  The  work  perforswd  for  this  section  was  sponsored  by  both  this  contract 
and  Contract  Ro.  AF  33(600) -426S7 . 


The  values  of  for  n  up  to  9  and  all  allowed  values  of  i  were 
computed  in  the  above  manner,  and  subsequently  machine  confutations 
were  also  performed.  Table  B-6  lists  the  Q 

Til 

(b)  The  Coefficients 

The  d|^^*s  are  the  coefficients  occurring  in  the  hypergeosmtric 
function  ^(n-i>l),  2i-f2:x^  .  They  are  given  by  Equation  (B-15). 


“nf  ^  (n-f-l-i)!  (2f+l+i^!  i  ! 

(B-15) 

For  i  >  n-f-1,  d^jj  vaalshat. 

The  following  recursion  relations  are  useful 
the  table  of  d^: 

in  extending 

nf  (21+1+1)1  ®nl 

(B-97) 

A  , _ A 

nJi  (n-lH21+i)(21+l+i)  n.1-1 

(B-98) 

,1  n-1-1  ,i 

“nl*  n-l-l-i  “n-l,l 

(B-99) 

**01  “  (21+l)fle+l+l)  **0-1, 1-1 

(B-lOO) 

A  few  special  cases  of  d^  are  of  interest 

d®,  -  1  (B-101) 


13S 


TABLE  B-6 


-1/2 

«  2  ‘  ^  for  all  n 


n.i 

Q 

n,i 

Q 

2P 

.2041241^ 

8p 

.935414»^ 

3p 

.333333B+00 

8d 

.3622848+00 

3d 

.3726781-01 

8f 

.6397088-01 

4p 

.45643SB-K)0 

88 

.6155598-02 

4d 

.790569B-01 

8h 

.3494708-03 

4£ 

.4980121-02 

81 

.1185408-04 

5p 

.57735QB4-00 

8k 

.2186208-06 

5d 

.132288B-(-00 

9p 

.1054098+01 

5f 

.125988B-01 

9d 

.4624818+00 

ss 

.524950B-03 

9£ 

.9343538-01 

n 

.697216B-K)0 

9g 

.1046258-01 

D 

.197203B-H)0 

9h 

.7117658-03 

.24397SB-01 

91 

.3060688-04 

B 

.1315401-02 

9k 

.8244681-06 

6h 

.456911B-04 

9i 

.1249778-07 

7p 

.816496B+00 

7d 

.2738611^ 

7£ 

.4123931-01 

7t 

.3290301-02 

7h 

.1465371-03 

71 

.3386851-03 
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(B-102) 


2(^+1) 


(B-103) 


A  table  of  for  values  of  n  up  to  5  was  prepared  using  the  above 
recusion  relations  and  subsequently  machine  computations  were  perforond 
for  n  up  to  9.  These  are  given  in  Table  B-7 

(c)  The  Screening  Parameters 

Equations  (B-56)  define  a  system  of  simultaneous  non*linear 

equations  in  which  the  right  hand  side  is  a  function  of  all  the  screening 

parameters  ic^^,  of  a  given  problem.  A  Fortran  program  for  use 

on  digital  computers  has  been  developed  for  solving  these  equations. 

Ue  work  along  an  isoelectronlc  sequence.  Ue  use  as  our  starting  point 

the  case  when  Z  -*  oo .  The  screening  parameters  for  this  case  can  be 

computed  easily  with  given  by  Equation  (B-61).  As  an  intermediate 

step  we  perform  iterative  calculations  for  Z  -  Z  -*40,  where  Z  is  the 

o  o 

lowest  amad)er  of  the  isoelectronlc  sequence,  using  n^'s  for  Z  «  oo  ceee 
as  our  starting  point.  The  iterations  arc  carried  until  the  relative 
difference  between  the  results  of  two  successive  iteration  is  less  then 
a  prescribed  quantity  e  (for  computations  reported  here  €  *  0.0001). 

The  process  is  repeated  successively  for  the  cases  Z  ■  Z^  +  19,  18, 

....Z^.  The  screening  paramsters,  expressed  as  s(a)(*  ^  n^)  are  given 
in  TAbles  B-8  to  Bol9. 
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I 


Jk 


n 


CM 
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TABLE  B-7  (continued) 


I 


Ion 

TABLE  B-8 

8CUBMIMG  PAXAMBTBBS 

Be  I  laoelectronlc  Sequence 

2  2 

Configuration:  la  2s 

Tern:  ^S 

8(ls) 

8(28) 

Ll‘ 

3.2761  B-1 

1.7902  BfO 

Be  I 

3.1861  B-1 

1.8599  B+0 

B  II 

3.1369  E-1 

1.8917  BfO 

C  III 

3.1184  E-1 

1.9025  BfO 

M  IV 

3.1164  B-1 

1.9036  EfK) 

0  V 

3.1221  B-1 

1.9006  Bf-0 

P  VI 

3.1312  E-1 

1.8961  BfO 

He  VII 

3.1414  B-1 

1.8912  B+0 

Na  VIII 

3.1517  E-1 

1.8863  BfO 

1^  IX 

3.1616  B-1 

1.8817  BfO 

A1  X 

3.1708  E-1 

1.8775  BfO 

81  XI 

3.1794  B-1 

1.8736  BfO 

P  XII 

3.1873  B-1 

1.8700  BfO 

8  XIII 

3.1946  E-1 

1.8668  BfO 

Cl  XIV 

3.2013  E-1 

1.8639  BfO 

A  XV 

3.2074  E-1 

1.8612  BfO 

K  XVI 

3.2130  B-1 

1.8587  BfO 

Ca  XVII 

3.2182  B-1 

1.8565  BfO 

8e  XVIII 

3.2230  B-1 

1.8544  BfO 

T1  XIX 

3.2274  B-1 

1.8525  BfO 

Z  ■  00 

3.3353  B-1 

1.8078  BfO 
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TABLE  B-9 

8CBBBNIM6  PABAMBTBIS 
BoT  Iao«l«ctronic  Sequence 

2 

Conflguretion:  Is  2s2p 
Tens: 


Ion 

8(1S) 

»(2s) 

•  (2p) 

Li’ 

3.1832 

B-1 

1.5178 

BfO 

3.0030 

BfO 

Be  I 

3.2321 

B>1 

2.1229 

BfO 

2.0640 

BfO 

B  II 

3.2444 

B-1 

2.1965 

BfO 

1.9749 

BfO 

C  III 

3.2590 

B-1 

2.2165 

BfO 

1.9390 

BfO 

N  IV 

3.2766 

B-1 

2.2186 

BfO 

1.9180 

BfO 

0  V 

3.2949 

E-1 

2.2138 

BfO 

1.9037 

BfO 

P  VI 

3.3126 

B-1 

2.2065 

BfO 

1.8932 

BfO 

Me  VII 

3.3289 

B-1 

2.1987 

BfO 

1.8849 

BfO 

Me  VIII 

3.3437 

E-1 

2.1911 

BfO 

1.8783 

BfO 

Ms  IX 

3.3570 

B-1 

2.1840 

BfO 

1.8728 

BfO 

A1  X 

3.3690 

B-1 

2.1775 

BfO 

1.8682 

BfO 

81  XI 

3.3797 

B-1 

2.1715 

BfO 

1.8642 

BfO 

P  XII 

3.3893 

B-1 

2.1662 

BfO 

1.8608 

BfO 

8  XIII 

3.3980 

B-1 

2.1614 

BfO 

1.8578 

BfO 

Cl  XIV 

3.4059 

B-1 

2.1570 

BfO 

1.8551 

BfO 

A  XV 

3.4130 

B-1 

2.1530 

BfO 

1.8528 

BfO 

K  XVI 

3.4195 

B-l 

2.1494 

BfO 

1.8507 

BfO 

Ce  XVII 

3.4254 

B-1 

2.1441 

BfO 

1.8488 

BfO 

8c  8VIII 

3.4308 

B-1 

2.1430 

BfO 

1.8471 

BfO 

T1  XIX 

3.4356 

B-l 

2.1402 

BfO 

1.8456 

BfO 

X  ■  SB 

3.5M7 

B-1 

2.0773 

■fO 

1.8129 

BfO 
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TAILE  B-10 
SCRBBMlllG  PAIUMBTBBS 
B«  X  Iso«l«etronlc  Scqucnca 

2 

Configuration:  la  2a  3p 
Tana: 


Ion 

•  (!•) 

a(2p) 

a(3p) 

Ll" 

3.1396  E-1 

1.9799  BfO 

2.5502  BfO 

Ba  1 

3.1969  B-1 

1.8484  BH) 

2.7110  BfO 

B  II 

3.2529  B>1 

1.7602  BfO 

2.7775  BfO 

C  III 

3.2957  B-1 

1.7074  BfO 

2.7953  BfO 

N  IV 

3.3280  B-1 

1.6774  BfO 

2.7864  BfO 

0  V 

3.3527  B-1 

1.6598  BfO 

2.7677  BfO 

P  VI 

3.3720  B-1 

1.6490  BfO 

2.7471  BfO 

Ma  VII 

3.3876  B-1 

1.6419  BfO 

2.7276  BfO 

Ha  VIII 

3.4003  B-1 

1.6370  BfO 

2.7102  BfO 

Mg  IX 

3.4110  B-1 

1.6334  BfO 

2.6950  BfO 

A1  X 

3.4200  B-1 

1.6307  BfO 

2.6818  BfO 

81  XI 

3.4276  B-1 

1.6285  BfO 

2.6703  BfO 

P  XIX 

3.4343  B-1 

1.6268  BfO 

2.6603  BfO 

8  XIXI 

3.4401  B-1 

1.6253  BfO 

2.6514  BfO 

Cl  XXV 

3.4453  B-1 

1.6240  BfO 

2.6437  BfO 

A  XV 

3.4498  B-1 

1.6230  BfO 

2.6368  BfO 

K  XVI 

3.4539  B-1 

1.6221  BfO 

2.6306  BfO 

Ca  XVII 

3.4575  B-1 

1.6212  BfO 

2.6251  BfO 

8c  XVIII 

3.4609  B-1 

1.6205  BfO 

2.6201  BfO 

T1  XIX 

3.4639  B-1 

1.6199  BfO 

2.6156  BfO 

X  ■  00 

3.5261  B-1 

1.6072  BfO 

2.5266  BfO 
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VMM  B-11 

SCtmXVG  POAMETBS 
0  I  lio«l«etroiiie  S«qu«iic« 
2  2  4 

Configuration:  Is  2s  2p 
tern:  ^8 


Ion 

■  (!•) 

■  (2s) 

■  (2p) 

h" 

3.3288  B>1 

2.9181  BfO 

3.6713  BH) 

0  X 

3. 4078  1-1 

3.0238  8f0 

3.5794  BH) 

F  II 

3.4825  B-1 

3.0779  BH) 

3.5248  BH) 

Ms  III 

3.5524  B-1 

3.1080  BH) 

3.4873  BfO 

Hs  IV 

3.6170  B-1 

3.1253  BH) 

3.4595  BH) 

Hg  V 

3.6763  B-1 

3.1354  BH) 

3.4378  BH) 

A1  VI 

3.7303  B-1 

3.1411  BH) 

3.4204  BH) 

Si  VII 

3.7794  B-1 

3.1442  BH) 

3.4061  BH) 

P  VIII 

3.8241  B-1 

3.1455  BH) 

3.3940  BH) 

8  DC 

3.8648  B-1 

3.1458  BH) 

3.3837  BH) 

Cl  X 

3.9019  B-1 

3.1453  BH) 

3.3748  BH) 

A  XI 

3.9358  B-1 

3.1445  BH) 

3.3671  BH) 

K  XII 

3.9668  B-1 

3.1433  BH) 

3.3602  BH) 

Cn  XIII 

3.9953  B-1 

3.1420  BH) 

3.3541  BH) 

8e  XIV 

4.0216  B-1 

3.1406  BH) 

3.3487  BH) 

T1  XV 

4.0458  B-1 

3.1392  BH) 

3.3438  BH) 

V  XVI 

4.0682  B-1 

3.1377  BH) 

3.3394  BH) 

Cr  XVII 

4.0890  B-1 

3.1362  BH) 

3.3353  BH) 

Mn  XVIII 

4.1083  B-1 

3.1348  BH) 

3.3316  BH) 

P«  XDC 

4.1263  B-l 

3.1334  IH) 

3.3283  BH) 

Z  ■  00 

4.6176  B-1 

3.0813  BM) 

3.2504  BH) 
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6.1«54 


The  hydrogenlc  wave  functions  using  the  screening  parameters 


are  compared  with  the  self ‘consistent 
oxygen  (Hartree,  Hartree  and  Swlrles^ 
(Roothaan,  Sachs,  and  Weiss In 
genic  wave  functions  are  given  below: 

0  I  (ls^2a^2p^  ^S)  Is  electron: 

0  I  (ls^2s^2p^  ^S)  2s  electron: 

0  I  (ls^2s^2p^  ^S)  2p  electron: 

LI  (ls^2s^)  2s  electron: 

The  agreesmnt  between  the 
remarkable  for  0  1  and  Is  good  for  LI 


fleldwave  functions  for  neutral 
‘  and  negative  Ion  of  lithium 

Figures  B-1  to  B-4.  The  hydro- 

42.394  r  exp  (-7.6592  r) 

7.8491  r  exp  (-2.4881r)  (l-2.4881r) 
8.3864  r^  exp  (-2.2103  r) 

0.94079  r  exp  (-0.66S24r)  (l-0.66524r) 

two  sets  of  wave  functions  Is 


Total  Ionization  energies  based  on  hydrogenlc  wave  functions 

are  compared  with  experimental  values  In  Table  B-20.  The  agreement  Is 

excellent;  for  Be  I  It  Is  as  good  as  that  for  the  multlconflguratlon 

(B*32) 

Hartree-Fock  wave  functions  of  Klbartas,  Kavetskls  and  lutsls; 
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Figure  B-1.  01  (l«^2»^2p^  ^S);  V«v*  Function  of  It  Bloctrma 
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O/C040-4SO 


Figure  B-3. 


01  (l«^2*^2p^  ^g); 


Have  Function  of  2p  iloctron 
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15« 


Figure  B-4.  Ll"  (le^2e^) :  Wave  Function  of  2s  Electron 


TABLE  B>20 


TOTAL  IONIZATION  ENERGIES  (Atomic  Unit) 


Ifydrogenlc  H-F 

Wave  Function  Nave  Function  Experimental 


0  I  (^S) 

-  74.62 

-  74.95 

F  II  (^S) 

-  98.61 

>  98.92 

Be  I  (^8) 

-  14.70  -  14.64* 

-  14.67 

B  II  <^S) 

-  24.38 

-  24.35 

C  III  (^8) 

-  36.55 

•  36.55 

N  IV  (^8) 

-  51.23 

•  51.25 

0  V  (^8) 

•  68.40 

-  68.45 

*BaMd  upon  multlconf juration  Hartree-Fock  calculatlona  of 
Kavatskle  and  luteleJ®"^*' 

Kibartaa* 

(d)  The  Dipol«  Integrals 


The  electric  dipole  integral  is  given  Equation  (B*90) 
in  terme  of  a  coe££ieientY'**(nj^f defined  by  Equation  (B»89). 
Coaputations  on  an  IBM  7090  have  been  perfornsd  for  the  coefficients 
T”  for  all  t«o<*electron  coe^inatione  of  practical  interaet*  ai^  the 
results  are  given  in  Table  B-21  to  B-26.  In  order  to  reduce  huoan 
errors  in  the  reproduction  of  a  large  collection  of  nunbers,  the  print* 
out  fron  the  digital  cooputor  has  been  photocopied.  P8I  7*  for  exesvle* 
means'f^^.  Nuadyers  such  as  0.23999999  or  0.24000011  should  be  rounded 
off  to  0.24. 

2 

A  few  sasiple  values  of  the  transition  integral  cr  cosiputed 

with  the  screened  hydrogenic  wave  function  are  given  in  Table  B*27» 

(B*29) 

and  compared  with  calculations  based  on  Coulosd)  approxisNition'  ' 
and  H-F  wave  functions.  The  agreesHnt  is  satisfactory. 
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(n.p.n.s)  (continued) 
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TABLE  B-27 


TRANSITION  INTEGRAL 


Atom 

2 

Transition  Integral  cr 

Ifydrogenlc 

Coulomb 
Approx . 

Har tree -Foe k 

Be  I 

2.27 

2.03 

1.86 

B  II 

0.977 

1.010 

C  III 

0.544 

0.573 

N  IV 

0.348 

0.368 

P  XII 

0.0522 

0.0540 
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APPENDIX  C 
D.  Layzer 


NEW  GENERALIZATION  OF  THE  SCREENING  THEORY 
AND  THE  HARTREE-FOCK  THEORY 

C.l.  Introduction 


The  distinguishing  feature  of  the  Z-dependent  theory  of 

(C- 1 ) 

atomic  spectra  (Layzer,  )  henceforth  referred  to  as  A) 
is  the  fact  that  it  allows  for  certain  effects  of  configu¬ 
ration  mixing  in  the  first  approximation.  This  ensures  that 
the  predicted  term  energies  have  the  correct  linear  depen¬ 
dence  on  Z.  Experiment  shows  that  the  structure  of  term 
spectra,  in  both  simple  and  complex  systems,,  depends  strong¬ 
ly  on  the  linear  Z-dependence  of  the  energies.  A  knowledge 
of  the  energy  coefficients  and  in  the  expansions 

2  -1 

W  =  W^Z  +  W^z  +  +  0(Z  )  (C-l) 

enables  one  to  describe  quantitatively  the  most  conspicuous 
features  of  term  spectre)  in  isoelectronic  sequences  (see 
especially  Figures  1  and  2  in  reference  C-l). 

In  order  to  calculate  the  coefficients  (W^  is  given 

by  the  well-known  hydrogenic  formula)  one  needs  to  know 

nothing  about  the  radial  dependence  of  the  atomic  charge 

distribution  beyond  what  is  provided  by  the  theory  of  the 

hydrogen  atom.  That  is,  all  the  necessary  calculations 

utilize  hydrogenic  wavefunctions .  One  can,  of  course,  use 

less  simple  wavefunctions  in  these  calculations,  as  was 

(C-2) 

done  in  the  first  calculations  of  this  type  (Layzei  ' 
but  such  refinements  only  improve  the  accuracy  of  higher 
terms  in  the  expansion. 


For  some  purposes,  however,  one  does  require  a  reason¬ 
ably  accurate  knowledge  of  the  radial  charge  distribution 
and  hence  of  the  radial  wavefunctions .  One  needs  this  in¬ 
formation,  for  example,  in  order  to  calculate  electric- 
dipole  transition  probabilities. 

The  various  procedures  in  common  use  for  calculating 
accurate  radial  wavefunctions  are  all  based  on  the  varia¬ 
tion  principle.  They  differ  through  the  restrictions  they 
impose  on  the  form  of  the  variational  functions.  Hartree- 
Fock  wavefunctions  (i.e.,  exact  solutions  of  the  Flartree- 
Fock  integro-dif ferential  equations)  fall  at  one  end  of 
the  gamut  of  variational  wavefunctions,  screened  hydro- 
genic  wavefunctions  at  the  other.  The  latter  involve  only 
a  single  adjustable  parameter,  while  the  former  are  re¬ 
stricted  beforehand  only  by  the  usual  conditions  of  inte- 
grability  and  smoothness.  Other  types  of  variational  wave- 
functions  seek  to  achieve  some  of  the  accuracy  of  Hartree- 
Fock  wavefunctions  without  sacrificing  all  the  simplicity 
of  screened  hydrogenic  wavefunctions.  Whicli  type  of  varia¬ 
tional  wavefunction  one  chooses  to  employ  will  be  governed 
by  the  purpose  in  view.  The  Hartree-Fock  description  is 
certainly  the  most  accurate;  from  a  qualitative  standpoint, 
however,  the  screening  description  offers  unrivalled  ad¬ 
vantages  of  simplicity  and  flexibility. 

Any  meaningful  comparison  between  different  variational 
wavefunctions  presupposes  a  common  approximation  scheme. 
Thus  we  may  compare  wavefunctions  calculated  with  neglect 
of  configuration  mixing,  or  with  certain  specified  kinds 


i;a 


of  configuration  mixing  taken  into  account.  (In  practice « 
the  calculations  that  allow  most  fully  for  configuration 
mixing  are  those  based  on  the  simplest  wavefunctions,  so 
that  strict  comparisons  are  difficult.) 

According  to  the  Z~dependent  theory,  the  moat  is^rtant 
kind  of  configuration  mixing  is  that  between  configurations 
belonging  to  the  same  complex  (i.e.,  the  same  set  of  princi- 
pal  quantum  numbers  and  the  same  parity) .  On  the  other  hand, 
the  theory  of  radial  wavefunctions  has  usually  been  devel¬ 
oped  within  the  framework  of  the  central-field  approximation, 
though  attempts  have  occasionally  been  made  to  allow  for 

effects  of  configuration  mixing.  A  notable  example  is  the 

(C-  3) 

work  of  Hartree^  Hartree>  and  Swirles'  \  which  dealt 

2  ^  n 

with  the  mixing  of  configurations  Is  2s^2p  with  configure- 
2  n^"  2 

tions  Is  2p  .  However,  such  attempts  have  lacked  an  ade¬ 
quate  theoretical  framework,  comparable  in  simplicity  and 
generality  with  that  provided  by  the  central-field  approxi¬ 
mation.  The  present  paper  seeks  to  supply  such  a  framework. 

C.2.  The  Theoretical  Framework 

Let  Y  denote  an  antisymmetric  N-electron  wavefunction 
pertaining  to  a  definite  state  in  a  definite  electronic  con¬ 
figuration.  The  quantity 

W  =  <THY>  ,  (C-2) 

where  H  is  the  N-electron  Hamiltonian  and  the  brackets  <C  y 
indicate  an  integration  over  all  3N  coordinates  and  a  sum¬ 
mation  over  all  the  spin  variables,  is  the  variational  ener- 
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gy.  It  is  a  functional  of  the  radial  wavef unctions .  If 
these  are  specified  beforehand,  except  for  the  values  of 
certain  parameters,  then  W  reduces  to  a  function  of  these 
parameters,  W  exceeds  the  true  energy  of  the  state  under 
consideration  provided  that  'if  is  normalized  and  is  orthogo¬ 
nal  to  the  true  eigenfunctions  of  all  states  lower  than  the 
one  under  consideration.  The  condition  6W  «  0  leads  to  a 
set  of  simultaneous  integro-differential  equations  for  the 
radial  wavef unctions,  or  to  a  set  of  simultaneous  algebraic 
equations  for  the  adjustable  parameters. 

One  can  easily  generalize  the  theory  just  sketched  so 
as  to  include  the  effects  of  mixing  among  a  finite  number 
of  configurations.  Let  'i'^^  denote  variational 

wavefunctions  for  states  with  fixed  values  of  the  quantum 
numbers  SLp  in  different  configurations.  Suppose,  for  the 
sake  of  definiteness,  that  corresponds  to  the  deepest  of 
the  states  under  consideration.  Then  the  corresponding 
variational  energy  W'*'  is  the  lowest  eigenvalue  of  the  matrix 
II  II  .  As  before,  the  condition  6W^  »  0  yields  a 

set  of  simultaneous  equations  for  the  unknown  functions  or 
parameters.  (Note  that  all  the  radial  wavefunctions  or 
parameters  that  figure  in  the  complete  set  may 

occur  in  a  single  equation.) 

An  accurate  solution  of  the  secular  equation  for 
would  lead  to  a  set  of  variational  equations  of  prohibitive 
complexity.  The  procedure  usually  adopted  consists  in 
retaining  only  first-  and  second-order  contributions  to 
in  the  sense  of  ordinary  perturbation  theory.  Thus  is 
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given  by 


+  )  -E-^r —  ' 

i/i  ^  ^ 

where  denotes  the  mean  energy  of  the  ith  configuration. 
Even  with  this  siflqplificatlon*  the  variational  equations  with 
configuration  mixing  taken  into  account  present  considerable 
more  difficulty  than  the  usual  equations.  And  frequently  the 
effects  of  mixing  are  so  large  that  second-order  perturbation 
theory  does  not  describe  them  adequately. 


Another  important  difficulty  inherent  in  the  usual  methods 
of  dealing  with  configuration  mixing  concerns  the  choice  of 
interacting  configurations,  it  is  of  course  true  that  the 
more  interactions  one  takes  into  account >  the  more  accurate 
will  be  the  theoretical  description.  It  has  been  amply  dem¬ 
onstrated,  however,  that  in  order  to  obtain  predictions  of 
term  spacing  in  neutral  and  slightly  ionized  systems  that 
are  accurate  to  within  a  few  percent,  one  needs  to  allow  for 
a  very  large  number  of  interactions.  Hence  a  calculation 
that  allows  for  mixing  between  a  small  number  of  configura¬ 
tions  would  be  significant  only  if  the  choice  of  the  configu¬ 
rations  considered  could  be  justified  on  grounds  other  than 
ma  thema  t ica 1  conven ience . 


Hie  problem  of  configuration  mixing  takes  on  a  differ¬ 
ent  complexion  if  we  adopt  the  Z-dependent  approach.  Con¬ 
sider  the  expansion 
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(C-4) 


i 


where  the  refer  to  distinct  configurations  and  the  sum 
runs  over  all  cc»ifigurations  connected  with  configuration  1 
(say)  by  nonvanishing  energy  matrix  elements,  (we  assume 
that  the  state  corresponding  to  has  been  nominally  assigned 
to  configuration  1) . 

The  numerical  values  of  the  coefficients  a^  depend  on 
the  choice  of  radial  wave funct ions .  For  example,  if  the 
radial  wavefunctions  in  are  solutions  of  the  Hartree- 
Fock  equations,  then  the  coefficients  a^  pertaining  to  con¬ 
figurations  with  the  same  set  of  azimuthal  quantum  numbers 
as  configuration  1  will  all  be  very  small,  because  the  vari¬ 
ation  principle  ensures  that  most  of  the  effects  of  mixing 
between  configurations  characterized  by  the  same  set  of  azi¬ 
muthal  quantum  numbers  have  already  been  incorporated  in  the 

wavefunctions.  Now,  every  coefficient  a.  is  either  0(1)  or 

-1  ^ 

0(Z  ) .  The  coefficients  of  the  first  type  are  associated 

with  configurations  that  have  the  same  set  of  principal 

quantum  numbers  as  configuration  1;  those  of  the  second  type 

with  configurations  that  have  a  different  set  of  principal 

quantum  numbers. 

The  coefficients  of  the  first  type  can  be  written  in 
the  form 

ai  =  a°  +  0(z'^)  ,  (C-5) 

0 

where  a^  is  independent  of  Z.  It  was  shown  in  A  that  the 
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constants  do  not  depend  on  the  form  of  the  radial  wave- 
functions  .  The  coefficients  and  parts  of  coefficients  that 
do  depend  on  the  form  of  the  radial  wavefunctions  are  0(Z~^) . 

We  may  conveniently  split  off  the  terms  on  the  right  side 
of  (  C-4)  that  have  constant  coefficients t 

■  1  “Ih  *  1 

i  i 

where 

aJ  =  a^  -  a?  *  0(z”^)  .  (C-7) 

The  nonvanishing  coefficients  a^  are  finite  in  number.  They 
can  be  calculated  with  arbitrary  precision,  being  essentially 
by-products  of  the  calculation  that  yields  the  energy  coeffi¬ 
cients  (see  A) . 

We  propose  to  replace  the  usual  approximation,  which  con¬ 
sists  in  neglecting  configuration  interaction  entirely,  by  an 
approximation  that  represents  the  wavefunction  for  a  given 
state  by  the  finite  sum 


y  -  ^  .  ,  (C-8) 

where  the  coefficients  may  be  regarded  as  Known  numbers.  (For 

2  cx 

spectroscopic  terms  belonging  to  complexes  of  the  type  (1  2^) 
and  (1^2*^)°  they  are  given  in  Table  C*l. 

The  approximation  just  described  becomes  increasingly 
accurate  as  z  increases  along  an  isoelectrMtic  sequence. 
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Even  for  moderately  ionized  and  un-ionized  systems,  however, 
it  should  prove  superior  to  the  central-field  approximation, 
since  experiment  has  shown  that  the  linear  Z-dependence  of 
the  term  energies  usually  dominates  the  structure  of  the  term 
spectrum. 

C.3.  The  Generalized  Screening  Parameters 

Notation.  The  letters  r,t,t'  will  denote  configu¬ 

rations;  the  symbol  P  will  distinguish  different  states  in  the 
same  configuration;  a  will  denote  an  orbital  nl;  and  a  will 
denote  a  complex.  Thus  the  label  a(r)r  distinguishes  the 
state  r  nominally  assigned  to  the  configuration  r  of  the  com¬ 
plex  a.  In  the  present  notation,  equation  (  C-4)  takes  the 
form 


y(arr) 


^  <arrltr>'i'(tr)  . 

t 


(C-9) 


The  Hamiltonian  H  for  a  many-electron  system  is  given, 
in  the  nonrelativistic  approximation,  by 

H  =  IIq  »  <C-10) 


where 


Hq  =  ^  {-p^/2m  -  ZeVr^)  , 
i 


(C-li) 


V  = 


(C-IZ) 
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The  variational  energy  in  the  present  approximation  is 
defined  as 

W  =  <arr'|n|  arr>  .  (C-13) 

We  now  evaluate  W  using  screened  hydrogenic  wavefunctions 
for  the  undetermined  radial  wavefunctions  that  figure  in 
this  formula. 

For  the  one-electron  part  of  W  we  obtain 

<arrlHQlarr>  =  <arri  tr>  <trlHQ|  t  •  r>  <t  Tl  arP) 
t.f 

=  KarritOl^  <trlHpltr> 

t 

t  a  (C-14) 

where  denotes  the  number  of  orbitals  of  type  u  in  the 
configuration  t,  and  s^^  denotes  the  screening  parameter 
that  corresponds  to  the  orbital  a  in  configuration  t.  A 
more  complete  notation  for  this  screening  parameter  is 
s  (at|rr).  The  suffix  Z  indicates  that  the  screening  para- 

z 

meter  depends  in  general  on  Z;  s^  differs  from  s  by  terms 
0(Z  .  The  labels  rl'  indicate  that  the  value  of  the  screen¬ 

ing  parameter  depends  on  the  state  under  consideration. 


For  the  two-electron  contribution  to  W  we  have 


Hence 


<arr|vlarr> 


y^<arr|tr>  <trlvlt* rt>  <t' r|arr> 


(C-15) 


<arr|H|arr>  »  ^Karrjtr)!^  ( 

t  a 

^<tr  arr> 


The  variational  equations  are 


(C-16) 


3w  Sxt^at  b 


2  <32  ^ 

at  n^  at 


«trlvltr>  +  2^Ctf;rr  °  * 

tyt  (c-17) 


where 


V  "  ^  "  ®at  '  ‘'tfjrr  “  <tT|arr>/  <tr|arr>  . 

(C-18) 

The  diagonal  matrix  element  ^tr|v|tr^  may  be  written  in 
the  form 


<tr|v|tr>  -  +  v^j,  , 


(C-19) 


where  the  contribution  represents  the  average  value  of 

<tr|v|tr>.  Explicitly, 


^<taSL|v|taSL>(2S  +  1)(2L  i-  1) 


aSL 


(2S  +  1)  (2L  1) 


(C-20) 


aSL 
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The  term-dependent  contribution  vanishes  when  the  con¬ 
figuration  t  consists  of  complete  groups  only,  or  of  complete 
groups  plus  or  minus  a  single  electron.  The  contribution  to 
v^p  from  interactions  between  electrons  in  a  group  l”  satis¬ 
fies  the  symmetry  condition 

v(l"r)  =  .  (C-21) 


In  a  first  approximation  one  may  wish  to  neglect  the  term- 
dependence  of  the  screening  parameters.  We  accordingly  define 


n^  <^V 

ot 

If  the  configuration  t  contains  no  more  than  one  inccxnplete 
group,  can  be  put  in  the  form 


where 


(a  =  p) 
(a  /  p) 


(C-23) 


(C-24) 


The  two-electron  interaction  (ap)  is  defined  in  A;  in  a  first 
approximation  (apl  ~  F®{ap) .  The  coefficient  denotes  the 

rmaiaer  of  ways  in  which  two  electrons  labelled  by  the  indicated 
quantum  numbers  can  be  selcted  from  the  configuration  t. 


With  the  help  of  these  definitions  we  can  write  s(at)  in 
the  form 


IM 


where 


s(alp)  =  (1  -  2  (aP)  •  <C-26) 

P  °*a 

The  two-electron  screening  parameter  s(a|p)  provides  an  unam¬ 
biguous,  qunatitative  measure  of  the  screening  of  an  a  elec¬ 
tron  by  a  p  electron. 

The  screening  parameters  s  (a|^)  are  independent  of  the 
configuration  t,  and  can  therefore  be  evaluated  once  and  for 
all. 


The  exact  screening  parameters  are  given  by 

2 

.(atlm  .  .(at)  t  ^  3I-  ,rrV'l 


tyt 


)  , 

(C-27) 


where 


Vtt  r  -  <tr|v|tT>  . 


(C-28) 


In  this  equation  one  is  not  in  general  justified  in  neglecting 

the  contribution  from  configuration  mixing.  In  some  instances  the 

second  term  in  the  brackets  (  )  may  be  considerably  larger  than 

the  first  term.  This  is  always  the  case,  for  example,  when  the 

configuration  t  is  composed  of  complete  groups  but  contains  one 

2  2 

or  more  incomplete  shells,  e.g..  Is  2s  .  Then  v^p  >  0,  but  ^ 


Ttie  aet  of  equations  (  C>27)  for  all  orbitals  a  in  the 
mixture  (  C-9)  determine  the  corresponding  screening  para¬ 
meters  s^iotlrP).  In  spite  of  their  complicated  appearance « 
these  equations  are  comparatively  easy  to  solve  numerically. 

I:,  practice  one  uses  an  iterative  procedure,  first  evaluating 
the  right  sides  for  Z  «  <».  This  gives  the  screening  constants 

s  (at|rr).  In  the  next  ai^roximation  one  evaluates  the  ri^t 
00 

sides  for  finite  values  of  Z>  replacing  the  parameters  s  by 
the  constants  s  ,  and  so  on. 

M 

One  can  take  advantage  of  the  dominant  contribution  of 

s(at)  to  s(at|rr)  to  shorten  the  iterative  evaluation  of  Z- 

dependent  screening  parameters.  One  begins  by  solving  the 

system  of  equations  (C-22)  for  s_(at).  One  then  evaluates 

z 

the  term-dependent  correction  to  s_(at|rr)  using  S-(at)  in 

z  z 

place  of  S2(at|rr)  in  equation  (  C>27).  With  this  method 
one  needs  to  go  through  the  iterative  procedure  only  once 
for  a  given  configuration. 


C.±i _ Generalised  Hartree-Fock  Theory 


The  Hartree-Fock  equations  may  be  written  in  the  form 


f  .  1  +  Ml  ^  1)  .  ^  “  ^nl 

‘  dr*  2r* 


nl 


E  P  -f 
nrm 


JU 

r 


<C-29) 


The  quantities  Y  ,  and  X  ,  are  functionals  of  all  the  radial 
nl  ni 

wavefunctions  appropriate  to  the  configuraticm  being  consid¬ 
ered.  The  exchange  term  has  its  origin  in  the  antisymmetry 
of  many-electron  wavefunctions. 


If  we  take  as  our  starting  point  the  a^roximation  (  c-9) 
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instead  of  the  central-field  approximation ^  we  arrive  at  an 
equation  that  is  formally  identical  with  (  C-29)  This  equa¬ 
tion  differs 4  in  fact,  from  the  ordinary  Hartree-Fock  equation 

only  through  the  exchange  term  X  , ,  which  contains  additional 

n  i 

contributions  arising  from  configuration  mixing „  These  con¬ 
tributions  are  of  precisely  the  same  type  as  the  ordinary  ex¬ 
change  contributions.  Thus  the  increase  in  accuracy  that  re¬ 
sults  from  allowing  for  the  dominant  effects  of  configuration 
mixing  has  not  been  bought  at  the  price  of  added  complexity. 

There  exists  a  very  close  analogy  between  the  Hartree- 
Fock  theory  and  the  hydrogenic  screening  theory,  which  it  is 
convenient  to  exploit  because  the  screening  theory  is  so  much 
simpler  than  the  Hartree-Fock  theory.  The  generalized  Hartree- 
Fock  equations  are  analogous  to  the  screening  equations  (C-27). 
Just  as  in  the  screening  theory  one  is  not  justified  in  neglect 
ing  the  mixing  contributions  of  beside  the  term-dependent 

contributions  so  in  the  generalized  Hartree-Fock  theory 

one  is  not  justified  in  neglecting  the  mixing  contributions  to 

X  ,  , 

nl 

Calculations  based  on  the  generalized  Hartree-Fock  theory 
are  planned  for  the  near  future. 
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APPENDIX  D 


APPLICABILITY  OF  THE  SCREENING  THEORY  TO  HEAVY  ATOMS 
D.I.  Introduction 


In  the  following^  the  term  "screening  theory"  will  be 
used  to  designate  the  Z-dependent  approximation  scheme  descri¬ 
bed  by  Layzer  and,  in  a  generalized  form  that  allows 

fully  for  relativistic  effects,  by  Layzer  and  Bahcall 
The  term  "Slater  theory"  will  be  used  to  designate  the  con¬ 
ventional  theory  of  atomic  spectra,  whose  basic  approximations 
were  introduced  by  Slater  in  his  classic  paper  of  1929p'^^ later 
recent  two-volume  monograph,  "Quantum  Theory  of  Atomic  Structure", 
describes  the  modern  form  of  the  theory. 

The  screening  theory  was  devised  to  account  for  certain 
systematic  discrepancies  between  the  Slater  theory  and  experi¬ 
ment,  pertaining  to  the  centers  of  gravity  of  spectroscopic  terms. 
These  discrepancies  are  sometimes  large  enough  to  reverse  the 
predicted  order  of  term  energies  in  a  configuration. 

Now,  any  reasonably  general  theory  of  atomic  structure  muLt 
incorporate  some  fairly  drastic  simplifying  assumptions  The 
Slater  theory  rests  on  the  assumption  that  the  effects  of  con¬ 
figuration  mixing  can  be  neglected  in  a  first  approximation. 

The  screening  theory,  on  the  other  hand,  rests  on  the  assumption 
that  the  energy  of  an  atomic  system  can  be  expressed  as  a  rapidly 
convergent  power  series  in  (Z  -  s)~^,  where  Z  is  the  nuclear 
charge  and  s  is  a  screening  parameter.  This  assumption  leads 
naturally  to  a  theory  that  is  in  most  respects  simpler  than  the 
Slater  theory,  but  which  allows  for  certain  effects  of  configura¬ 
tion  mixing  in  the  first  approximation. 
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As  yet,  the  screening  theory  has  been  applied  only  to 
light  atoms  —  specifically  to  atoms  with  ten  or  fewer  elec¬ 
trons.  In  these  applications  it  has  proved  very  successful, 
accounting  in  a  quantitative  way  for  features  of  term  spectra 
that  the  Slater  theory  describes  only  semi-quantitatively  or 
qualitatively.  The  question  now  arises  whether  the  screening 
theory  will  also  provide  a  useful  description  of  the  spectra 
of  heavy  atoms.  The  answer  hinges  on  how  rapidly  expansions 
of  the  term  energies  in  inverse  powers  of  Z  -  s  converge . 

The  fact  that  Z  -  s  increases  with  N  (the  number  of  atomic 
electrons)  argues  for  the  applicability  of  the  theory  to 
heavy  atoms.  On  the  other  hand,  one  would  also  expect  the 
coefficients  in  the  energy  expansion  to  increase  with  N, 

Since  it  is  far  from  obvious  to  what  extent  the  second  effect 
will  offset  the  first,  the  question  cannot  be  decided  with 
confidence  on  general  grounds;  specific  calculations  are 
required. 

The  iron-group  elements  provide  ideal  subjects  for  such 
calculations.  They  are  sufficiently  heavy  to  afford  a  defini¬ 
tive  test.  At  the  same  time,  they  are  sufficiently  light  to 
enable  relativistic  contributions  to  the  energy  to  be  treated 
as  perturbations.  Finally,  the  spectra  of  many  atoms  and 
ions  in  the  iron  group  have  been  well  analysed  experimentally. 

In  evaluating  the  merits  of  a  theoretical  description, 
one  must  compare  it  not  only  with  the  data  it  sets  out  to 
describe  but  also  with  other  available  theoretical  descriptions. 
During  the  past  few  years,  attempts  have  been  made  to  patch  up 
the  Slater  theory,  as  it  applies  to  the  spectra  of  elements  in 
the  iron  group,  by  adding  to  the  theoretical  energy  formulae 
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an  empirical  correction  term  of  the  form  aL(L  -t-  1)  .  That  a 
term  of  this  type  could  account  for  a  fairly  large  part  of  the 
discrepancy  between  the  Slater  theory  and  experiment,  in  spectra 
of  the  iron  group,  was  first  pointed  out  by  Layzer  and, 

independently,  by  Racah.^’'"  '  The  Slater  theory,  amended  in 
this  way,  has  since  been  used  extensively  by  Trees,  Rohrlich, 

and  Racah  and  his  pupils  in  semi-empirical  analyses  of  term 

(r-7) 

spectra  in  the  iron  group.  Racah  has  suggested  that  a 

theoretical  basis  for  the  L(L  -f  1)  correction  is  afforded  by 
the  work  of  Bacher  and  Goudsmit.^^’**^ 

In  the  following,  we  shall  re-examine  the  theoretical 
basis  of  the  L(L  t  1)  correction.  We  shall  also  present  the 
results  of  new  calculations,  based  on  the  screening  theory, 
for  complexes  with  a  partly  filled  M  shell.  We  shall  find 
that  the  modified  Slater  theory  cannot  account  quantitatively 
for  the  effects  of  configuration  mixing  described  by  the  screen¬ 
ing  theory.  We  shall  also  find  that  the  latter  are  of  the  same 
order  of  magnitude  as  the  observed  discrepancies.  This  suggests 
tluit  i:hc  screening  theory  will  provide  an  adequate  description 
of  these  spectra.  Unfortunately,  sufficient  experimental  data 
are  not  yet  available  to  test  the  theoretical  calculations 
presented  here,  the  complexes  here  studied  having  been  >liosen 
for  their  comparative  mathematical  simplicity  rather  than  for 
comparison  with  experiment.  However,  work  now  in  progress 
should  lead  before  long  to  a  definitive  experimental  test  of 
the  screening  theory  as  applied  to  atmns  and  ions  in  the  iron 
group . 
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D.2.  Theoretical  Results  for  Complexes  with  a  Partly-Filled 
M  Shell* 

2  8  2 

We  have  considered  the  following  complexes:  (123), 

(1  2  3  )  4  (1  2  3  ),  (1  2  3  )  .  The  configurations  belong¬ 

ing  to  these  complexes  are  listed  in  Table  D-I.  These  are  the 
simplest  complexes  with  partly-filled  M  shells ^  other  than 
those  consisting  of  closed  shells  plus  or  minus  a  single  elec¬ 
tron.  Complexes  containing  from  three  to  seven  3d  electrons 
present  a  somewhat  more  difficult  theoretical  problem. 

The  evaluation  of  the  energy  coefficients  and  the 

unitary  matrices  that  diagonalize  the  energy  submatrices 

V/  \  proceeds  along  the  lines  laid  down  in  the  reference 
(n)pSL  fn-l) 

cited  earlier  (Layzer  '  *  is  the  coefficient  of  Z 

in  an  expansion  of  the  nonrelativistic  part  of  the  energy  in 

inverse  powers  of  Z;  V,  .  is  a  finite  matrix  whose  elements 

( n  /  psit 

connect  states  labelled  by  the  indicated  quantum  numbers 
is  an  eigenvalue  of  the  radial  integrals  that  figure 

in  the  matrix  elements  being  evaluated  with  straight  hydro- 
genic  wavefunctions . ] 


The  results  of  the  calculation  are  present  in  Tables  D-2  thru 

D-9.  Tables  D-2  thru  D5 give  the  coefficients  Wj^  as  well  as  the 

values  that  would  have  been  obtained  for  these  coefficients 

if  configuration  interaction  had  been  neglected  (i.e.,  the 

values  of  according  to  the  Slater  theory)  .  Tcbles  D-6  t’  ru  D-9 

give  the  unitary  m-' trices  that  diagonalize  the  various  energy 

subroatrices  V ,  .  . 

(n)pSL 


*  The  work  summarized  in  this  subsection  was  supported  in  part 
by  a  grant  from  the  National  Science  Foundation. 
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TABLE  D-1 


TBIMB  AND  CONPIGUIATXONS  FOR  COMPLEXES  NXTH  A 
PARTIALLY  PILLED  M  SHIU 


I 


19/ 


TABLE  D-2 


wj  In  (iVa^). 

Is  the 

term*dependent 

part  of  Wj^; 

-  Wj  +  Wj 

1 

(W, ,  denotes 

1  C  .  I . 

the  value 

1 

of  obtained 

when  configuration  mixing 

is  neglected. 

1 

1 

f  1 

Configuretlon 

Term 

W 

"l 

<«l>c.f. 

V<“l>c.£. 

3s^ 

-.00357 

0 

-.003S7 

3s3d 

-.00228 

-.00228 

- 

.01581 

.00684 

.00897 

3P^ 

-.00481 

-.00432 

-.00049 

-.00560 

.00432 

-.00992 

.01748 

.01727 

.00021 

3d^ 

-.00564 

-.00564 

- 

.00295 

.00201 

.00094 

.00373 

.00322 

.00051 

.00616 

.00616 

- 

's 

.02719 

.02382 

.00337 

TABLE  D-3 


In  (1^3^)°  (See  Table  D-2) 


•  t  II 


Config. 

Term 

U 

"l 

(Wj)  c.f. 

W^-(Wj)c.f 

3r.  3p 

.01801 

.02116 

-.00315 

V 

-.00753 

-.00705 

-.00048 

3p3d 

.02215 

.01900 

.00315 

.00255 

.00207 

.00048 

'd° 

-.00869 

-.00869 

- 

o 

a 

.00086 

.00086 

- 

.01934 

.01934 

- 

-.00859 

-.00859 
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TABLE  0-4 


'  in  (1V3I®: 

Configuration  Term 

”1' 

3p  3d 

's 

.00262 

38  3p^  3d® 

^0 

-.00228 

^0 

.01926 

3s^  3p^  3d^° 

-.00406 

lo 

-.00763 

's 

.01625 

38^  3p^  3d® 

-.00564 

^0 

.00154 

. 00296 

.00616 

's 

.02226 

(See  Table  0-2) 


<Wl’>cf 

0. 

+.002620 

-.00228 

.... 

.00684 

+.012424 

-.00432 

+.000258 

.00432 

-.011951 

.01727 

-.001024 

-.00564 

.00201 

-.000472 

.00322 

-.000258 

.00616 

.... 

.02382 

-.001595 
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TABLE  D-5 


Wj[  In  (1^2®3^®).  (See  Table  D-2) 


Configuration  Term 

w  ' 

“l 

<^l'>cf 

3 

1 

SB 

,  ,5  -.10 

3s  3p  3d 

.02281 

.02116 

.001649 

-.00678 

-.00705 

.000274 

15  9 

38^  3p^  3d^ 

.01735 

.01900 

-.001649 

.00180 

.00207 

-.000274 

-.00869 

-.00869 

____ 

.00086 

.00086 

.01934 

.01934 

-.00859 

-.00859 

I 

I 


cf 
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TABLE  D-6 


THE  UNITARY  MATRICES  |  |<^rr|tl>|  |(lV3^) 
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TAIU  D-7 

TW  UMlTAtY  MAItlCIS  |  l<^rr|t^j  (,(1*2*3*)® 


1, 


TABLE  D-8 


I 


THE  UMITARY  MATRICES  |  |<trr|tf>j  ( 
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TABLE  D-9 


THE  UNITARY  KATRICES  |  IcirF  |  tr>!  |  ,(1^2®3^^)® 

IP® 

383p^3d^° 

3®^3p^3d® 

3s3p^3d^° 

1 

.99770 

-.06782 

'259 

,3s^3p^3d’' 

1 

.06762 

.99770 

3P° 


3s3p^3d^° 

2  5  9 

38^3p  3d’^ 

3s3pV° 

.99960 

-.02818 

38^3p^3d^ 

.02818 

.99960 

We  call  attention  to  the  following  features  of  the  results. 


(1)  In  the  complexes  considered  here,  the  effects  of  con¬ 
figuration  mixing  within  a  given  complex  are  comparable  in  mag¬ 
nitude  with  the  corresponding  effects  in  complexes  of  the  form 
(1  2^)  studied  previously.  Thus,  the  mixing  coefficients,  which 
indicate  the  magnitude  of  departures  from  the  central  field 
approximation,  are  of  the  order  of  .1  and  the  changes  in  the 
intervals  between  adjacent  terms  are  of  the  order  of  10-30%. 

(ii)  In  a  few  instances  the  changes  may  be  much  greater 

than  the  values  mentioned  in  the  last  paragraph.  For  example, 

2  2  8  2 

the  configurations  3s 3d- and  3p  in  the  complex  (12  3)  are 

very  strongly  mixed,  the  mixing  coefficients  being  about  ,44. 

The  same  ia  true  of  the  corresponding  pair  of  configurations 

2  8  16 

in  the  complex  (12  3  ),  where  the  mixing  coefficients  are 

about  .55.  In  these  cases  the  resulting  change  in  the  value 
of  W.  is  great  enough  to  change  the  original  order  of  the 
spectroscopic  term  (for  large  Z) .  Thus  W^(  D)  actually  lies 
below  the  value  of  W^(^P)  in  the  configuration  3p^  of  (1^2^3^) 
when  allowance  is  made  for  configuration  mixing  within  the 
complex.  The  change  in  for  this  term  thus  exceeds  100% 
of  the  separation  »  Wj(^D)  -  Wj^(^P)  as  predicted  by  the 
Slater  theory. 

(ili)  It  is  necessary  in  general  to  allow  for  all  effects 
of  configuration  mixing  within  a  given  complex,  though  in  parti¬ 
cular  cases  some  ccMitributions  may  be  negligible.  In  particular, 
.interactions  connecting  configurations  that  differ  in  the  quantum 
numbers  of  a  single  electron  cannot  in  general  be  neglected  be- 
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aid*  interactions  that  connect  configurations  that  differ  through 

the  quantuHi  nundDers  of  two  electrons.  The  configurations  3s3d 
2 

and  3p  discussed  in  the  last  paragraph  Illustrate  this  point. 

(iv)  The  effects  of  configuration  mixing  are  smaller  in 
the  complex  (1^2^3^^)  than  in  (1^2^3^).  This  beheirior  is  similar 
to  that  observed  in  complexes  of  the  type  (12^) t  The  effects 
of  configuration  mixing  tend  to  diminish  with  increasing  q.  The 
reason  is  that  the  off-diagcmal  elements  of  electrostatic  inter- 
action  are  essentially  independent  of  the  number  of  electrons  in 
the  system,  whereas  the  differences  between  diagonal  elements 
of  the  energy  matrix  increase  with  increasing  H  (total  nusiber  of 
electrons) .  Since  the  mixing  coefficients  are,  in  a  first  approxi¬ 
mation,  equal  to  the  ratios  bet%ireen  off-diagonal  interacticm 
elements  and  diagonal  energy  differences,  they  will  tend  to  diminieh 
with  increasing  N. 

P.l.  The  i^^^-ftonirical  Method  of  Bacher  and  Goudsmit 

The  method  of  Bacher  and  Goudsmit  seeks  to  relate 

the  term  energies  of  an  N-electron  atom  with  the  term  energies 
of  its  ions,  so  that  if  the  appropriate  term  energies  of  systMM 
(N,Z)  with  N  <  N  are  known  (from  either  theory  or  experiment) 
one  can  calculate  the  term  energies  of  the  system  (li,Z)  to  a 
certain  approximation. 

The  theoretical  relations  connecting  the  N-electron  ener¬ 
gies  with  the  1-,  2-,  ...  ,  (N  -  1) -electron  energies  are  based 
on  the  ii^pendent-orbltal  approximation.  That  is,  each  electron 
is  assigned  a  definite  pair  of  orbital  quantum  numbers  nl.  We 


207 


now  derive  these  relations,  using  a  more  explicit  notation 
than  that  of  Bacher  and  Goudsmlt. 

Let  denote  the  nundser  of  Is  orbitals,  the  number  of 
2s  orbitals,  q  the  number  of  2p  orbitals,  and  so  on  up  to  q  , 

J  IT 

so  that  for  an  N~electron  atom 

**1  **■  **2  “  **  * 

Let  denote  the  set  of  quantum  numbers  that  are  required.  In 
addition  to  the  q^,  to  specify  the  state  of  an  N-electron  atom. 

Consider  a  given  state  An  M~electron  subsystem, 

consisting  of  the  nucleus  and  M  <  H  electrons,  will  not  In 

1  M  1 

general  be  In  a  definite  state  (a  ...  a  g  denote 

pairs  of  orbital  quantum  nunbers,  not  necessarily  all  distinct, 

selected  from  the  set  of  N  pairs  that  characterise  the  N- 

electron  state  under  consideration.  We  shall  use  the  abbrevl- 
1  M 

at  Ion  (a  . . .  a  )  «  Thus  wo  may  use  the  symbol  In 

place  of  (qj^  ...  With  every  state  we  may  associ¬ 

ate  a  definite  probability  P(a„r  |a„r  ) .  Hore  generally,  we 
consider  the  set  of  all  probabilities  of  the  type  > 

where  A  <  B  <  N  and  the  set  Is  contained  In  the  set  a^. 

These  probabilities  are  not  all  Independent,  since  they  evi¬ 
dently  satisfy  all  relations  of  the  form 

^f(r^|rg)»tr,|rc)  -  Kr^ir^)  •  (b-D 

B 

(For  sls^llclty  we  have  omitted  the  labels  a.) 

By  means  of  what  are  now  standard  methods  In  the  theory  of 
complex  spectra,  one  can  evaluate  any  probability  • 
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Let  W(a„r  |a„r„)  denote  the  energy  of  the  M-electron  sub- 
M  N  N  N 

system  in  state  P  .  considered  as  part  of  the  N-electron  system 

M 

in  state  r„.  Then 

N 

''m 

is  the  mean  energy  associated  with  the  subconfiguration  a 

N 

when  the  N-electron  system  is  in  state  P  .  Note  that  the  ener- 

N 

gies  W(a  P  |o„P„)  do  not  coincide  with  the  energies  of  the  cor- 
responding  M-electron  ions  of  the  N-electron  system  which  would 
be  denoted  by  W(a  P  ja  P„),  because  the  electrons  in  the  sub- 
system  move  in  the  combined  field  of  a  nuclear  charge  Z  and  N-1 
orbital  electrons,  while  the  electrons  in  the  ion  move  in  the 
combined  field  of  the  same  nuclear  charge  and  M-1  orbital  elec¬ 
trons.  Thus  the  energy  levels  of  the  subsystem  lie  above  the 
corresponding  levels  of  the  ion.  However,  in  practical  appli¬ 
cations  of  the  method  to  be  described,  Bacher  and  Goudsmit  used 
empirical  ionic  energies  as  approximations  to  the  corresponding 
subsystem  energies.  The  latter  cannot,  of  course,  be  obtained 
directly  from  experiment.  We  shall  return  to  this  question 
later . 


Following  Bacher  and  Goudsmit,  we  write  the  energy  of  an 

M-electron  subsystem  as  a  weighted  sum  of  contributions  « 

K  <  M,  where  a„  stands  for  a  set  of  K  labels  selected  fr<m 
^  1  M 

the  set  of  M  labels  a  =  a  . . .  a  .  The  energy  contributions 

N 

w  are  defined  implicitly  by  the  set  of  equations 


“'Vm*  -  I""*’  *  *  •••  * 


ft 


ft. 
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•  •  • 


(0-2) 


Th«  Man  energy  contributions  W  are  defined  in  analogy  with 
the  Man  energies  Wt 

"k 

fhe  set  of  all  these  equations «  for  all  states  with  M  <  N, 

has  a  unique  solution  for  the  energy  contributions  w  in  terms 
of  the  energies  Wt 

‘^M-1  “m-2 


-  (-1)' 


r* 


2^W(a) 


(D-4) 


In  this  equation  all  energies  pertain  to  subsystems  of  the 
H-electrM)  system  in  the  state  The  proof  of  the  formula 

is  by  MtheMtical  induction.  For  M  -  1,  the  formula  reduces  to 

W(a)  »  w(a)  ,  (D-5) 

which  coincides  with  equation  (D-2)  for  M  ■  1.  Suf^se  now 
that  the  formula  holds  for  N  >  1 ,  . . .  ,  n .  We  now  show  that 
it  must  then  hold  for  M  >  n-i-l.  Prom  equation  (D-4)  we  have. 


a  r 
n  n 


a  ,r  , 

n-1  n-1 
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a  ,  r  , 
n-1  n-1 


-  I‘“K-/n-2lVn’*'‘“n-2''n-2>  *  •••’■••  ' 

“11-2^11-2  (D-6) 

Using  the  probablllty-ccMBposItlon  rule  (I^IV,  we  car.  put  the 
right  side  of  this  relation  into  the  form 


n 


''*“n+l''n+l’  *  Z‘^k^"*“n+l-k'“n+l'^n+l*  ' 

•‘-I 

The  number  of  individual  contributions  to  the  kth  summand  is 
C^(  ^  ) .  This  number  is  also  given  by 


(D-7) 


k  k-t-1 

j»0  j«0 

Using  the  addition  theorem  for  binomial  coefficients,  we  obtain 


k+1 

z 

j-0 


,  . »k+l/  0  .  - 

<k+i>  ■  ®  • 


(D-9) 


Hence 


<=k  -  ' 

which  cosipletes  the  proof. 


(  -10) 


Bacher  and  Goudsadt  assume  that,  for  a  given  state  an^n* 
the  energy  contributions 

^w(o^lojjr^) 
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diminish  rapidly  with  k.  Consider  the  expansion 


■  1"'“’  *  I  "SIVn’  I  '"<“kl“NV 

“  “2  “k 

+  «  .  *  +  * 

which  is  a  special  case  of  equation  (D-2) .  By  omitting  the 
last  term  on  the  right  we  obtain  an  approximate  formula  for 
the  energy  of  the  N-electron  system  in  terms  of  quantities 
pertaining  to  k-electron  subsystems  with  k  ^  N.  From  equa¬ 
tion  (D-4)  it  follows  that  the  approximation  w(a„r  )  »  0  is 
equivalent  to  the  approximation 


"<Vh>  ■  • 


a. 


a 


*N-1  N-2 

Ulus  the  energy  of  the  M-electron  system  is  expressed  as  a 
linear  condoination  of  energies  of  k-electron  subsystems  with 
k  <  N. 


More  generally,  we  obtain  the  kth  approximation  in  the 
sense  of  Bacher  and  Gioudsmit  when  we  retain  the  first  k-fl 
terms  on  the  right  side  of  equation  ('P-11).  Ii  particular, 
the  first  approximation  (which  is  exact  for  2-electron  sys¬ 
tems)  is  defined  by  the  equation 


(/P-11) 


(D-12) 


('P-13) 
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/V) 

with  the  help  of  equation  (D-4)  we  can  express  W  ,  for 
any  k  <  Nj  as  a  linear  combination  of  energies  pertaining 
to  (k+l) -electron  subsystems. 

In  practical  applications  of  the  Bacher-Goudsmit  approxi¬ 
mation  method  one  uses  experimental  values  of  the  ionic  ener¬ 
gies  W(a.r.  la.r.  )  in  place  of  the  subsystem  energies  W(a.r.  |a„r  ). 
In  order  for  this  to  be  a  good  approximation^  the  electrons 
that  are  present  in  the  M-electron  system  but  not  in  the  k- 
electron  subsystem  must  make  only  a  minor  contribution  to  the 
electrostatic  field  seen  by  the  remaining  electrons.  Since 
the  removal  of  "core  electrons"  (electrons  in  closed  shells) 
does  substantially  alter  the  effective  electrostatic  field  of 
the  atom,  one  should  always  regard  the  core  electrons  and  the 
field  they  produce  as  fixed.  That  is,  one  should  reckon  all 
energies  from  the  core  energy  as  a  fixed  zero-point.  The 
formulae  already  developed  apply  without  change  to  these 
relative  energies,  since  they  do  not  depend  on  any  assump¬ 
tions  regarding  the  nuclear  contribution  to  the  electrostatic 
field  except  the  assumption  of  radial  symnetry.  Since  the 
core,  by  definition,  is  conposed  of  coeqplete  groups  of  elec¬ 
trons,  its  c<mtributi«i  to  the  field  will  have  spherical 
symmetry  and  may  therefore  be  lusqped  with  the  nuclear  con¬ 
tribution  . 

_ The  Bacher-Goudsmit  Method  and  Configuration  Mixing 

The  independeht-orbital  approximation  is  just  as  basic 
in  the  theory  of  Bacher  and  Ooudsmit  as  in  the  theory 'of 
Slater.  The  basic  equations  of  the  theory,  (D-2)  and  (D-4) , 
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presuppose  the  possibility  of  assigning  every  subsystem  of  an 
M-electrcm  system  to  a  definite  electronic  configuration,  and 
the  coefficients  in  these  equations  are  to  be  evaluated  by 
means  of  a  theory  based  on  the  same  approximation.  Neverthe¬ 
less,  the  Bacher-Goudsmit  approximation  does,  in  a  sense, 
allow  for  certain  effects  of  configuraticMi  mixing,  because 
it  makes  use  of  empirical  ionic  energies,  since  these  include 
contributions  from  configuraticm  mixing,  so  roust  the  "predicted" 
energies.  Now,  the  approximation  of  greatest  practical  and 
theoretical  interest  in  the  Bacher-Goudsmit  scheme  is  the  first 
one,  which  yields  the  energies  defined  by  equaticm  (B-13) . 

The  question  therefore  arises t  Which  effects  of  configuraticm 
mixing  are  included  in  W^^^  when  the  two-electron  energies 
wCa^r^la^r^)  are  taken  from  experiment? 

Bacher  and  Goudsmit  discuss  this  question  briefly  and 
schematically,  using  ordinary  perturbation  theory.  They  con¬ 
clude  (incorrectly,  as  we  shall  see)  that  the  first  approxi¬ 
mation  includes  all  effects  of  configuraticm  mixing  with 
doubly  excited  configurations  (i.e.,  configurations  differ¬ 
ing  from  the  original  onm  in  two  pairs  of  one-electron  quan¬ 
tum  numbers) .  They  also  conclude  (correctly)  that  the  first 
apfsroximation  does  not  include  all  effects  of  configuration 
mixing  with  singly  excited  cmif igurations .  Before  turning 
to  a  detailed  consideration  of  these  questions,  however,  let 
us  examine  the  problem  from  another  angle. 

As  has  previously  been  pointed  out  (Layser  t\tm 

Slater  theory  predicts  certain  qualitative  regularities  in 
the  structure  of  term  spectra,  which  do  not  accord  with  the 
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observed  regularities.  For  example^  the  Slater  theory  pre¬ 
dicts  that  all  configurations  consisting  of  n  or  6-n  equiv¬ 
alent  p  electrons,  together  with  any  number  of  complete 
groups,  have  the  same  structure  apart  from  a  scale  factor 
Specifically,  if  n  =  2,  3,  or  4,  the  two  independent  inter¬ 
vals  between  the  spectroscopic  terms  with  L  -  0,  1,  2  are 
predicted  to  satisfy  the  relation  SD/DP  =  3/2.  The  Z- 
dependent  approximation  scheme  predicts  different  quali¬ 
tative  and  quantitative  regularities.  It  predicts,  for 
example,  that  the  interval  ratio  SD/DP  approaches  a  con¬ 
stant  value  with  increasing  Z  in  a  given  isoelectronic 
sequence,  and  that  this  constant  (which  can  be  calculated 
from  the  theory)  depends  both  on  the  principal  quantum 
numbers  of  the  p  electrons  and  on  the  complete  groups  pres¬ 
ent  in  the  electronic  configuration  How  does  the  Bacher- 
Goudsmit  theory  stand  on  this  matter? 

2  2  2 

Consider  a  simple  example,  the  configuration  Is  2s  2p 
in  the  carbon  sequence.  We  regard  the  K-electrons  as  com¬ 
posing  a  fixed  core,  and  proceed  to  calculate  the  term 
energies  relative  to  this  core  in  the  approximation  (0-13). 

We  obtain 

W'^'(2s^2p  SL)  «  W(2s  )  +  4W(2s2p)  +  W(2p  SL)  -  4W(2s)  -  4W(2p) ,  (D-14 

where  we  have  omitted  unnecessary  labels  on  the  right  side 

Note  that  the  mean  two-electron  energy  W(2s2p)  does  not  depend 

on  the  quantum  nundjers  SL,  because  the  two  2s  electrons  form 

a  complete  group.  From  this  equation  we  see  at  once  that  the 

2  2 

configuration  2s  2p  is  predicted  to  have  the  same  structure 

2 

as  the  configuration  2p  .  Actually,  the  interval  ratio  SD/DP 
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is  substantially  less  than  3/2  in  all  systems  belonging  to 

2  2 

the  beryllium  sequence  (i.e.,  in  all  configurations  Is  2p  ) . 
The  difference  between  the  two  limiting  ratios  is  due  to 
the  mixing  of  configurations  belonging  to  the  same  complex; 
both  limiting  ratios  are  correctly  predicted  by  the  Z-depen- 
dent  theory.  We  conclude «  therefore «  that  the  Bacher*>6oudsmit 
method  does  not  in  general  allow  correctly  for  this  type  of 
configuration  mixing  —  the  most  important  type  according  to 
the  Z-dependent  theory  —  in  the  first  approximation. 

Now,  in  the  example  just  given  the  "perturbing"  con¬ 
figurations  that  need  to  be  considered  if  .the  limiting  inter- 

4 

val  ratios  are  to  be  correctly  predicted  are  2p  (interacting 

2  2  2  2 
with  2s  2p  )  and  2s  (interacting  with  2p  ) .  These  are 

"doubly  excited"  configurations  in  the  sense  explained  above 

2 

(though  the  configuration  2s  should  perhaps  be  described 

2 

as  doubly  de-excited  relative  to  2p  ) .  It  follows  that  the 
first  approximation  in  the  Bacher-Goudsmit  scheme  does  not, 
in  fact,  allow  fully  for  the  effects  of  mixing  with  doubly 
excited  configurations. 

The  argument  given  by  Bacher  and  Goudsmit  fails  because 
it  does  not  allow  properly  for  the  Pauli  exclusion  principle. 
Bacher  and  Goudsmit  tacitly  assume  a  strict  correspcmdence 
between  doubly  excited  N-electron  configuraticms  and  doubly 
excited  configurations  in  which  the  principal  quantum  num¬ 
bers  of  the  excited  electrons  differ  from  those  of  the  orig¬ 
inal  electrons.  It  does  not,  however,  extend  to  doubly 
excited  configurations  with  the  same  principal  quantum  num¬ 
bers  as  the  original  configuration.  Thus,  the  de-excitation 
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2 

of  the  2p  electrons  in  the  configuration  2p  to  form  the 

2 

configuration  2s  is  permitted  by  the  exclusion  principle; 

2  2 

but  in  the  configuration  2s  2p  ,  de-excitation  of  the  2p 

4 

electrons  would  yield  the  non-permitted  configuration  2s  . 

The  predictions  of  the  first-order  Bacher-Goudsmit 
theory  differ  from  those  of  the  first-order  Z-dependent 
theory  in  other  important  respects.  For  example,  we  have 
seen  (see  subsection  2)  that  the  effects  of  mixing  with 
singly  excited  configurations  may  be  very  large.  We  have 
also  seen  that  in  some  instances  the  effects  of  configura¬ 
tion  mixing  are  too  large  to  be  adequately  described  by 
second-order  perturbation  theory  Final ly^  we  have  seen 

that  the  effects  of  configuration  mixing  are  systematically 
greater  in  almost-empty  shells  than  in  the  corresponding 
almost-full shells .  All  these  features  of  the  Z-dependent 
theory  are  at  variance  with  the  first-order  Bacher-Goudsmit 
theory . 

We  conclude  that  the  Bacher-Goudsmit  theory  does  not 
provide  a  satisfactory  freunework  for  dealing  with  the  most 
important  effects  of  configuration  mixing,  even  semi-em-. 
pirically . 

D.5.  The  Modified  Slater  Theory 

The  theory  of  Bacher  and  Goudsmit  is  held  by  Racah 
to  provide  a  theoretical  justification  for  the 
L(L  +  1)  correction  in  the  following  sense.  By  means  of 
equation  (D-13)  one  can  esqpress  the  energies  of  an  H-elec- 
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tron  system  in  terms  of  the  energies  of  two-electron  sys¬ 
tems.  Let  us  admits  for  the  sake  of  the  argument^  that  the 
dominant  effects  of  configuration  mixing  for  the  N-electron 
system  are  included  in  the  energies  calculated  with  equation 
(D-13)  if  one  uses  the  true  values  of  the  energies  W.a_r  la„r  J 
in  evaluating  the  right  side.  Instead  of  using  experimental 
values  for  the  last-named  energies,  one  could  use  an  interpo¬ 
lation  formula  containing  a  number  of  adjustable  parameters 
Suppose  we  are  considering  configurations  of  the  type  d”™ 

The  configuration  d  has  five  terms,  so  that  if  one  uses  an 
interpolation  formula  with  five  adjustable  parameters,  the 
values  of  the  parameters  will  be  uniquely  determined  by 
experiment.  Equation  (D-13)  then  enables  one  to  express 
the  energy  of  any  term  belonging  to  the  configuration  d*' 
in  terms  of  these  five  parameters.*  According  to  Slater's 
theory.^  the  mean  energy  of  all  the  terms  in  a  given  con¬ 
figuration  labelled  by  given  values  of  S  and  L  is  a  linear 
combination  of  certain  radial  integrals,  the  Slater  P  and  G 
integrals.  For  configurations  d"  there  are  three  such 
integrals;  P  ,  p^  and  P^.  However,  even  if  one  regards 
these  as  adjustable  parameters,  one  cannot  obtain  a  very 
good  representation  of  the  experimental  data.  The  addition 
of  two  empirical  correction  terms,  aL(L  •i’  1)  4-  PQ,  the  secMid 


*In  practice  the  parameters  are  not  regarded  as  having 
fixed  values,  corresponding  to  the  configuratim  d  ,  but 
as  functions  of  the  effective  central  field,  (hie  can 
evaluate  the  parameters  for  n  >  2  by  equating  the  pre¬ 
dicted  energies  to  their  experimental  counterparts  and 
solving  the  resulting  equations  by  the  method  of  least- 
squares  . 
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of  which  was  suggested  by  Racah^^'^^  makes  it  possible  to  represent  the 

2 

term  energies  in  the  configuration  d  exactly  and  to  represent  the  term 
energies  in  other  configurations  that  have  been  examined  with  fairly  good 
accuracy . 

From  the  discussion  in  the  preceding  subsections,  however,  we  know 

that  equation  (D-13)  does  not  provide  a  satisfactory  framework  for  dealing 

with  the  dominant  effects  of  configuration  mixing.  It  is  true  that  the 

empirical  description  of  spectra  in  the  iron  group  based  on  the  Slater- 

Racah  interpolation  formulae  has  been  remarkably  successful.  However, 

guided  by  the  considerations  of  the  present  paper,  one  can  easily  choose 

examples  of  configurations  for  which  the  method  will  give  poor  results. 

2  2  6  2 

A  simple  example  is  the  configuration  Is  2s  2p  3p  . 
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APPENDIX  E 
D.  Layzer 

RELATIVISTIC  EFFECTS  IN  HEAVY  ATOMS 

E. 1  Introduction 

In  heavy,  as  In  light,  atoms  the  chief  groupings  of  energy  levels 

have  their  source  in  the  nonrelativistic  (electrostatic)  interactions 

between  the  electrons.  These  groupings  are  characterized  by  sets  of 

prinicpr  I  quantum  numbers  (  Adjacent  groups  are  separated  by 

2  2  2 

intervals  of  the  order  of  Z  a  me  ,  where  Z  is  the  effective  nuclear 

e  e 

charge  for  the  outermost  electrons. 

The  structure  of  a  given  group,  however,  has  an  entirely  different 

character  in  heavy  atoms  .'nd  in  light  atoms.  In  the  latter,  each  group 

(or  complex)  splits  up  into  spectroscopic  terms  oSL  consisting  of  levels, 

disti  iguished  by  the  quantum  number  J.  The  spacing  between  terms  is  of 
2  2 

the  order  of  Z^a  me  ,  while  the  spacing  between  levels  in  a  given  term 

4  4  2  2  3 

is  of  the  order  of  Z^a  me  •  In  light  atoms,  by  definition,  a  Z^  «  1, 

so  that  the  fine  structure  is  clearly  differentiated  from  the  coarse 

structure.  The  fine  structure  originates  in  relativistic  interactions 

between  the  orbital  electrons  and  the  nucleus,  t*  c  course  .strnv  ture  in 

nonrelativistic  interactions  between  the  electrons. 

2  3 

In  heavy  atoms,  on  the  other  hand,  the  coupling  parasMter  a  Z^  is 
not  seutll.  Consequently  the  relativistic  interactions  cannot  adequately 
be  treated  by  first-order  perturbation  theory,  as  they  can  be  in  light 
atoms.  The  structure  for  the  swjor  groupings  (complexes)  Is  correspondli^ly 
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complicated,  in  general,  there  Is  no  clear-cut  distinction  between 
coarse  structure  and  fine  structure;  the  two  kinds  of  interactions  con¬ 
tribute  comparably  to  the  spacing  between  levels  in  a  given  complex. 

3  2 

When  Z  a  »  1>  the  relativistic  contributions  become  dominant; 
the  complex  breaks  up  into  sub  complexes,  characterized  by  quantum 
numbers  (n^j^).  Within  a  given  sub  complex  the  spacing  of  levels  is 
due  primarily  to  the  nonrelativistic  (electrostatic)  interaction. 

This  state  of  affairs  is  characteristic  of  x-ray  spectra.  Here 
may  be  comparable  to  Z  itself  because  the  missing  electrons  belong  to 
inner  shells. 

The  conventional  method  for  describing  the  structure  of  energy 

2  3 

spectra  in  the  intermediate-coupling  regime  (q  1)  is  set  forth  in 

(E-1) 

Theory  of  Atomic  Spectra.  It  is  based  on  the  central -fie Id  ap¬ 

proximation  and  allows  for  the  leading  relativistic  contribution  to 
the  energy  -  the  spin-orbit  interaction.  This  is  evaluated  for  electrons 
moving,  not  in  the  Coulon^  field  due  to  the  nuclear  charge  Z,  but  in  an 
effective  central  field  due  to  the  nuclear  charge  plus  the  electrons. 

The  conventional  method  has  several  drawbacks. 

It  is  better  suited  to  a  seml-eaq>irical  description  than  a 
purely  theoretical  one.  This  is  because  the  radial  integrals  (  ()  . 
that  determine  the  swgnitude  of  the  spin-orbit  interaction  cannot  be 
accurately  evaluated  by  purely  theoretical  smans.  They  are  best 
evaluated  by  treating  them  as  adjustable  parameters,  to  be  assigned 
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value*  that  will  make  the  theoretical  formulae  fit  the  data  as  well 
aa  poaslhle.  Thla  la  nearly  alwaya  the  adopted  In  practice. 

The  apln-orblt  Interaction  appear*  aa  a  aeparate  term  In  the 
Haadltonlan  for  a  many -electron  atom  only  In  the  Pauli  (low  Z)  ap¬ 
proximation.  When  Z  la  not  aauill,  It  la  by  no  maana  clear  that  the 
apln-orblt  Interaction  adequately  approximate*  to  the  relatlvlatlc 
part  of  the  Hamiltonian. 

An  alternative  deacrlptlon  of  relatlvlatlc  effect*  In  the 

atructure  and  apectra  of  many-electron  atom*  1*  provided  by  recent  work 

(E-2) 

of  Laycer  and  Bahcall  .■  Thla  la  an  extenalon  of  the  nonrelatlvlatlc 
Z-dependent  theory^^  The  principal  feature*  of  the  Z-dependent 

relatlvlatlc  theory  are  aa  follow*: 

(a)  The  Hamiltonian  employed  1*  ahown  to  contain  all  known 

relatlvlatlc  Interaction*  up  to  (but  not  Including)  thoae  of  order 
4  2  2 

a  Z  me  .  The  largeat  known  Interaction*  not  Included  In  the  Haadltonlan 
are  thoae  reaponalble  for  the  lead>  ahlft. 

(b)  The  apln-orblt  Interaction  la  allowed  for  exactly,  In 
the  xeroth  approximation  -  Juat  aa  In  the  Dirac  theory  of  the  hydrogen 
atom. 

(c)  The  theory  enable*  one  to  evaluate  all  contribution*  of 

.X 

order  Z  coa^ared  with  the  apln-orblt  Interaction.  Thla  la  not  poa- 
alble  with  the  conventional  theory  -  and  herein  Ilea  the  crucial  dif¬ 
ference  between  the  old  and  the  new  methoda.  It  1*  the  poaalblllty 
of  praclaely  evaluating  the  next-to-the-leadlng  relatlvlatlc  contrlbutlonv 
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that  renders  Che  new  neChod  quenCiCeClve  where  Che  old  one  wee  nerely 
seal -quenCi caC Ive . 

AlChough  Che  nexC-Co-Che  leading  relecivlacic  conCribuCions  can  be 

evelueCed  In  principle,  Che  ecCuel  celculacions  ere  long  end  difficulc. 

2 

Even  for  such  eoaperecively  siaple  elecCronic  configur scions  as  Is  2s 

2  (E-2) 

end  Is  3b,  considered  by  Leyser  end  Behcell  ,  Che  eveluacion  of  Chase 

Cema  required  a  long  celculecion.  For  heavy  aCoas  Che  corresponding 

ealculeCiona  are  ec  presenc  prohibicively  difficulc i  IC  is  Cherefare  of 

inceresc  Co  inquire  whecher  chare  is  an  approxiaecion  cnac  will  enable 

one  Co  esciaece  Che  conCribuCions  in  quesCion  wich  fair  accuracy.  Vm 

deal  wich  Chis  quaacion  in  Saccion  B.2.  In  SecCion  E.2,  we  consider  a 

relaCed  problea  (as  ic  curns  ouc  Co  be) -Che  relaciviscic  genera lisac ion 

(E*3) 

of  Che  nonralaciviscic  screening  paraaecere  defined  by  Layser' 


1.2  Annrexiaace  Fora  of  Che  aelacivlacic  Correcciona  for  Heavy  Atnaa 


2  4  2 

The  coefficienc  of  a  Z  ac  in  cha  sKprasaion  for  Che  energy  of  an 

(B-2) 

acoaic  scace,  according  co  Che  Cheory  of  Layser  and  Bahcall' 

2  3 

is  a  funccion  of  cha  couplii^  paraascer  X  •  a  Z  .  In  cha  excraaa  rola< 
Civlscic  case,  Chis  coefficienc,  which  we  denoce  by  a^,  is  given  by  a 
direcC  generalisacion  of  Soaaarfeld's  foraula: 


1^1  “l  ^  ' 


as  X  -•  00 


(B-1) 


In  cha  opposico  liaic  (X  0)  a^  does  |i2£  i^ioach  the  value  tbac  would 
be  prodicced  on  cha  basis  of  cha  eencral-fiald  approxiaacioni  co  Acain 
Che  correcc  liaicing  valuer  one  ausc  allow  for  cha  effaces  of  configuraCion 
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mixing.  (This  can  be  done  most  simply  by  using,  in  place  of  the  central* 
field  wave  function  \|r,  a  finite  linear  combination  of  wave  functions 
^a^\|r^,  the  constants  a^  being  given  by  the  screening  theory;  see 
Appendix  C. 

In  the  general  case  a^(X)  can  be  calculated,  without  excessive 
labor,  to  any  desired  accuracy  by  methods  described  in  the  paper  by 
Layser  and  Bahcall  cited  earlier.  We  may,  therefore,  consider  it  to 
be  known. 

4  3  2 

It  is  the  coefficient  a^(X)  of  the  term  a  Z  me  that  presents  the 
real  problem,  and  which  we  would  like  to  approximate. 


Now  in  the  standard  theory  of  one-electron  spectra  (including  x-ray 
spectra),  the  relativistic  contribution  to  the  energy  is  usually 
written  in  the  form 


a^(Z  -  a)  a  me  +  0(a  Z  ) 


(E-2) 


where  a^  is  given  by  Sommerfeld  's  formula. 


1  /  r. 

.2n" 


3 

4 


and  a  is  a  suitably  chosen  constant.  Evidently 

a. 


aj  -  -4oa^;  a  - 


4a, 


(E-1') 


(E-3) 


Thus  a  knowledge  of  o  is  equivalent  to  a  knowledge  of  a^. 


Mow  soaw,  though  not  all,  of  the  contributions  to  a^  can  be  inter¬ 
preted  physically  in  terms  of  screening.  For  example,  tlM  spin-other- 
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orbit  interactions,  which  contribute  to  a^,  tend  to  reduce  the  magnitude 
of  the  (true)  spln-orblt  Interaction  In  much  the  same  way  as  the  Inter- 
electronlc  electrostatic  Interactions  reduce  the  magnitude  of  the  electro¬ 
static  Inter  ctlons  between  the  nucleus  and  the  electrons.  On  the  other 
hand,  such  two-electron  interactions  as  the  spin-spin  interaction  have 
no  one-electron  counterpart.  The  interpretation  of  a  as  a  screening 
parameter  is  therefore  an  oversimplification.  Nevertheless,  it  may 
furnish  a  useful  means  of  estimating  the  value  of  a^  through  Equation  (L-3). 

The  simplest  assumption  one  can  make  about  a  is  to  equate  it  to  S, 
the  screening  parameter  that  figures  in  the  screening  theory.  A  com¬ 
parison  of  emprical  values  of  a  for  hydrogenic  spectra  with  theoretical 
values  of  S  shows  that  this  a  a  fair  approximation,  though  not  as  good 
as  one  might  wish. 

Fortunately,  the  fractional  errors  due  to  inaccuracies  in  the 
estimated  value  of  o  tend  to  become  progressively  smaller  with  increasing 
Z  because  the  quartic  term  becomes  progressively  larger  in  comparison 
with  the  cubic  term.  Thus,  the  screening  approximation  to  a^  should 
be  considerably  better  for  Uranium,  say  than  for  Lithium. 

A  semi -empirical  study  to  improve  the  estimate  of  a  might  proceed 
along  the  following  lines: 

The  two-electron  screening  parameter  S(u/3)  [representing  the 

screening  of  an  a  electron  by  a  ^  electron]  is  given  in  a  first  ap- 

fE-3) 

proximation  by  the  formula'  ' 


226 


dr 


(E-4) 


OD 


p;  <r)  f 

j-v- 


S(a/p)  - 


(rS  dr^ 

P 


00 

I 


P^(r) 


dr 


This  can  be  written  as 


S(a/^) 


i.e.  as  a  certain  weighted  average  of  the  probability 
a  0  electron  will  be  found  closer  to  the  nucleus  than 
of  an  cx  electron.  The  weight  w  is  given  by 


(E-5) 


r  I  the  distance 
a 


w<r^ 


P^<r  )r  ■' 

a  a  a 


(E-6) 


This  weighting  factor  is  appropriate  for  electrostatic  screening,  but 
not  necessarily  for  other  kinds  of  screening.  A  more  general  law  than 
(E-6)  Is: 


"a 


(B-7) 


It  would  be  Interesting  to  test  whether,  by  an  appropriate  choice  of  k, 

one  could  arrive  at  a  formula  for  o  In  better  agreement  with  experiment. 

[The  connection  between  the  total  screening  constants  and  the  two- 

(E-3) 

electron  screening  parameters  Is  discussed  In  the  paper  by  Layxer' 
cited  earlier.  ]  In  the  meantime,  the  old  definition  (S  ■  o)  can  be 
used  with  a  reasonable  expectation  of  fair  accuracy. 
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E.3  Relativistic  Screening  Parameters 


The  screening  parameters  that  we  have  dealt  with  so  far  are  defined 
by  a  nonrelativlstic  theory.  The  question  naturally  arises  as  to  what 
constitutes  the  proper  relativistic  generalization  of  the  screening 
parameters. 

To  begin  with,  the  relativistic  screening  parameters  must  depend 
on  the  coupling  parameter  X  as  well  as  on  Z.  They  must,  in  fact,  have 
the  form 

S(X,Z)  -  S  (X)  +  S  ,(X)Z‘^  +  0(Z’^).  (E-8) 

In  the  limit  X  -•  0,  S(X,Z)  -*  the  nonrelavlstlc  screening  parameter  S(0,Z). 
The  coefficients  S^(X)  are  all  bounded  from  above  and  below.  In  the 
relativistic  limit  X  -•  oo, 

S(X,Z)  -*  S^(oo)  (E-9) 

since  the  remaining  terms  on  the  right  side  of  (a)  vanish  in  this 
limit. 

It  is  possible  to  develop  a  detailed  theory  of  the  X-dependent 
screening  parameters  along  the  lines  of  the  variational  nonrelativlstic 
theory.  The  main  result  is  that  the  contributions  to  the  screening 
paranmters  from  electrons  in  closed  shells  are  independent  of  X« 

Now  the  X-dependence  of  the  screening  parameters  could  be  significant 
only  in  heavy  atoms.  Since  most  of  the  electrons  in  a  heavy  atom  are 
norswlly  found  in  closed  shells,  the  amplitude  of  variation  of  the 
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screening  parameters  with  X  must  be  small,  so  that  the  nonrelativlstlc 
approximation  to  the  screening  parameters  should  be  adequate  for  most 
practical  purposes. 
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APPENDIX  F 

D.  Layzer  and  A.  Naqvl 


STATISTICAL  DESCRIPTION  OF  COMPLEX  SPECTRA 


F. 1.  Mean  Energies  of  Configurations  and  Complexes 

In  any  partial  description  of  complex  energy  spectra  the 
mean  energies  of  configurations  and  complexes  play  a  central 
role.  Even  if  no  other  information  about  the  energy  spectrum 
is  available  one  can  still  obtain  from  the  mean  energies,  with 

i 

the  help  of  Saha's  equation,  useful  estimates  for  the  concen¬ 
trations  of  the  various  atomic  and  ionic  species  in  a  mixture 
of  gases.  Again,  in  constructing  radial  wavefunctions  for 
use  in  the  evaluation  of  absolute  transition  .probabilities  — 
either  by  the  screening ■ approximation  or  by  a  more  refined 
method  —  one  may  not  wish  to  use  a  different  set  of  wave- 
functions  for  each  distinct  level.  In  a  first  approximation 
one  could  then  use  a  single  set  of  wavefunctions  correspond¬ 
ing  to  the  mean  energy. 

If  we  neglect  entirely  the  effects  of  configuration  mix¬ 
ing,  we  can  write  the  mean  interaction  energy  of  any  configu¬ 
ration  as  a  sum  of  two-electron  contributions  [nl,n'l').  If 
the  atom  has  N  electrons  the  mean  interaction  energy  consists 

of  (")  s  N(N  -  l)/2  such  contributions.  For  example,  the 

^  2  2  5 

mean  energy  of  the  configuration  Is  2s  2p  is 

(Isis)  •••  (2s2s)  ^  10(2p2p)  +  4(ls2s)  10(ls2p]  -i-  10(2s2p)  .  (F-1) 

Here  N  >  9  and  the  number  of  two-electron  interactions  is 
9  X  8/2  »  36. 

The  two-electron  interactions  can  be  expressed  in  terms  of 
the  Slater  F  and  G  Integrals.  The  appropriate  formulae  are 


231 


given  by  Layzer  and  Slater  but  are  In^llclt 

in,  earlier  work  by  Condon  and  Shortley  and  Shortley. 


To  calculate  the  total  mean  energy  of  a  configuration, 
including  the  one-electron  part,  we  may  conveniently  eBg>loy 
the  Z-dependent  screening  theory.  The  mean  energy  is  then 
given  by  the  formula 

Here  is  the  number  of  electrons  of  the  type  nl  in  the 
configuration  under  consideration  and 

*nl  -  *  -  ‘nl  • 

The  screening  parameters  s^^  are  defined  in  terms  of  the 
t%#o-electron  interaction  energies: 

•  ,  -  y  (q  .1.  -  *  I  ,)B(nl|n'r)  ,  (P-4) 

nV  ^n'l*  nl,n'l'  »  •  »  ^ 

n'r 

n^d(nl,n‘l')/Sz  ,  (nl  i<n'l') 

n  i 

s(nl|n’l')  -  (P-5) 

in^^{nl,nl)/^z^^  (nl  -n'l') 

The  two-electron  screening  parameters  depend  on  Z«  approach¬ 
ing  CMistancy  as  z  •»  «. 
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The  preceding  discussion  can  easily  be  extended  to  com¬ 
plexes.  The  advantage  of  doing  this  is  that  the  resulting 
formula  for  the  mean  energy  of  a  complex  will  have  the  cor¬ 
rect  linear  dependence  on  Z.  On  the  other  hand,  if  one  is 
dealing  with  a  complex  that  contains  a  large  number  of  con¬ 
figurations,  the  extra  work  required  to  allow  properly  for 
all  the  interactions  may  offset  the  resultant  gain  in  accuracy. 

F.2.  The  Coefficients  for  Complex  Configurations 

Instead  of  working  with  the  energies  themselves,  one  may 
wish  to  consider  the  statistical  properties  of  the  energy  co¬ 
efficients  W.,  defined  as  the  coefficient  of  Z  in  an  expansion 
of  the  energy  in  powers  of  Z  .  In  this  subsection  we  discuss 
briefly  some  relevant  properties  of  the  W^,  before  returning, 
in  the  following  subsection,  to  the  main  stream  of  our  discus¬ 
sion.  The  considerations  of  the  subsequent  subsections  will 
apply  equally  to  the  energies  themselves  and  to  the  coeffici¬ 
ents  W^. 

The  coefficient  is  composed  of  two  parts.  The  two- 

electron  part  can  be  calculated  from  the  two-electron 

interactions  and  is  the  same  for  all  levels  belonging  to  a 

given  configuration.  The  remainder  of  W^,  which  we  denote 

by  W^,  results  partly  from  LS-dependent  interactions  within 

the  configuration  and  partly  from  interactions  connecting 

terms  belonging  to  different  configurations  in  the  same  com- 

(F~l) 

plex.  It  has  been  shown  (Layzer  )  that  in  atoms  with  a 
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partly-filled  L  shell  these  two  contributions  to  are  often 
comparable  in  magnitude  Some  recent,  as  yet  unpublished,  cal¬ 
culations  show  that  the  same  is  true  of  atoms  with  partly-filled 
M  shells  Hence  the  effects  of  configuration  mixing  within  a 
complex  cannot  in  general  be  neglected  beside  the  LS-dependent 
contributions  from  interactions  within  a  given  configuration. 

As  a  preliminary  to  calculations  of  that  allow  properly 
for  the  effects  of  configuration  mixing,  and  as  a  temporary  sub¬ 
stitute  for  such  calculations,  we  may  conveniently  evaluate  the 
contributions  to  from  interactions  within  the  configuration. 
This  is  by  no  means  a  trivial  matter  for  complex  configurations, 
since  such  configurations  may  contain  several  distinct  spectro¬ 
scopic  terms  characterized  by  given  values  of  S  and  L.  In  the 
representations  normally  used  the  energy  operator  is  diagonal 
with  respect  to  S  and  L  but  not  with  respect  to  the  quantum 
numbers  that  serve  to  distinguish  spectroscopic  terms  with  the 
same  S  and  L  Hence  the  energy  matrix,  instead  of  being  dia¬ 
gonal.  will  in  general  consist  of  multi-dmi'jnsional  blocks. 

This  IS  the  form  in  which  energy  calculations  have  usually  been 
left  in  the  past.  In  order  to  complete  the  diagonalization 
one  needs  to  know  the  ratios  of  the  .various  Slater  integrals 
that  figure  in  the  energy  matrix  elements  In  the  standard 
theory  of  complex  spectra  these  ratios  are  regarded  either  as 
parameters  to  be  determined  semi-empir ically .  or  as  theoretical 
quantities  whose  evaluation  requires  the  solution  of  a  set  of 
Hartree-Fock  equations  On  the  other  hand,  in  the  Z-dependent 
scheme  one  evaluates  the  Slater  integrals  using  straight  hydro- 
genic  wavefunctions  The  contributions  to  fzom  interactions 
within  a  given  configuration  have  definite  numerical  values 
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which  can  be  calculated  once  and  for  all,  as  can  the  elements 
of  the  transformation  matrices  that  diagonalize  the  blocks  re¬ 
ferring  to  spectroscopic  terms  of  given  S  and  L.  The  trans¬ 
formation  to  a  representation  in  which  has  no  nondiagonal 
elements  connecting  terms  in  the  same  configuration  is  the 
first  step  toward  the  ultimate  objective  of  a  representation 
in  which  the  energy  matrix  has  no  nondiagonal  elements  connec¬ 
ting  terms  belonging  to  the  same  complex . 

P.3.  The  Variance  of  the  Energy 

The  variance  of  the  energy  (in  a  given  configuration)  is 
defined  by 

var  w  =  <(W  -  W)^>  g  0^  .  (P-6) 

/p.c) 

Moszkowski  '  '  has  written  down  formulae  that  give  a  for 

w 

an  arbitrarily  complex  configuration  in  terms  of  the  o^s  for 
two-electron  configurations.  The  latter  can  easily  be  evalu¬ 
ated  by  direct  calculation.  Proofs  of  Moszkowski 's  formulae 
have  not  as  yet  appeared;  but  for  the  purposes  of  this  report 
we  shall  assume  that  the  formulae  are  correct  —  or  at  least 
that  0^  can  be  found  for  any  configuration  of  interest. 

f 

The  parameter  measures,  roughly,  the  width  of  a  con¬ 
figuration.  Knowing  o^,  one  can  estimate  (from  Boltzmann's 
formula)  the  variation  in  the  population  of  levels  over  a 
given  ccmf iguration .  More  generally,  one  can  estimate  the 
error  that  results,  in  any  particular  calculation,  from 
assigning  the  same  energy  to  all  levels  in  a  given  configu¬ 
ration.  From  one  can  also  derive  a  rough  estimate  for 
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the  average  spacing  between  levels  in  a  configuration.  Thus 
if  one  assumes  that,  say,  50%  of  the  levels  lie  in  the  range 
W  ±  0^,  one  can  reasonably  expect  the  resulting  estimate  of 
the  level  spacing  to  be  accurate  to  within  a  factor  of  2 . 

F.4.  The  Distribution  Function  of  the  Levels  in  a  Configuration 

Knowing  o  we  know  the  width  of  the  probability  distri¬ 
bution,  but  not  its  shape.  The  question  now  arises:  What 
is  the  best  assumption  as  to  the  shape  of  the  distribution 
that  is  consistent  with  our  limited  knowledge  of  the  energy 
spectrum? 

In  addition  to  W  and  o  we  know  from  direct  calculations 

w 

that  the  distribution  of  energy  levels  in  a  configuration  is 
strongly  skewed.  This  fact  has  suggested  (unpublished  Convalr 
report,  1961)  that  the  distribution  can  be  represented  by  an 
exponential  function  of  the  energy.  Once  the  exponential  form 
has  been  chosen,  the  two  parameters  on  which  such  a  distribu¬ 
tion  depends  are  uniquely  determined  by  the  values  of  W  and  o^. 

In  fact,  the  distribution  is  given  by 

N  exp  [-(w  -  W)/o^]  (W  >  W  -  o^) 

P(W)  •=  (P-7) 

0  (W  <  W  -  o^) 

Thus  the  energy  spectrum  is  supposed  to  have  a  sharp  cut-off 

a  units  below  the  mean  energy  W.  The  probability  distribution 
w 

achieves  its  maximum  value  at  the  cut-off  and  then  falls  off 
exponentially . 
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Before  this  law,  or  any  such  law,  can  be  compared  with  experiment 
we  need  to  have  some  reasonable  and  consistent  way  of  passing  from  the 
theoretical  or  experimental  energy  spectrum,  which  is  discrete,  to  a  con¬ 
tinuous  probability  distribution.  One  such  way  is  the  following.  We 
assign  to  each  level  the  appropriate  statistical  weight  (2S  +  1).  Then 
we  divide  the  energy  range  into  a  convenient  number  of  equal  intervals, 

bounded  on  the  left  by  the  lowest  energy  in  the  configuration,  and 
determine  the  total  statistical  weight  of  all  levels  falling  in  each 
energy  interval.  We  then  plot  a  histogram  of  the  running  means  of  the 
statistical  weights  against  the  energy.  This  corresponds  to  smearing 
out  the  discrete  energy  levels,  which  occur  in  bundles,  into  a  somewhat 
continuous  distribution  of  energy.  As  long  as  the  intervals  used  to 

smooth  the  data  are  reasonable  small  compared  with  the  total  energy 
range,  the  exact  manner  of  smoothing  is  of  no  importance.  One  should 
bear  in  mind  that,  as  discussed  in  Subsection  F.2,  the  distribution  of 
energy  levels  as  calculated  with  neglect  of  configuration  mixing  is  at 
best  semi -quantitative ,  so  that  much  of  the  detail  present  in  an  un¬ 
smoothed  plot  of  the  distribution  of  energy  levels  is  spurious  anyway. 

Figures  F-1  thru  F-6  illustrate  the  smoothed  distributions  of 

energy  levels  in  configurations  consisting  of  equivalent  d  and  f  elec- 

(F-8) 

trons.  Experimental  data  from  the  compilation  by  Miss  Moore'  ' ,  are 

taken  for  the  d -electron  configurations  and  theoretical  energies  are 

taken  for  the  f -electron  configuration  from  the  paper  by  Elliot,  Judd 
(F-9) 

and  Runcisian.  If  in  the  experimental  data  all  theoretical 

terms  are  not  identified,  the  missing  terms  are  given  in  each 
figure.  These  plots  exemplify  the  skewed  character  of  energy  dis¬ 
tribution,  but  the  skewness  appears  to  be  much  less  marked 
than  that  associated  with  the  exponential 
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Figure  F-1.  Smoothed  Distribution  of  Experimental  Energy 

4 

Levels  in  V  II  (Configuration  3d  ) 
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Figure  F-2.  Smoothed  Distribution  of  Experimental  Energy 
Levels  in  Cr  II  (Configuration  3d^) 
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Figure  F-3.  Saoothed  Dletrlbucion  of  Experimental  Energy 
Levels  in  Nb  II  (Configuration  4d^). 
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Level  in  the  Configuration 


distribution.  In  every  case  the  main  peak  of  the  distribution 
lies  significantly  to  the  right  of  the  lower  cut-off. 

Another  defect  of  the  exponential  distribution  consists  in 
its  not  having  the  lower  cut-off  in  the  right  place.  The  lower 
cut-off  has  in  fact  been  arbitrarily  placed  units  below  the 
mean  energy  W.  This  assignment  is  usually  in  error  by  a  quan¬ 
tity  not  much  smaller  than  o  itself.  Thus  the  main  features 
of  the  distribution  are  not  as  well  represented  by  the  expon¬ 
ential  law  as  one  would  like. 

We  have  sought  to  improve  the  statistical  description  in 
two  ways.  Firsts  we  have  added  a  third  parameter,  which  can 
in  every  case  be  calculated  exactly  with  comparatively  little 
labor.  Secondly,  we  have  chosen  the  distribution  function  in 
accordance  with  an  objective  principle  —  the  principle  of 
maximum  entropy  inference  —  which  ensures  that  the  probability 
distribution  incorporates  the  maximum  possible  uncertainty 
consistent  with  having  the  correct  mean,  variance,  and  lower 
cut-off.  These  two  improvements  are  discussed  in  detail  in  the 
following  subsections. 

F.5.  The  Lower  Cut-off 

The  lower  cut-off  is  simply  the  energy  of  the  ground 
state  of  the  configuration.  Now,  according  to  a  well-known 
rule  of  Hund,  the  lowest  spectroscopic  term  in  a  given  con¬ 
figuration  is  the  one  with  the  highest  value  of  S;  or.  If  two 
or  more  terms  have  the  same  value  of  S,  the  one  of  these  with 
the  highest  value  of  L.  One  of  the  levels  belonging  to  this 
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term,  namely,  the  level  with  M  »  S  and  M_  »  l,  is  charac- 

S  L 

terized  not  only  by  definite  values  of  M_ ,  S,  M-,  and  the  one- 

L  S 

electron  azimuthal  quantum  numbers  i  ,  but  also  by  definite 
values  of  the  one-electron  magnetic  quantum  numbers  m^.  Hence 
the  corresponding  wavef unction  is  a  Slater  determinant,  and 
the  corresponding  energy  can  be  obtained  by  a  straightforward 
calculation.  As  an  example  we  calculate  the  energy  of  the 

7 

ground  state  of  the  configuration  f  .  (See  Section  F>9) 


F.6.  Maximum  Entropy  Inference 


Consider  the  problem  of  determining  the  probability  dis¬ 
tribution  p(x)  associated  with  a  given  random  variable  X. 
Suppose  that  the  possible  values  of  X  are  known  to  lie  in  the 
range  (a,b)  and  that  the  first  k  moments  of  X  are  also  known t 


b 

J  x^p(x)dx  {i  ■  1,  . . .  ,k)  ,  (F-8) 

a 

p(x)  must  satisfy  the  further  conditions 


P(x)  >  0, 


b 

J  p(x)dx  -  1  . 


(F-9) 


a 

According  to  the  principle  of  maximum  entropy  inference 
the  best  choice  for  the  function  p(x)  is  the  one 
maximises  the  entropy 


(Jaynes 

that 


b 

H  ■  -  J  p(x)  in  p(x)dx  ,  (F-10) 

a 
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subject  to  the  k  +  1  auxiliary  conditions  expressed  by  the 
above  integral  equations.  Using  the  method  of  Lagrange  multi¬ 
pliers,  one  easily  finds  that  p(x)  must  have  the  form 


2  Ic 

p(x)  =  N  exp  [ax+ax  +  ...  +ax'|  . 

X  2  K 

The  constants  a  ...  a  ,  N,  can  be  found  from  the  k  +  1  inte- 

X  Iv 

gral  conditions. 


(F-11) 


The  entropy  of  a  probability  distribution  is  a  measure  of 

(F-7) 

the  uncertainty  inherent  in  the  distribution.  Shannon 
showed  that  no  other  function  can  consistently  bear  this  inter¬ 
pretation.  Thus  the  principle  of  maximum  entropy  inference 
may  fairly  be  described  as  the  correct  method  for  choosing  a 

probability  distribution  when  partial  information  about  the 

(F-6) 

distribution  exists.  As  Jaynes  has  pointed  out,  it  is 

» 

a  quantitative  embodiment  of  Laplace's  principle  of  insuffi¬ 
cient  reason. 


F. 7.  Application  to  the  Distribution  of  Energy  Levels  in  a  Con¬ 
figuration 

We  now  have  three  parameters  at  our  disposal  —  the  mean 
energy  W,  the  standard  deviation  of  the  energy  o^,  and  the 
energy  of  the  ground  state,  which  we  denote  by  W^.  Let 

X  =  (W  -  WQ)/e  ,  (F-12) 

where  c  is  an  as  yet  unspecified  unit  of  energy.  The  variable 
X  is  a  dimensionless  measure  of  the  energy  relative  to  the 
ground  state. 


Applying  the  principle  of  maximum  entropy  inference,  we 
find  that  the  probability  distribution  has  : he  form 

p  =  N  exp  -  b)  '^  (x  >  0)  . 

The  energy  unit  e  has  been  chosen  so  as  to  make  the  coefficient 
2 

of  X  in  the  exponential  -1.  By  integrating  the  last  equation 
over  the  range  of  x  we  obtain  a  formula  for  the  normalization 
constant  N; 

-V  -1 

N  =  27r  [  1  +  erf  (b)  ] 

where  erf(x)  is  the  error  function 

1  ^  _  2 

erf(x)  =  2tt  ^  e  ^  dy  . 

0 

Thus  p(x)  depends  on  the  single  parameter  b. 

Let 


.  -  (W  -  Wq)A  ,  a  =  a^€  ; 

u  and  0  are  the  mean  energy  and  the  standard  deviation  in 
dimensionless  units. 

The  integral  conditions  which  ensure  that  the  probability 
distribution  has  the  required  mean  and  variance  yield  two  inde¬ 
pendent  relations  fr<»n  which  b  and  e  can  be  evaluated  in  terms 
of  u  and  0.  The  following  form  of  these  raiatioQS  is  conveni¬ 
ent  for  calculation! 

•  '  • 
u  «  b  +  ip(0)  , 

(o/,i)  •  -  II  -  HP(0)  ]/2s?  . 


F-13) 


F-14) 

F-15) 


F-16ab) 


F-17) 
( F-18) 
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By  substituting  the  first  equation  into  the  second,  we  obtain  a  trans¬ 
cendental  equation  for  b  in  terms  of  the  ratio  (a/^) .  The  numerical 
solution  of  this  equation  for  a  range  of  values  of  a/^  can  be  obtained 
easily.  Table  F-1  gives  or/u  for  various  values  of  b,  and  tlM  sans  is 
plotted  in  Figure  F-7.  Having  obtained  the  value  of  b  from  this 
table,  one  can  evaluate  the  energy  unit  e  from  equations  (F-17)  and 
(F-16a). 

The  agreement  between  the  predicted  and  observed  distributions 
seems  to  be  satisfactory.  If  future  studies  should  disclose  a 
systematic  discrepancy  between  the  predicted  and  observed  distribu¬ 
tion,  then  the  very  existence  of  such  discrepancies  would  represent 
fresh  information  about  the  probability  distribution  and  could  be 
incorporated  into  the  conditions  of  the  problem.  Horeover,  it  is 
likely  that  any  such  systematic  discrepancy  would  have  an  assignable 
physical  cause.  Thus  the  principle  of  maximum  entropy  Inference 
appears  to  be  a  powerful  method  for  Isolating  those  properties  of  com¬ 
plex  spectra  that  are  susceptible  to  interpretation  in  reasonably 
siiiq>le  terms. 

tht  99iblLniU9ff  Sitttrf 

Consider  a  transition  array,  defined  as  Che  aggregate  of  lines 
resulting  from  all  possible  transitions  connecting  two 


TABLE  F-1 


^  and  b 


b 

S. 

Jti-. 

0  0 

.75122 

C.l 

.70515 

0.25 

.63495 

0.5 

.59937 

0.75 

.bini 

1.0 

.53512 

1.25 

.48832 

1.5 

.43851 

1.75 

.39107 

2.0 

.34742 

2.5 

.28235 

3.0 

.23567 

3.5 

.20203 

4.0 

.17677 

f 


subject  to  the  k  +  1  auxiliary  conditions  expressed  by  the 
above  integral  equations.  Using  the  method  of  Lagrange  multi¬ 
pliers,  one  easily  finds  that  p(x)  must  have  the  form 

2 

p(x)  =  N  exp  [a^x  +  a^x  +  ...  +  1  • 

The  constants  a.  . . .  a,  ,  N,  can  be  found  from  the  k  +  1  inte- 
1  k 

gral  conditions . 

The  entropy  of  a  probability  distribution  is  a  measure  of 

(F-7) 

the  uncertainty  in^ierent  in  the  distribution.  Shannon 
showed  that  no  other  function  can  consistently  bear  this  inter¬ 
pretation.  Thus  the  principle  of  maximum  entropy  inference 
may  fairly  be  described  as  the  correct  method  for  choosing  a 
probability  distribution  when  partial  information  about  the 
distribution  exists.  As  Jaynes  ^  has  pointed  out,  it  is 

a  quantitative  embodiment  of  Laplace's  principle  of  insuffi-- 
dent  reason. 

F. 7.  Application  to  the  Distribution  of  Energy  Levels  in  a  Con¬ 
figuration 

We  now  have  three  parameters  at  our  disposal  —  the  mean 
energy  W,  the  standard  deviation  of  the  energy  o^,  and  the 
energy  of  the  ground  state,  which  we  denote  by  W^.  Let 

X  »  (W  -  W^)/e  , 

where  e  is  an  as  yet  unspecified  unit  of  energy.  The  variable 

* 

X  is  a  dimensionless  measure  of  the  energy  relative  to  the 
ground  state. 


(F-11) 


(F-12) 
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given  configurations.  What  can  we  say  about  the  probability 
distribution  of  line  frequencies  in  such  an  array,  given  the 
probability  distributions  for  the  energy  levels  of  the  two 
configurations? 


Let  W  .  denote  the  energy  of  the  ith  level  in  the  configu- 
ration  a.  The  frequency  corresponding  to  a  transition  ai  bj 
is  given  by  Bohr's  frequency  condition, 

hv(ai,bj)  =  W  .  W-  .  . 

*  ■'  ai  b] 

We  assume  for  the  sake  of  definiteness  that  configuration  a 
lies  above  configuration  b,  so  that  the  frequencies  defined 
by  this  equation  are  always  positive. 

If  departures  from  LS  coupling  are  not  too  great,  the 
combination  spectrum  [the  spectrum  of  frequencies  v(ai,bj)] 
will  have  a  clumpy  structure.  Transitions  between  levels 
belonging  to  the  same  pair  of  spectroscopic  terms  give  rise 
to  multiplets.  and  in  LS  coupling  "the  separation  bet%ifeen 
multiplets  is  large  compared  with  the  separation  between 
lines  in  a  multiplet. 

The  distribution  of  lines  in  a  multiplet  and  the  dis-. 
tribution  of  multiplets  in  a  transition  array  present  distint 
problems.  In  a  first  approximation  we  may  assume  that  the 
two  distributions  are  stochastically  independent  •—  that  is. 


(F-19) 


251 


that  the  form  of  the  probability  distribution  for  lines  In  a 
multlplet  does  not  depend  on  the  position  of  the  multlplet  in 
its  transition  array.  Similar  remarks  apply  to  the  case  of 
pure  jj  coupling. 

On  the  other  handj  if  the  coupling  scheme  is  not  close  to 
either  LS  or  jj  coupling,  neither  the  Individual  energy  spectra 
nor  the  combination  spectrum  will  exhibit  any  marked  fine  struc¬ 
ture.  In  this  respect  the  case  of  intermediate . coupling  is 
easier  to  deal  with  than  either  extreme  LS  coupling  or  extreme 
jj  coupling. 

There  is  a  second  important  respect  in  which  the  stochastic 
description  of  the  combination  spectrum  is  simpler  for  inter¬ 
mediate  coupling  than  for  LS  coupling.  In  LS  coupling,  transi¬ 
tions  connecting  levels  belonging  to  terms  of  different  multi¬ 
plicity  (different  S)  are  strictly  forbidden.  This  well-known 
selection'  rule  substantially  reduces  the  number  of  lines  in  a 
transition  array  that  occur  with  nonvanishing  intensity.  More¬ 
over,  it  forces  one  to  deal  separately  with  the  various  classes 
of  transitions  corresponding  to  the  various  possible  values  of  S. 
In  jj  coupling  transitions  connecting  terms  differing  by  more 
than  one  quantum  number  j  are  forbidden.  In  intermediate  coup¬ 
ling,  on  the  other  hand,  intersystem  transitions  are  not  only 
allowed,  but  often  have  high  probability,  since  intermediate 
coupling  prevails  in  heavy  atMtis,  we  may  reasonably  assume  that 
the  choice  of  an  initial  level  and  the  choice  of  a  final  level 
are  ato<diastloally  independent  events.  The  remainder  of  this 
discussicm  proceeds  from  this  assumption. 


T 


If  we  regard  the  quantitiea  v,  W  .  W  as  random  variables 

A  D 

[whose  possible  values  are  the  numbers  v(al,bj) ^bj' 
spectlvely]^  then  the  aaaun^tlon  just  stated  means  that  the 
random  variables  W  and  W.  are  stochastically  Independent^  so 

A  D 

that  the  random  variable  v  Is  the  difference  between  two  stoch¬ 
astically  Independent  random  variables,  nils  fact  enables  us 
to  write  down  the  probability  distribution  of  v  In  terms  of  the 
probability  distributions  of  and  According  to  a  well- 

known  theorem  In  probability  theory,  the  first  distribution  Is 
just  the  convolution  of  the  probability  distributions  of  W  and 
Let  Pj^Cx)  denote  the  probability  distribution  of  and  p^ 
denote  the  probability  distribution  of  Then  the  probability 

distribution  p^^  of 


W 


ab 


W  -  W, 


(P  -20) 


la  given  by 


00 

—00 

From  this  equation,  or,  more  directly,  from  the  definition  of 
stochastic  Independence,  one  can  show  that  the  expectation  and 
the  variance  of  W  .  are  given  by 


W 


a 


(f-22) 


Var  W 


ab 


Var 


+  Var  W.  . 
D 


(F-23) 


As  an  Illustration  of  these  rules  we  consider  the  example 
of  two  normal  distributions.  By  an  appropriate  choice  of  sero- 
points  we  can  always  ensure  that  ■  0,  so  that  we 
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may  put 


p  (x)  *  (2Traf)“^  exp  [~x^/2a^]  , 

S3  S 

p^(x)  =  (2vo^)~^  exp  [-xVza^]  . 

Since  the  normal  probability  distribution  is  syimnetric^  the 
probability  distribution  Pj^  of  coincides  with  p^  in  the 
present  example;  but  this  is  not  true  in  general.  A  short 
calculation,  in  which  one  performs  the  integration  by  ex¬ 
piating  the  square  in  the  exponent,  gives 

1 

p^^(x)  =  exp  l-x‘/2c^^l  . 


ab 


2^2 
a  +  o. 

a  b 


the  probability  distribution  of  the  difference  is  also 
normal  and  i€s  variance  is  the  sum  of  the  variances  of  and 
w^.  ihis  is  the  well-known  "reproductive  property"  of  the 
normal  distribution. 


We  assume  that  the  probability  distributions  p^,  have 
the  form  derived  earlier  in  the  present  report.  It  is  conveni¬ 
ent  to  introduce  the  function 

1  (X  >  0) 

€(x)  « 

0  (X  <  0) 

With  the  help  of  this  function  we  can  write  the  probability 

distribution  for  W  in  the  form 

s 

P^(x)  «  9(x;a^,P^)€(x)  , 


(F-24) 

(F-25) 


(F-26) 

(F-27) 


(P-28) 

(F-29) 
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where 


! 


<p(x;o  6  )  =  N  exp  [-(x  -  ^  )^/2a\]  .  ( 

cl  o  ^  S  A 

The  definition  of  is  exactly  analogous.  Note  that 

Pj^(x)  -  Pj^("X)  ( 

We  therefore  obtain  for  the  probability  distribution  of  W^^ 


00 

Pj^jj(y)  »  J  <P(x;o^,3^)<p(x  -  yj0j^,Pj^)6(x)€(x  -  y)dx  .  ( 

—00 

From  the  definition  of  the  function  €(x)  it  follows  at  once,  that 

e(x  -  y)  (y  >  0) 


Hence 


€(x)€(x  -  y) 


€(x)  (y  <  0) 


Pjjjj(y)  -  J  <P(x;o^,p^)<p(x  -  y;a^,Pjj)dx  (y  <  0)  ;  ( 


‘P(x;o^,P^)<p(x  -  y;Ojj,Pjj)dx 


y^9(y  +  d*  (y  >  0) 


The  integrations  are  straightforward  (ms  the  Macheaaticsl  Mote  at 
the  end  of  this  appendix).  The  final  result  is: 


F-30) 


-31) 


F-32) 


F-33) 


F-34a) 


F-34b) 
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1  -  erf 


/|3  ^  +  p.  O’  ^  +  y  O’ 

[  a  *>>  fe-1— 

W  o.  \2  W  2  +  0-  i) 

'  a  0  a  D 


)] 


-  P.,(y) 


\Vb  V2  ((r^2  +  a^2)  ^ 


(  y  <  0) 
(F-35a) 

2  vT 

(y  >  0) 
(P-35b) 


where  P^^(y)  .  ^ 


1/2 

K  ■  Pb  -  >-/  1 

2  ".''b 

r.  *''"b 

and  erf(x)  has  been  defined  in  e'quation  (F‘15) 

F.9.  Ground  State  Bnarav  for  a  Configuration  of  Kouivalent  Electrons 

Hs  consider  a  configuration  of  equivalent  electrons  The  energy 
of  the  cosftletely  filled  shells  is  conaon  to  all  the  teras  of  the  con* 
figuration*  It  can  be  calculated*  if  desired,  but  is  of  no  conse¬ 

quence  for  our  purpose.  FOrtheraore,  the  average  energy  of  such  a 
configuration  can  be  calculated  In  a  straightforward  manner  (Slater^' 


(i^)  (f-36) 

2i 

and  (i.i)  (oi,  nX)  (F-37) 

k 

where  the  coefficients  f  for  configurations  of  p*  d*  f  and  g  electrons 

|r 

are  given  In  Table  J*1  and  F  are  the  Slater  P  Integrals  defined  In 
Equations  (B-26)  and(B-27). 


2S& 


Us  shall  describe  below  a  method  of  calculating  the  energy  of  the 
ground  term. 

It  Is  well  known  that  a  subshell  i  can  contain  at  most  2(2i  +  1)  electrons 
and  Is  therefore  exactly  half  filled  when  q  ■  2i  +  1.  It  is  also  well 
known  that  the  term  energies  (with  respect  to  the  average  energy)  of  a 
configuration  with  n  electrons  are  the  same  as  those  of  a  configuration 
with  [2(2f  +  1)  -  q]  electrons.  It  Is,  therefore,  necessary  only  to 
consider  configurations  which  are  half  or  less  than  half  filled,  l.e. 

<  2i  +  1. 

The  ground  term  of  any  configuration  has  the  largest  value  of  S 
for  that  configuration.  If  several  terms  have  the  same  largest  value  of 
S,  say  It  Is  the  one  with  the  largest  value  of  L,  say 

among  them.  Each  term  Is  highly  degenerate,  but  one  of  the  SL 
states  belonging  to  the  term  la  such  that 

^  OMUC  wmX 

“s- 

Ud  “l  • 

TIm  sero*order  state,  which  beloi^s  uniquely  to  the  above  set  of  values 
of  the  quantum  nwriiers  can  be  found  easily.  Mg  >  8^^  arises  idien  all 
spins  are  parallel  and  In  this  case  of  all  n  electrons  Is  +1  (.'i:  -z) 
Therefore, 

and  tha  Mjltlpllcity  of  tha  ground  tom  la  (2  S  +  1)  or  (q  +  1). 

■MX 
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Tlie  largest  value  of  consistent  with  the  above  requirement  on  the 

m  and  M.  is  obtained  by  a  simple  application  of  the  exclusion  principle, 
s  S 

The  m.  values  of  the  q  equivalent  electrons  must  be 
£ 

=  i,  f  -  l,  . .  (i  -  q  f  1)  (F-39) 

For  this  case 

M,  -  L  -  f  (2f  -  q  H  1)  (P-40) 

L  max  i 

The  corresponding  zero  order  state  is  designated  in  the  usual  way  by 
(f)'  (i-1)^  .  (/  -  q  +  1)' 

where  the  quantity  in  parenthesis  is  the  value  of  m^  of  individual 

electrons  and  the  +  sign  indicates  that  m^  -■  »■  1/2.  Tlie  ground  terms  of 

a  few  configuration  of  type  are  given  in  Table  P-2.  For  a  shell  wliich 

is  exactly  half  filled  (q  -  2i  >  1),  it  is  readily  seen  from  Equation 

(F-5)  that  L  =0  and  the  term  is  ^ 
max 

TABLE  F-2 

Ground  Terms  for  Configurations  of  Equivalent  Electrons 


Configuration 

Multiplicity 

1=  q  +  1] 

L 

1-  1  (2i-q+l)J 

Ground  Term 

2 

P 

3 

1 

3 

P 

d' 

4 

3 

S’ 

f' 

8 

0 

«S 

8^ 

6 

10 

258 


I" 


The  energy  of  the  ground  term  in  the  absence  of  external  fields  is 

the  same  as  that  of  any  one  of  the  zero  order  states  of  this  term,  say 

of  the  zero  order  state  mentioned  above.  This  energy  can  be  found  by 

(P“3) 

standard  methods  (Condon  and  Shortley  '),  and  is  given  by, 


■EZ[ 

a  k  '• 


k  .  -  a 

a  (i  m 


i*  ^ 


b^  (i  m^,  i  (ni,  ni) 

(P-41) 


where  the  first  auimaatlon  la  over  all  two  electrons  combinations  and 

the  second  susmatlon  Is  over  k  which  takes  on  all  even  values  from  zero 

k  k  *  6 

to  2i.  The  coefficients  a  nnd  b  are  defined  by  Equations  8  (14)  and 

6  '  (F*3) 

8  (6)  of  Condon  and  Shortley^  and  are  tabulated  by  them  in  Tables 

6  6 

1  and  2  for  s,  p,  d  and  f  electrons. 


As  an  example  we  calculate  below  the  energy  of  the  ground  term 
of  configuration  f^.  The  zero  order  state  of  lowest  energy  is 


3"^  2*  I*  O'^ 


g 

and  uniquely  belongs  to  the  S  term.  Table  P-3  gives  the  appropriate 
k  k 

a  and  b  for  f  electrons.  From  the  sums  of  these  quentities,  given  in 

the  last  row  of  this  table,  it  readily  follows  chat  the  energy  of  the 

8  7 

ground  term,  8,  of  the  configuration  f  is 
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TABLE  F-3 


*n<i  ^  £ -electrons 

Mote:  a?  «  1  for  ell  esses 


"  Q  “  “i  ^  ■? 


Electrons 

a 

2 

s 

X  225 

X  1089 

6 

a 

X  7361.64 

;  b2 

1 

lx  225 

X  1089 

b^ 

X  7361.64 

3* 

2^ 

0 

-21 

-6 

j 

25 

30 

7 

1^ 

-15 

3 

15 

10 

54 

28 

0* 

I 

lO 

o 

IS 

-20 

0 

63 

84 

-1 

-15 

3 

15 

0 

42 

210 

-2^ 

0 

-21  * 

-6 

0 

0 

462 

-3^ 

25 

9 

1 

0 

0 

924 

2* 

1^ 

0 

-7 

-90 

15 

32 

105 

0^ 

0 

-42 

120 

20 

3 

224 

-1^ 

0 

-7 

-90 

0 

14 

378 

-2* 

0 

49 

36 

0 

70 

504 

-3^ 

0 

-21 

-6 

0 

0 

462 

1* 

0* 

12 

6 

-300 

2 

15 

350 

-1* 

9 

1 

225 

24 

40 

420 

-2* 

0 

-7 

-90 

0 

14 

378 

-3^ 

-15 

3 

15 

0 

42 

210 

0* 

+ 

1 

12 

6 

-300 

2 

IS 

350 

-2^ 

-42 

120 

20 

3 

224 

-3* 

18 

o 

1 

0 

63 

84 

-1* 

-2^ 

-7 

•90 

15 

32 

105 

-3^ 

-15 

3 

15 

10 

54 

28 

-2* 

f 

-3* 

0 

-21 

25 

30 

7 

2 

•42 

-77 

Bl 

168 

616 

5544 

F.IO.  Mathawatlcal  Votm 


From  Equations  (F*30)  and  F*34a)  we  obtain,  for  y  <  0, 

(x  -  (x  -  y  -  ** 


00 

’•‘.‘•if  “p 


2  CT. 


2  a. 


dx  <F-43a) 


00 

J  exp  |-{8K^  +  2qx  +  r)j  dx 


(P-43b) 


where 


2  2 
(T  *■  +  a.  ‘ 

^ _ b _ 

2  2 

2  a/  ^  0^2 


-P,  ^b^  ‘^a 

<1  2  2 

2  \ 


r  ■ 


2  2 
2  cr  ^  a.  ^ 

a  b 


(P-44a) 


(P>44b) 


(P-44c) 


The  integral  on  the  right  hand  aide  of  Equation  (P*43b)  can  bn  evaluated, 
by  means  of  formula  (4b),  Table  314,  page  65  of  the  ceepllation  of 
definite  integrals  by  GrSboer  and  Hofreiter^^'^^^  Applying  this  formula 
we  obtain  Equation  (P-35a). 


Por  y  >  0,  «M  combine  Equations  (P*30)  and  (P’*34b)  and  get 


00 

",  "b  /  «p 


(b  -  (-y  ^  g.)  (Z  -  P|,)^ 


2  CT. 


2  <7. 


(lx 


(P-45) 


which  Is  the  same  as  Equation  (P-42a)  If  i^e  Interchange  the  subscripts 
a  and  b  and  replace  y  by  (-y).  Performing  the  integration  as  before 
we  obtain  Equation  (F-3Sb). 


\ 
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APPENDIX  G 

H.  E.  Stubbs  and  G.  A.  Victor 

THE  EFFECT  OF  LINE  WIDTH,  STRENGTH,  AND  SHAPE 
ON  RADIATIVE  TRANSFER  IN  AN  EXPANDING  GAS 

G.l  Introduction 

The  details  of  line  width,  shape,  and  strenQth  have  a  very  important 

role  in  the  calculation  of  radiative  transfer  phenomena  in  non-expanding 

gas  clouds.  This  dependence  on  details  of  line  shape  introduces  very 

{G-1  2) 

serious  complications  in  most  astrophysicsl  applications'  *  since 
reliable  information  on  line  shapes  is  available  only  for  simple  gases 
as  hydrogen  end  helium The  object  of  this  appendix  is  to  investi¬ 
gate  a  simple  ouxiel  of  a  rapidly  expanding  gas  which  results  in  expressions 
for  line  absorption  which  are  largely  independent  of  detail  line  shape. 

The  origin  of  the  simolification  in  this  model  is  found  in  the  fact  that 
the  Doppler  shift  due  to  the  velocity  of  expansion  becosms  the  dead.annt 
factor,  and  that  the  details  of  line  shape  due  to  collision  and  thermal 
Doppler  broadening  become  unimportant. 

Tliis  appendix  will  be  divided  into  three  sections.  Part  G.2  will 
be  devoted  to  soaw  simple  estiawtions  of  line  widths  due  to  the  therauil 
Doppler  broadening,  and  to  compare  it  with  the  Doppler  shift  ^e  to  the 
expansion  velocity.  Part  G.3  will  show  that  resonance  absorption  is 
independent  of  line  shape.  Part  G.4  will  present  a  simple  calculation 
for  an  ei^nding  hydrt^en  gas  for  the  transmission  of  the  82,  260  em~^ 
line. 


M4 


G.2  Doppler  Shift  Due  to  V«loclty  of  Exp<nilon 

We  shell  construct  the  following  ideellsed  model  for  our  expending 
cloud: 


The  nueiber  density  of  hydrogen  etoms  is  given  hy: 

p  (r,t)  ■  1.5  X  10^/t^  perticles/cm^  r  < 

p  (r.t)  -  0  r  >  R^(t) 

R«(t)  -  V^t 

o  o 

V  -  1  X  10*^  c  (6-1) 

o 

t  is  the  time  in  seconds 

10  . 

c  is  the  speed  of  light  ■  3  x  10  cm/sec. 


the  reletive  velocity  between  two  perticles  is  elong  their  line  of 
centers  end  is  given  by: 


(6-2) 


Where  r^2  perticles  seperetion: 

The  reletivistic  Dozier  formule  is: 


X  •  X_J  cos  6  (6-3) 

o 

Since  the  velocity  is  elong  the  line  of  centers  we  cen  set  cos  dal. 
Since  «  C  we  heve  the  reletion  >  .01,  therefore: 

Sua  - 

The  fmmtle  for  the  thenMl  Doppler  width  is 
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ax_  -  X 
T 


(G-5) 


U2 


for  Hydrogen  we  can  construct  the  following  table : 


3/2  kT 

5X_ 

T 

100  ev 

3.2  X  10*\ 

1.6  X  10*^ 

10  ev 

io‘\^ 

5.0  X  lO'^ 

1  ev 

3.2  X 

1.6  X  lO"^ 

.1  ev 

io’\ 

•4 

5  X  10  ^ 

t 

o 

From  Table  6*1  we  see  that  the  Doppler  width  due  to  theraal  aotion  is 
always  sHich  ssaller  than  the  sMxiwiai  Doppler  shift  due  to  the  velocity 
of  expansion.  The  results  of  the  rest  of  this  appendix  depend  on  this 
fact. 

G.3  '  Independence  of  Resonance  Absorption  on  Line  Shene 

Let  a  photon  of  frequency  w  be  enitted  et  r  ■  0.  Let  the 
probability  of  emissimi  of  a  photon  of  energy  between  liu  end  IT  (w  ^  du) 
per  unit  tisa  per  particle  be  giyen  by  P(wX  The  probability  that  a  given 
atom  will  absorb  a  photon  of  energy  Hu  is  given  by, 

p(w)  ■  k  P(w)  (C-6) 

from  the  principle  of  detail  balance. 
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The  probability  that  a  photon  of  frequency  u  will  tranaverae  a  dis¬ 


tance  r  before  being  absorbed  is  given  by: 

r 

A(|,w)dJ 


I 


1  (u.r)  ■  e 


(G-7) 


with 


A(6,w)  -  p(t)  p  (u*[r,u^]) 


(G-8) 


where  p(t)  is  the  density  of  particles  which  can  abaorb  this  resonance 
radiation.  I(u,r)  can  therefore  be  written  as: 


p(l)  p  (<*>[1“ 
l(u,r)  ■  e  o 


where  ve  have  used  the  Doppler  fonaula: 


(G-9) 


ct 


(C-IO) 


and  evaluated  this  at  a  particular  tisM  t^. 


We  intergrate  this  over  all  frequencies,  weighted  with  P(u)  to  get  the 
probability  that  any  photon  will  go  a  distance  r. 


T(r) 


■/. 


du  P(u)  e 

00 


-jf  P 


which  will  be  re-written  as 


’I 


♦00 

I(r)  -  j  d«  P(w) 
•00 


.k  /  d|o(|)  P((j(l.  ^1) 


(G-11) 


(0-12) 
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If  p(|)  It  •  constant  out  to  a  distance  r,  which  is  consistent  with  our 
cloud  nodal,  we  have: 


where:  K' 

Finally 


Using  the 


if  wo  let 

then  dr)  ■ 
so: 


r"" 

T(r)  ■  I  du  P((it)  e 
•'-00 


b)*Au(r) 
r(i)dt 


-  kp  Aa)(r)  •  ^ 

o 


U(l- 


r  *L  K*p(u')du* 

T(r)  ml  P(w)du  e 
-00 

result  of  part  G.2  we  can  calculate  this  for  r  ■  oo 

•00 


■t-oo 


T(od)  ■  I  P(w)dwe 
•00 


F(w*)dw' 


^00 


1]  ml  P(a)')du' 


F(u)du 


we  see  that  for  an  infinite  cloud 


2M 


(G-13) 


(C-14) 


(C-15) 


(G-16) 


Iron  Equation  (6-16) 


resonance  absorption  is  Independent  of  line  shape.  For  a  finite  cloud, 
it  is  only  the  photons  thatsre  eaitted  near  the  edge  that  will  show  resonance 
absorption  that  is  dependent  on  the  form  of  P(u) .  The  ratio  of  the  thick* 
ness  of  this  edge  to  the  cloud  radius  is  essentially  Ve  there* 

fere  conclude  that  resonance  absorption  is  controlled  priasurily  by  line 
strength,  for  this  simple  model  of  an  expanding  cloud.  The  object  of 
this  appendix  is  to  investigate  a  siiqple  model  of  a  rapidly  expanding  gas 
\diich  results  in  expressions  for  line  absorption  which  are  largely 
independent  of  detail  line  shape. 

Arguments  similar  to  thoae  .given  above  for  resonance  absorption  can 
be  given  for  line  absorption  due  to  two  close,  but  non*overlapping  lines. 

It  is  easy  to  see  that  an  absorption  line  centered  at  can  only  absorb 
radiation  from  lines  centered  at  such  that: 

It  is  again  obvious  that  this  absorption  is  independent  of  line 
shape  except  for  lines  where  the  equality  holds  in  Equation  (G*17).  flio 
effect  of  a  finite  cloud  on  the  absorption  of  non*overlapping  line  radia-  . 
tion  does  little  but  introduce  some  simple  geometrical  factors. 

G.4  Transmission  of  Hvdroaen  82.260  cm*^  Line 

The  line  of  Interest  Is  the  S2,260  cm*^  line  which  Is  nsltted  by  a 
transition  from  the  first  excited  state  to  the  ground  state.  We  will 
asaume  that  the  gas  has  a  temperature  history  as  shoun  In  Figure  6*1. 

Xn  addition,  we  assume  that  the  gas  la  at  thermal  equilibrium  at  all  times. 


TEMPERATURE  (•¥) 


OIO*40~I5SO 


TIME  ($0CJ 
» 

Figure  G-1.  Teaperaturc  During  Expansion. 
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Ac  equilibrium,  the  population  of  etetea  is  given  by  the  Sahe  equation 
“n  “  exp 

■  number  deneiCy  of  protcma 

■  number  density  of  electrons 

■  -onization  potential  (  r  degree) 
a  ■  principle  quantum  number 

a  nuiri)er  density  of  atoms  in  state  n 

'he  number  density  of  particles  which  can  absorb  the  82,260  cm'^  line  is 
given  by 

P(t,T)  -  (t)f(t,T)  (C-19) 

N,(T) 

where  f(t,T)  »  TTTT 

/ 

Istiaiations  of  f(t,T)  from  Equat ion <C.X8) lead  to  the  results  shown  in  Che 
following  table. 


SI 

1 

j  2  [cm] 

*t 

r 

00 

10‘® 

1.5  X  10*^  i  2  X  10^ 

_ ! _ 

7  X  lo’f 

5.7  X  10® 

2.5  X  10'* 

10'* 

1.5  X  10^* 

'  3 

i  2  X  lO"* 

1 

( 

7  X  10*f 

QQHQIII 

1.3  X  lO'^ 

10'^ 

1.5  X  10^* 

<  A 

1  2  X  10^ 

7  X  10*f 

.26 

1.4  X  10'® 

10'* 

1.5  X  10^* 

2  X  10* 

7  X  lO^f 

1 

1.4  X  10'* 

1.5  X  10® 

2  X  10® 

.X 

7  X  10 
■ . . 

1 
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From  Equation  (G>16)  tra  hava: 


where  the  subscript  Indicates  time  dependence.  is  given 
ing  eiq^ression: 

j..  .  P  .  P  (tjT)  A  t 

*■13  3 

8jt  V  V  8«  V 

o 

&  *1 

A  is  the  Sinstein  coefficient  5  x  10°  sec*  .  The  results 
at  various  times  is  shown  in  Table  G*2.  FiguresG>2  end  G-3 
variation  of  T(ao)  with  time  for  this  simple  siodel. 


(G-20) 

by  the  follow* 

(G-21) 

for  K* 
show  the 
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TRJkNSmiTTED  TO  OOHT) 


0I0A40-ISI00 


to'*  to'*  to'*  to'* 


TtmtE  (seconds) 

Figure  G-2.  TrensaiMion  in  Expending  Mndiua. 
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0IDA40- 13130 


TIME  <9ec) 


Figure  G*3.  Fractitm  of  Photoae  oscaping 
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